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SYNOPSIS 


Electricity, vhich is a vital input to almost all 
the sectors of the economy, provides an important basic 
infrastructure for economic developnent of a nation. The 
future developments of a country rely heavily on the rate of 
growth of power generation and their efficient disposal. A 
challenging situation is emerging out of the present in- 
adequacies of supply and fast increasing power demands. To 
avoid serious repercussions on the economy due to shortage 
of power supply it is necessary that careful planning in the 
power sector is undertaken. Over all planning in the power 
sector involves forecasting of futiire demands for electricity 
and the development of optimal strategies for expansion of 
electricity generation capacities to meet the future demands , 

In this thesis an attempt is made to develop quan- 
titative models for forecasting of electricity demand and 
expainsion of electricity generation capacities to meet the 
ititure demands. The thesis comprises of two parts. Part I 
deals in. til the development of econometric and univariate 
stochastic time series models for forecasting the future 

energy and peali power demands of electricity for India 

* 

till 2000 A,D, Statistical control limits of forecasts 
for specified levels of significance have been estimated. 
Econometric models have also been developed for emalysing 



electricity demand in the industrial sector. The indus- 
tries considered are : Iron and Steel, Hon-ferrous metals, 
Chemicals, Vehicles, Textiles, Paper, lining and Quarrying, 
Engineering and Food iiidustries. 

In part II of the thesis an analytical approach 
for planning for capacity expansion is presented. Planning 
for capacity expansion involves .solving of two sub-prohlems 
in an integrated fashion, to yield minimum cost capacity 
expansion programme. The t\;o sub -problems which pertain to 
the determination of capacities and operational strategies 
are solved in an iterative manner using Bender's decompo- 
sition technique to yield an overall minimum cost plan for 
capacity. Capacity expansion sub-problem is formulated as 
an integer programming model. Demand is assumed to be 
concentrated at load centres between which capacity and 
energy is to be transmitted. The size, location and type 
of i^ants are determined to meet a projected demand schedule. 
The operational planning sub-problems which finds minimum 
cost operation of the system is formulated as a non-linear 
programming problem with linear constraints . The feedback 
between tlie capacity expansion sub-problem and the opera- 
tional j^anning sub-problem is repeated until the expansion 
sequence stops improving. The methodology has been applied 
to a case study and the numerical results have been 
presented. 



CHAPTER I 


IKTRODUCTIOi: 

Throug}iout human Mstorj’’, the foundations of 
civilisation have rested heavily on supplies of various 
fonris of energy. The nineteenth ceiitury industrial 
revolution -was highly energj'' depandent, as is the proces 
of industrialisation nov; unden-ra-y in India. Energy’’ has 
played a decisive role in the economic development of a 
nation. To tlie extent of its availability' it has stimu- 
lated or hindered economic groirth. 

Energy is of urdversal use, being not only a 
component of productive process but also an element 
fundamental to i-relfare. Productivity is directly influ- 
enced by the amount of energy which can be incorporated 
into the productive process . Variations among the 
countries in their levels of consumption of energy is 
attributed to their degree of industrialisation, product 
mix in industry, relative importance of energy'’ intensive 
processes in industrial structure, efficiency of utilise 
tion, climatic factors and costs . 

The ^lationships underlying economic growth 
and energy consui25)tion vrfL’tH reference to ’the Indian 
context has been exhaustively discussed in the Energy 



2 


Survey Committee report of Government of India (1) . 

Planning in India is set to accomplish the oTj,1ectives of 
removal of econor.ic bach-wardness and attainment of self 
reliance. The s’agyested strategies for achievins these 
goals heavily relj^ upon the usage of energy. Davar(2) 
has emphasised the distinct role of energy’- as an effec- 
tive catalyst and driving force for economic development. 

The exponential gro^/tb rates pf consumption of 
the exhaustible and non-renei^?able resource like energy, 
the four-fold increase in price of crude oil in the 
world market have compelled the world community to focus 
their attention on energy. Shortage of energy lias 
created economic and political problems of enormous 
magnitude. Prompt and effective solution of these prob- 
lems have become imperative to permit development to 
continue and to avoid threat to national stability. 

1 .1 Role of Electricity 

Electricity is an advanced form of energy and 
it constitutes a major component of total energy require- 
ments, Development of electricity has made possible to 
talie advantages of renewable i-zater resources. Indus- 
trial development, automation, mechanisation and urban 
progress have virtually become a function of electricity 
supply. 
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In the Indian context, electricity provides 
an important and basic infrastructure for development. 

It is a vital input to industrj^, agriculture, transport 
and is of particular importance to the developing rural 
sector vliich needs more and more pouter for its agricul- 
tural operations, small scale and agro-industries. The 
futirre development of the country depends to a great 
extent on the rate of groirtli of po\;er generation and their 
efficient disposal. A challenging situation is emerging 
out of the present inadequacies of supply and fast 
increasing power demands. This calls for a careful 
plsuining in the po\;er sector. Reliable estimates of 
future demands of electricity is a first step in this 
process . 

1,2 Heed For Lonr Term Forecasting 
Of Electricity 

Electricity is accepted as a basic commodity 
and hence planning for electricity production should 
precede or at least closely be in harmony i;ith planning 
in other sectors. The creation of production facilities 
for electricity involves long gestation periods. 
Especially this is so in developing countries v;ith in- 
adequate facilities for manufacturing of power plant 
equipment and scarcity of finances. 



If 

Thc investment in powr sector is huge. This 
calls for great care and attention in nalzing important 
decision'- in pov/cr sector cspec:?all" ’/ith limited dores- 
tic and foreign capital resources. Fc'>" developing 
countries tliis necessiates projection of their future 
demands for electricity on a long term, basis and then 
plan for the production facilities, k knoidedge of the 
magnitude of demand obtained through appropriate fore- 
casting helps decisions on planning indigoneous manufac- 
ture of eauipnent and complimentary’' facilities. For a 
perspective plan, demand forecasting for a time horison 
of ti/entyfive years (spanning five five-year plans in 
India) uould be referred to as long term projection for 
purposes of this study, 

1 *3 Forecasting Electricity Demand For 
Specific Povrer Systems 

A higlily aggregated demand for electricity for 
an entire countiy can be obtained by projecting the 
demands for individual power systems existing in a region 
or grid. Individual utility systems have to devise better 
ways to anticipate aggregate and coincident demand of all 
customers for planning and successful operation. 

A utility S 3 rsten needs forecasts for a variety 
of purposes. The lead time and frequency of forecasts 



5 


depend on their uGage. For day to day operation hourly 
forecasts for at least a day ahead is needed. Devising 
of naintenance policies reoioire i/eckly forecasts uith a 
lead time of six to t^^enty four months . Generation 
plannj.ng requires monthly forecasts u-pto sixty months in 
advsnce and long term corporate planning require annual 
forecasts for a load time of fifteen to tventy years. 

In general, for an elect ricitj’- supply undertal.ing involv- 
ing both generation and distribution it caii be stated 
that the secret of rendering better service at economical 
prices and reasonable profits to the organisation is to 
predict needs and plan ahead, 

1 ,4 Comnonents Of Electricity Demand Forecasts 

Power system load forecasters are basically 
concerned with the forecasting of total electric energy 
requirements and pea]-: power requirements , Energy fore- 
casts provide a basis to estimate future revenue. Anti- 
cipated peali power demands detemine company's investment 
in additional generation and transLiission facility to 
ensure adequate supply. 

A record of the weekly, monthly or annual peal: 
and energy demands indicate the basic components of 
demand as well as the salient features of load growth. 

The deiiand pattern normally has the following components ; 
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1, A continual grov/t]T of demand o^'er tine: This grov;th 
is aLtributed to an underlying component called base 
demand (trend caiiponent). Eventiiough short varia- 
tions exist, long terii grov^th is approximated by a 
smooth varj-ing trend curve. 

2. Variations in cleraand repeating themselves \/ith a 

period: These variations constitute the cyclic 

or seasonal component whose pattern may be preserved 
over the years but magnitude may gro^j ’.mLth time by 
random amounts , 

3« Random Variations : Superimposed on tlie above tvro 

components there are random ccsiiponents, with magni- 
tudes much smaller than seasonal components. This 
is designated as noise component. 

Pluctuations in the above three components 
depend on customer class, seasonal changes and economic 
activities . 

1*5 Important Techno-Economic Characteristics Of 
An Electricity Production-Distribution System 

An electric power production - distribution 
system has the following techno-economic characteristics, 

1. Electricity is a non-storable item. Hence instan- 
taneous power demand has to be met by instantaneous 
response. 
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. Failure to meet dewand may entail purchasing pov/er 
at higiicr costs from a neighbouring region or systen 
or cause serious conseoucnces to custcuers, 

3 . He serve caoacity is expensive a.nd 'unproductive . 

hence a pover svsteL: tries to keen a reserve as loi: 
as possible consistent uitlx a satisfactory level of 
reliability of service. 

Integration of po^'er systems bet\7een nei^.b curing 
states or regions for economic incercliange of energy 
and capacity is a ranidly emerging phenomenon. 

5 . It talies a gestation period of four to five years 
to plan and install generation and transriission 
requirerients . Future corporate and systems develop- 
laent requires long term forecasts! ifteen to tventy 
years in adirance. 

1 .6 Plaimina For Expansion Of 
Electric Power Systems 

With increasing automation and industrialisa- 
tion it is evident tliat demand for an advanced form of 
energy like electricity wo'old increase. V/liile reliable 
estimates of future demand is an important step in 
planning, the responsibility of power system planners 
does not end here. Power sj^stem planning can be defined 
as a rational program for development of an electric 





pov/er system so that it can evolve in an orderly and 
economic manner. It includes the rationalisation of 
standards of service, anticipation of trends in equinmen 
design and coordination of the various elements into a 
\/ell designed whole . It is particularly concerned \rith 
plans for changes anci additions to generation and 
transmission facilities . Surir'.arily planning of an elec- 
tric pov/er system has to solve the problems of v;hen and 
what facilities should be provided where to assure ade- 
quate service, 

1 ,7 Need For Long-range Planning Of 
Power Systems 

Systems planning is necessary because as the 
system grows existing capacity is insufficient to meet 
future demands. In absence of advanced plenning quality?- 
of service v:ill suffer due to long gestation periods of 
capacity additions . Fartiier the onornous amount of 
investment in power develoment, inflation of costs and 
increasing cariying charges demand a careful planning in 
this sector as the costs of an incorrect plan is very 
high. 
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1 .8 The Ifaturc and Complexities Of 
The Taslv Of Flanninfr 

Modern povTer systoris conFift r cou.plox of 
installations vitli iiitric:te coniiecticnj . The concept 
of fominy a grid has to come u? in currenL practice. 

The process of integrctio’i effecting c:cc]i.. n^c of capscit;' 
and power entails research into these systons as a vliolc , 

The factors of emit sizes, pattern of load 
groirth, fuel cost charges technical end econociic compe- 
titiveness of alternative plants, levels of reliability. 
Spinning reserves, transmission cost and inflation have 
strong implications for planning. 

The relative economics of po't7er resources in- 
fluences planning for expansion. Hydro potentials, 
fossil fuels and nuclear energy constitute the primary 
sources of power* Hydel resources are cheaper but limited 
in quantity and located at few places. Fossil fuel 
reserves are vast but thermoelectric power is economical 
near the source of fuel. Costs of nuclear plants arc 
compara]Dle but'its development is restricted due to avai- 
lability of nuclear fuel and radiation hazards. 

Technology for other forms of energy resources are far 
from a stage of economical exploitation. These consi- 
derations of relative economics suggest an optiEial mix 
of hjrdro, theimal and nuclear plants for economic power 
planning in a region or grid. 
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For future expansion prof'raninies decesions should 
be s\ich that they fit into an overall coordinated plan 
that balances sizes against reserve requirements, location 
ageinst transrlssion requirements and dates of installa- 
tion against risk of loss of load. 

Considering the factors stated above it is 
evident that determination of the optimum variants for the 
structural development of an electric power system is a 
complex and m\£Lti variate task. It is also a voluminous 
task on account of present day growth of interconnected 
systems involving large number of transmission lines and 
plants. The decisions pertaining to the location of 
plants, their type and capacity are too complex to be 
handled by conventional techniques , Further a planner is 
interested in the sensitivity of optimal plans to varietj^ 
of conditions and changes in parameters . It is therefore 
vitalljr important to form a theory of optimum development 
of an electric power system and devise procedures for 
optimisation throu^ algorithms. The determination of 
these optimal variants should take into account prevailing 
economic, technical, natural and social conditions. 

The solution of this complex and intricate 
problem can be arrived at by the application of mathema- 
tical programming and simulation techniques with the aid 
of modern computers , Tiie solutions vriLll provide decisions 
regarding the location and time phasing of plants, their 
typ)es and capacities. 
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1 .9 The Characteristics Of A Pov:er Systei.i 

anc" Tiielr Influence On Ilathenmtical flodels 

The main features of an electric pov;er system are 
characterised hy the follo\ang: 

1. Protabilistic nature of source data : Tlris comprises 

forecasts of demand, d^’anges i 2 i technical parameters, 
performance of equipment and future trends of economic 
policy. 

.2. Discrete nature of capacity'- additions : The capacity 

expansion taJces place in discrete steps according to 
unit capacities chosen- for plants, the standards in 
force and availability of machinery and equipment in 
market • 

3. Non-linearities of system interaction: The economic 

and technological link that exists betvreen individual 
components of a pover system are higlrly non-linear in 
nature . 

The aforementioned characteristics suggest that 
the mathematical model of the power system \dll be large 
scale, stochastic, non-linear, discrete and dynamic. 
Unfortunately mathanatical theory at the present stage 
of development cannot provide solution to tliis large and 
complex problem. Hence mathematical models to be used 
in practice have to be considerably simplified and sys- 
tematised by using the following strategies. 
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(a) Linearising of principal eqviations expressing 
technical and econoirdc relationshiiDS . 

(b) Disregarding discrete nature of systens develop- 
ment in linear ajid non-linear models. 

(c) Eliminating the stochastic dimensions by taking 
e^rpected values of parameters and variables . 

(d) Selecting basic strategies prior to solution. 

(e) Determining an optimui: zone of solution and 
ma]:ing sensitivity studies to consider effects 
of probable variations in parameters . 

Hie problem of planning can be considered from 
two points of view. 

(a) Optimisation of the system as a whole. 

(b) Optimisation of individual components like 
generation, transmission etc. 

The first approach involves the development of 
global model \diich is difficult to tacld.e due to problem 
of dimensionality and limitations of present day computers . 
The second approach, although, permits the determination 
of optimal or near optimal solutions for individual 
components of the system, the problem that underlies this 
approach is the difficulty in deteiroination of complex 
non-linear functional relationships linking individual 
coE5)onents. However an attempt can be made to optimise 
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the development of a power sj’-steni ty an "order of 
procedure" v;Iiich corresponds as closely as possible to 
the real process of development. 

1.10 Ilotlvation and Scope For 
The Present Study 

In previous sections we have emphasised the role 
of electricity in the economic development of a nation and 
the consequent need and importance of forecasting future 
requirements and long range planning in this vital sector 
of the economy. It needs to be pointed out that develop- 
ment planning is a continuous and sequential process 
involving mobilisation and efficient use of resources. In 
tliis context the important aspect of development planning 
in the electric energy sector is the formulation of 
suitable policies to carryout economic activities of gene- 
ration, transmission and distribution of electric energy 
over time and space. There emerges a need to study the 
problems arising in connection \7ith rational utilisation 
of investment in the power sector in view of heavy canital 
investments involved in the development of power systems, 

A survey of the literature reveals that very 
little work has been reported in the area of forecasting 
the future demands of electricity and planning. for 
expansion of capacities of poi/er systems in the Indian 
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context. Further, there exists very little evidence of 
the relevant applications of the sophisticated tools of 
ocono’netrics , stocliastic tir’e -series analysis and tlieoilcs 
of prediction and operational research in the field of 
fore casting and. planjiirn of power systems in India. These 
above nentionod concideratio.ns 5 have provided us a j.iotiva- 
tion for the present study. 

This dissertation concerns itself ’.n.th the deve- 
lopment of methodoloijies for forecasting electricity demand 
and for planning tlie optimal expansion of cap'.citios ox' 
power systems, T \*70 methodologies viz. - the econometric 
methods and stociiastic time series analysis and prediction 
notheds have been proposed for forecasting electricity- 
demand incorpors-ting tlie probabilistic characteristics of 
electricity demand. Integer and non-linear mathematical 
programming models to be used in an iterative nianixr for 
solving the co.pacity expansion sub-problem and operational 
planning sub-problem respectively have been proposed for 
solving the overall capacity expansion problems of an 
electric power system. The models proposed have been 
applied to case study of an actual povror sys'^iem. 

The present dissertation comprises of two parte. 
Part I of the study deals vrith the problem of forecasting 
future dejaands of electricity. Part- II of the disserta- 
tion is devoted to the problem of planning for capacity 



expansion in an integrated poi/er system. Chapters II to VI 
constitutes Part I and Chapters VII to IX constitute Part II 

Chapter II presents a brief revic of the litera- 
ture on forecasting the demands for electricity. The surve 
is intended to provide an overview of relevant existing 
work on statistical econometric and stochastic techniques 
of analysis and projection of der.and. 

Chapter III discusses the various methodologies 
of forecasting electricity demand and presents ,a relative 
evaluation of these methodologies . Reasons for selecting 
the methodologies used in this dissertation are also 
discussed. 

Chapter IV is devoted to econometric r’‘=>thods of 
analysis and forecasting of electricity demuid. Cimg’ -- 
macro nod^els for peak po\irer and energ" demand have bejn 
formulated and estimated. On the basis of the estimated 
equations point and inter\^ forecasts of annual peal: and 
energy demand have been obtained for India for a period of 
25 years, i.e,, till 2000 A.D, Tliis chapter also presents 
a few econometric models for analysing industrial elec- 
tricity demand. The industries considered are : Food, 
Chemicals, Textiles, VeMcles, Engineering, JtLning and 
Quarrying, Paper, Iron and Steel and Hon-fcrro'':s metals. 
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Chapter V comprises of a brief description of 
the theory and netnodolory of univariate stochastic time 
series anclysis and their forcco.stinr . Ti’ 'e-scrics rnodclvS 
have been fitted to d?tp on past electricity' d.ei-.anc’ , Fore- 
casts based on these models liave Leen ol Gained for India 
for r, tine snan of t'vX'ntyfive years. 

Chapter Mi", is devoted to discussion of the 
results obtained by econometric and stochastic time series 
methodolohies forecasting, electricity demand. Results 
for analysis of industrial electricity dejiand are also 
presented. Scopes for Airthar research into the area of 
electricity dera.nd forecasting have been suggested. 

In Chapter VII \k present a brief survey of the 
literatuire in the field of planning for canocity expansion 
%ath special emphasis on capacity expansion planning of 
electric po^/er systems. Mathematical prograiiming and other 
models of po-ier system planning have been reviewed. 

In Chapter VIII v;e present an approach to find 
the ninimum cost capacity expansion policies involving 
the determination of size, location, type of plants as 
well as imports and. exports and long distance transmission 
lines between demand centres of an electric power system, 

A methodology for solving the capacity expansion problem 
and opera. tional planning aspect of the problem in an 



17 


iterative manner is presented, Tlie overall capacity 
planning rjroMeri is dccoriposee into two sul.-problens by 
Bender’s decomposition principle. The capacity expansion 
sub-pro 'blc;:': is fon'iLlatcd cS an intej_oi' nrogrc-ming model 
and t}ie operational pla-'ining problem is for/rnlated as a 
non-linea.r prograsr’ing model , 

The nemierical results obtained by the applica- 
tion of the met’aod clo^.y presented in Chapter VIII to a 
case study are presented in Chapter Ih. Conclusions are 
drav/n on the basis of results obtained, and scone for 
further vork is also presented. 



P A R T ■» I 


FORECASTING THE DEMAND 
FOR ELECTRICITY 



CHAPTER II 


F0EEGA3TIRG ELECTRICITY DEITvlID - A 
LITERATURE I^E\T[EW 

This chapter is devoted to a brief revie\; of 
the literature in the field of forecastin:-, ^/ith particular 
emphasis on forecasting of electricity demand. For 
convenience the literature reviev is divided into three 
parts. Section 2.1 presents the relevant literature in the 
field of forecasting by statistical, and econometric methods 
and their applications to forecasting the demand for elec- 
tricity. Section 2.2 is devoted to the survey of the 
pertinent work in the area of stochastic theories of 
prediction and their applications to electricity demand 
forecasting. Section 2,3 d.eals with the review of other 
general method of projection as applied to the forecasting 
of electricity requirements , 

2»’1 Review Of Literature On Statistical 

and Econometric Methods Of Forecasting 
The Demand For Sliectricitv 

In forecasting techniques time-series and 
multiple Ingres Sion analysis plaj!" a very important role. 

One popular time-series model is exponential smoothing. 
Exponential smoothing is based on a weighted average of 
two sources of evidence, one the latest (most recent) 
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observation and the other is a v.-iltio computed one period 
before. As such it is an easy and quid: nethod since very 
little information is needed for obtaining, forecasts . 

Reference is made to Frown (3, 4, 5> 6 , 7 ) and Broun 
and Ileyer ( 8 ) ^dio have presented an extensive description 
of the theory, and application cf this technique to various 
problems of forecasting. Tlie exponential smoothing tech- 
nique aH-ong with its various applications have also been 
discussed by Winters ( 9 , 10) liith.(ll) Geoffrion 2egels(13) 

V/iener (1^) , Duffin and WlrLddin ( 15 ) , Duff in and 

Sclmidt (16) , Harrison (17^18), Harris on and Davies ( 19 ) 

Harrison and Scott ( 20 ) , Kirby ( 21 ) , Shiskin (22) , 

Holt ( 23 , 2 ^) , Cox ( 25 ) > Buff a ( 26 ) , Buff a and 

Taubert (27) , Arrow et. al. ( 28 ) , Vftiitin ( 29 ) , 

V/elch ( 30 ) , 1‘Iagee and Boodman (31) , 1-Ioore ( 32 ) , 

Greene (33) , Eilon (34) , Buchan and Koenigsberg ( 3 ^) 

and a host of others. 

Theil and Wage (36) have formulated a stochastic 
model underlying the procedure of adaptive forecasting of 
an economic time-series. For a case in \;hich the time 
scries to be forecast has no seasonal components, they 
determined the wei^ts to be used in adaptive forecasts 
Tidiich are optimal in the sense of minimum mean square 
error. The model presented by Theil and Wage is a good 
example of the exponential smootliing technique. The 
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model is expressed as 

+ Residual ... (2.1) 

“ ^t-1 €t (2.2) 

where at time t 

3!^ = Trend value 

= Trend change 

= Seasonal component. 

Exponential smootliing procedure predicts the trend value. 

At the end of the t period observation X^ is at the 
forecaster's disposal. From equation (2.1) it is equal 
to X ^ apart from S^ and the residual , 

In the absence of any Information on residual and 
S^, they are rejiLaced by their expected values, \diich are 
zero and respectively, where I is the length of 

the seasonal cycle. Now is the new evidence 

on trend level. The latest trend value X^„^ refers to 
period t - 1, X^ is obtained by adding to and 

^ is replaced by • ^Hais loads to the following 

exponential smoothing procedure for Xj.. 

^t = ^t - ^t-I >+<■>-<)( et-1 ^ 

( 2 . 3 ) 

Given X^ can be determined in terms of most recent 
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observation and previous computed values X^^l t 

t-i and 

Exponential smoothing provides procedures for 
detecting and adjusting to changes in forecast series 
rather than predicting these changes. These methods by 
and large do not predict in the behavioral sense. There 
is little attempt at explanation of causality and no 
information beyond the historical data of time series is 
used, 

Herlove and Wage (37) demonstrated that adaptive 
forecasting are optimal in a much wider sense, Ahtougli 
the series generated by the Thcil and V/age model is non- 
stationary there exists a simple transformation of the 
series which converts it into a stationary series , This 
observation permitted Nerlove and Wage to apply the 
Wiener - Hopf (1^) theory for stationary time-series 
to Idle "transformed series. Further contributions to the 
technique of adaptive forecasting was made by Chow (38) 
Dudman (39) , Wheelright and llabridabis (^0) ^ Jain and 

Batra (^1) , Packer (^2) , Trugg and Leach (43) y Box 

and Jenkins (44, 4^) , Mcclain (46) , Ifcclain and 

Thomas (4?) , Morris and Glassey (48) , Buff a (49) and 

Griffin (50) • Ibst of exponential smoothing time- 

series models including the one given in (2*1) and (2.2) 
were designed to break a tinK series into its components 
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namely trend, cyclicity and noise. This aD.lo’ircd the 
forecaster to gain insight into the past history of a 
scries through study of changes in the individual compo- 
nents , 

I^is cussing monblily average load forecasting on 
the "Tennesse Valle;/ Authority” system Uev (51) emphasised 
the h re aiding of a time -series into its components and use 
of seasonal indices for forecasting seasonal components, 

A seasonal index (obtained as a ratio of load for a given 
month to that of trend value for the same month) was 
utilised by I'Ie\;, The seasonal component of demand was 
obtained by extrapolating monthly seasonal indices and then 
multiplying those index of m.onth by corresponding projected 
value of trend component. The computation resulted in 
forecast of averoge montlily load. The trend component \ms 
projected by extrapolating an appropriate time function 
fitted to load data by least scuares technique. 

Doobie (52) has proposed a method for the 
computation of coefficients even when the fitting function 
consists of a finite number of sine and cosine teims , 
Without breaking the tim.o series into their components, 
suitable trigonometric functions in additions to time 
polynomial functions \/ere incorporated to include sea- 
sonality, Cliristianse (53) used the general exponential 
smoothing model for short term load-forecasting of an 
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electric power system. Hourly loads for a power system 
were obtained by Gupta and Yamada (5^) by adaptive short 
term forecasting. Weather information was utilised for 
determination of those forecasts. 

Berry and V/hiting (55) recognised the fact that 
trend curves must level at some stage and hence attempted 
to fit a logistics trend to the time scries. Inspite of 
all possible refinements in forecasting by trend extrapo- 
lation, tliis methodology was critised by Cowden (56) . 

Coi/don gave three reasons for his criticism: (1) Difficul- 
ties in finding a logical basis sufficient to justify the 
type of trend selected. (2) Thou^itless extrapolation can 
produce ridiculous results. (3) Economic and social 
conditions responsible for trend may not continue to be 
applicable outside the time -span of data and forecasting of 
trend component is subject to statistical errors which 
tends to increase as vre move awa 3 ^ from the centre of the 
time period to which trend curve is fitted. 

The second category’- of forecasting techniques, 
multiple regression analysis, predicts a change in the 
forecast series throu^ explanatory variables for a given 
time-series . The explanatory variables are selected on 
the basis of economic theory and forecaster's judgement. 
However unlike time-series models multiple regression 
models usually do not utilise infoimation contained in 



historical pattern of a time series. A regression model 
for forecasting electricity demand was used by Ilieneman 
et, al . (57) who recognised the responsiveness of system 
loads especially to seasonal factors. Chen and Winters 
(5^8) developed a model \f.ilch forecasted daily peak load 
of an electric power system tj combining exponential 
smoothing and multiple linear regression models, leading 
to a h 5 ?-brid forecasting model. 

Yaraada (59) extended Chen and VJinter's model 
to the forecasting problem of a more general class by 
replacing the linear regression part of the model by other 
regression models. His general model is of the foim 

Y^ = 

S - ^ \t * t 

n 

®t = ^ <=1 * ‘I -O 

St = It ‘ \ ^it > * '1 h 

VJhere 

Y^ = Total demand at time t 
:= Base demand at time t 
R^ = Residual at time t 
S. = Seasonal demand at time t 


25 


Xit, ^ 2 t’ '*• ^nt explanatory variables 

0 ^ is a zero neaiijconstant variance independent noise 

c(j are snoothing constants, 

L = Length of seasonal cycle, 

^t -1 ^t-L expected values of and esti- 

mated at period t -1 and t-L respectively, Yamada also 
presented a metiiod for identifying the parameters involved 
in tlie forecasting model. 

The literature reviewed above mostly dealt vrith 
statistical and regression analj’-sis techniques of forecast- 
ing, and their applications to electricity demand forecast- 
ing. In the following paragraphs we present a survey of 
the literature on econometric analysis and forecasting of 
electricity demand. 

The general econometric techniques of forecasting 
has been discussed by various authors and researchers. 

Some of the important contributions are due to Johnston 
(60) , ICLein (61) , Goldberger (62) , Theil (63,6lf,65) 

Dhrsrmes ( 66 ) , Christ (67) , Tinbergen ( 68 ) , 

Tintner ( 69 ) , Wannacot ( 70 ) , Leser(7l) , Cramer ! 

( 72 ) Hood and Koopmans ( 73 ) , Koonmans ( 74 ) and 
Wold ( 75 ) • Iwo good general references for forecasting 

economic time series are the books by Buttle r and Platt 
( 76 ) and Chisholm (77) • Ibe statistical aspects of 
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estima-ting econometric models have been discussed by 
Maulinvad (78) . Suits (79) has demonstrated the use of 
ail actual econometric model as a tool of forecasting for 
the U.S. economy. The utility of various available econo- 
metric models as instruments of forecasting v;ere evaluated 
by Stekler (80) , 

Enormous amount of rescarcli has been done on the 
development of input - output models , and their applica- 
tions to forecasting. Lconlief's (8l, 32, 83) 
pioneering research into this area has opened up many 
avenues for use of this methodology. Modern contributors 
to the development of this technique include Morgens tern 

(8^-) , Rasmussen (85) , Stone (86) , Dorfman (87) , 

* 

Dorfman, Samuels on and Solov (88) , Koopman ( 89 ) , 
Chenery and Clark ( 90 ) , and a host of other authors and 
researchers. The input - output tables for India have been 
discussed by Mathur and Bhardwaj ( 91 ) * 

The dLassical theories of demand for commodities, 
(both capital and consumer), have been presented by Wold 
and Jureen (92) . All the above mentioned references 
on econometric methodologies (including input-output 
analysis) provide us techniques for projecting the demand 
in future for the various sections of the econony inclu- 
ding electricity. 
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As far as the specific applications of econometric 
methods to analysis and forecasting the demand for elec- 
tricity are concerned, it is observed that most of the 
researchers have focussed their attention on sectoral 
demands and especially on the domestic (residential) sector. 
Very few studies exist for the cormercial and industrial 
sectors, Ta^dLor (93) provides an evaluation and critique 
of the few studies in the area of econometric analysis of 
demand for electricity. Taylor attributes many of the 
problems in modelling the demand for electricity to the 
existence of multi-step block pricing, the fact that demand 
for electricity is a derived demand and the existence of 
distinct short-run and long-run elasticities for each 
class of consumer. 

Classical theory of consumer demand sees the 
consumer as maximising a iitility function defined over all 
goods subject to his level of income. However, while 
recent years have seen rapidly increasing uses of theore- 
tically plausible demand functions (Refer Parks (9**-) > 

EOuthakker and Taylor (95) , Phlips ( 96 ) , Taylor and 

Weiserbs (97) , Brown and Ileien ( 98 ) , Christensen 

et al* (99) ) there does not exist a single econometric 

stu^ of the demand for electricity for which this is the 
case* Demand for electricity has usually been approached 
in isolation or else in conjmetion with the demand for 
its close substitutes. 
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Ilouthakker (100) has discussed the econometric 
implications of the existence of a price schedule. The 
literature has focussed very narrov;ly on the question of 
type of price - marginal or average that should be inclu- 
ded in the demand function. Eventhough the tJieoretical 
implications of quantity discounts and block tariffs have 
been stated by Buchnan (101) , Gabor (102,. 103) , and 

Oi (1C4) , tlie implications of the price sched-ule in the 

case of electricity for the equilibrium of the consumer and 
therefore for the demand function itself has not been 
systematical! 3 ’’ investigated, 

Houthakker’s (105) study focussed on residential 
electricity demand in the U.K, Using cross-section as 
^rell as time-series observations, he estimated linear and 
log-linear regression models, by including average raoney 
income per household, marginal price of electricity, 
marginal price of gas and average holding of domestic 
electric equipments as the explanatory variables. His 
results shoved an income elasticity of 1.17, price elas- 
ticit 3 ’- of -0,89, and cross -elasticity v/ith marginal price 
of gas as 0.21. Houthakker did not clarify as to whether 
these elasticities referred to*short-ran or long-run 
demands for electricity. 

The standard and the most ambitious reference 
on demand for electricity in the U.S. is the monograph 
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of Fisher and Kajsen (29) • Tlic authors anal 3 ’’sed residen-* 

tial and industrial demand and distinguished explicitl:^ 
tetucen short-run and long-run demands . According to the 
aut'iors the short-run demaiid was a function of iitilisation 
rate of e^jLsting equipment vhile demand in the long-run 
v;as influenced hy a choice of the size of the capital 
stock, Houthakker and Taylor (95) estimated an equation 
for personal consumption expenditure on electricity that 
was based on a state adjustment model of consximption. The 
short and long-run elasticities of Income were 0,13 and 
1,93» while that of price were -0,13 and -1.89 respectively, 

Bauxter and Rees (107) have developed models 
explicitly for the industrial demand for electricity. 

Various types of fuels along\d.th capital and labour were 
used as an input to the production functions . The demand 
function for electricity was derived from a Cobb - Douglas 
type production function. The main conclusions from their 
analj'^sis was that relative price changes are not unambi- 
gyously an important determinant of growth in electricity. 
The chief determinants were growth in output and changes 
in technology. There existed a marked responsiveness of 
demand to relative fuel prices in some indiis tries, while 
in others the price elasticity of denand was zero. 
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Wilson (108) analysed the residential demand for 
electricity. Tlie exogenous variables incorporated vere 
price of electricity, average price of gas, median family 
income, average number of rooms per household and climatic 
variables. The results of particular interest were subs- 
tantial negative price elasticity and negative income 
elasticity. Ilis results vis-a-vis price elasticity of 
demand are in conflict v;ith those obtained by Fisher and 
Kaysen (106) , The peal: load Tiienomenon has received great 
deal of attention by various authors such as Boiteux (IO9) 
Euchnan (101) , Gabor (102) y Houtlialcker (109) and Lewis 
(110) . Inspito of this the econometric. literature on 

pealc load demand appears to consist of a single study by 
Cargil and Meyer (111) . The regressors for estimating 

the demand function were the ratio of the average revenue 
per KWH to average price per therm of gas, real per capita 
income, employment of production workers in manufacturing 
and time. Their equation explained ^0% of the variation 
in monthly hourly demand and indicated that price increase 
has a negative effect on demand, v/here as income is of 
little consequence. 

In a study that in concerned more with develop- 
ment of methodology than obtaining results, Anderson 
(112) analysed producers demand for energy in the U.S, 
primary metals industry. Anderson* s analysis, was based 
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on tlie methodology of Fisher and kaysen, but uith the 
following important differences : 

(1) Focus was on demand for total energy rather than only 
electric power. 

(2) Allowance v/as made for quantity discount of inputs, 

(3) Price of input and effects of competing or related 
input prices were incorporated, 

(If) Effects of variation in industries were considered. 

The dependent variable vfas KWH of electricity purcliased per 
unit of value added. Explanatory variables included v;ere 
price of coal, price of electricity and oil, average wage 
rate of production workers in primary metals. Hie price 
elasticity of demand was negative, substantial, and highly 
significant. 

Mount, Chapman and Tyrel (113) analysed both the 
short-run and long-run demand for electricity for tliree 
classes of consumers, residential, commercial and industrial. 
The models were estimated by pooled cross-section and 
time -series data of annual observations , Lagged dependent 
variable was used as a predictor along -vatli population, 
income per capita, average price of electricity, lagged 
price of gas and appliances. Models for all the tliree 
sectors was estimated using ordinary least square (OLS) 
tedanique as ijell as ths instnmiental variable procedure. 
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The long-^^in elasticities demonstrated that electricity 
demand is generally price elastic in all sectors hut income 
inelastic,- Population exliihited approxir,r.tely unit 
elasticity. 

In another stud}' Anderson (1l4) analysed the 
residential demand for electricity. He estimated two 
different models - one for predicting the stocks of energy 
using equipment, and the other for predicting energ}' con- 
sumption. The second model, involving double -log form of 
equation, considered income, price of various sources of 
energy and several demograpliic ■'^ariahles . 

The demand for electric po\;er for tliree major 
classes of consumer, residential, commercial and industrial 
was analysed by Lyman (115) » Biis study contained some 
innovations including the use of firm data as opposed to 
state aggregates and use of non-linear demand functions. 

The non-linear demand function was in line \7ith those 
developed by Box and Cox (116) , Zarembka (II 7 ) , and 
Zellner (II 8 ) • Models were estimated using maximum 
likelihood methods . Using the variable transformation 
functional form Lyman suggested a linear semilogarthmic 
model for residential sector and double logarthmic form 
for ccamnercial and industrial sector. Demands were found 
to be price elastic in each consumer class and residential 
demand was found to have a positive correlation with 
income • 
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Ilouthakkcr, Veerleger and Sheehan (11*9) employed 
a logartlimic flo\/ adjustnent model of the type used by 
HouthsJdier and Taylor (95) in snaD.y£?ing the demand for 
electricity in the residential sector. Time-series and 
cross-scction sainuLes of sta.tc aggregates vere used. The 
model \ 7 as estimated by using the error component technique 
developed by Bales tra and iJerlove (120) , Llerlove 

(121) Ciccheti and Smith (122) examined the implication 
of selecting alternative price measures for statistical 
properties of the estimated demand relationship. The 
criteria for selecting the best price measure was defined 
in form of Ramsey’s tests (123,12^,125) • The findings for 
the residential sector suggested that with appropriate 
ad;3us 'jiient for simultaneity, the average price measure is 
preferable to a measirre based on typical electricity bills • 

Ashbury (126) examined the residential market 
with three cross-sections of fortyeiglit states, using 
average revenue as a price measure. Both OLS, and 2S1S 
were used. Results iidicated that price elasticity esti- 
mates were stable vath introduction of income, substitute 
prices, density and climate as explanatory variables. 

Halvore sen’s (128, to 131) econometric models 
were designed to analyse the demand for electricity in 
the residential sector of U.S. as well as whole of U.S. 



The estimated models indicated that the long run direct 
elasticity of demaaid vith regard to price was at least 
uniterj’-. 

Hav/kins’ (131) study dealt vith the demand for 
electricity in the residential, coi.imercial and industrial 
sector of Ne^; Southwales and Australian territory. Cross- 
section data vas used. There was little evidence that 
commercial and industrial sector demand was responsive to 
price. The commercial sector demand function was inter- 
preted as a production function with labour and services 
of electricity as factors . The production function showed 
increasing returns to scale. 

The other important studies worth mentioning 
on analysis of electricity demand are due to Wilson (132) 
Tyrel3(133) , Chapman, Mount and Tyrell (134-). 


2.2 Demand Forecasting Using Stochastic 
TTme-^eries Analysis and Theories 
Of Prediction 

The application of probability theory to 
forecasting is not new. Soe^ of the important contri- 
butions in the field are due to Hannan (13?) , Nerlove 
and Wage (37) , Theil et al., (36), Weiner Wold(136 

“^37) , Whittle (138) , Box and Jenkins (139 ,1^0, 1^1 j 1^2) 
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Bartlett (14-3) , Lailoy (144) , Doob (14^) , Parzen (146) 

Anderson (14?) , Ivakhenko (148) , Qiierouille (149) , 

Eartliolomew (150) , Rosenblatt (151) , IColraogrove (152) 

Pugachev (153) , Borde and Shannon (15^) , Gabor (155) ) 

Lubcecl: (156) , Cramner and Leadbctter (157) , Cox and 

ItLller (15^) } Zadeh and Ragazzini (155) , and numerous 

other researchers who developed the idea of stochastic 

prediction of a time series. Unfortunately this fast 

developing science of stochastic theory was never applied 

to forecast electricity demand. This v;as natural because 

except for a few research papers by Box (139), Nerlove (37)j 

Theil (36) and Wliittle (138) the stochastic piediction 

theory was almost universally applied to stationary time 

series. Modern stochastic approach to problem of filtering 

and prediction was launched by Wiener (l4) in his now 

classical work. He developed techniques for synthesis 

of optimal linear systems for filtering and prediction of 

stationary time series, V/iener showed that the linear 

filter was the absolute optimum filter for Gaussian noise 

under minimum mean square error criterion. Wold (136) 

carried out the idea of linear filter forward and stated' 

that any stationary time series could be considered as 

white 

an output of a linear filter ■vd.th/noise as input. Given 
a time series a linear filter given by its impilse res- 
pose function could be derived, which has its output the 
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desired forecasts, V7hittlc(136) used the concept of • Z’ 
transform in order to estimate the impulse response func- 
tion of the prediction filter based on Wold's theory. Use 
of Z transform simplified mathematical computation and made 
prediction scheme very elegant. Carrying fonmrd the idea 
of stationary time series ho formulated a prediction scheme 
for those types of non-statio2iai^^ scries \iiLch could be 
reduced to a stationary time series by finite linear trans- 
formation. I'ox and Uenkins (44) introduced a practical 
prediction scheme for handling non-stationaiy time-series 
and showed it to be optimal for a particular class of 
stochastic process. While most literature dealt \d.th the 
problem in the time domain, an occasional try vra-S made to 
define a Power spectrum of a non-stationary time series, 
Priestley (160) developed an approa-ch for spectral analysis 
of non-stationary time-series , which v;as based on the 
concept of evolutionary spectra i.c,, a spectral function 
\irhich was time dependent and has a physical interpretation 
as local energy distribution over frequency. Ii^pite of 
Priestley’s elegant mathematical analysis, the method was 
to divide the series into segments, compute spectra, con- 
sider each segment as stationary and then extrapolate the 
evolving value of spectra. In a paper Rao and Shapiro 
(161) proposed an alternative scheme for using the 
evolutionary spectra in adaptive smoothing. Ihe smoothing 
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constants vero dctorniincd as a function of niaximurn change 
in various frequency' conponents of successive spectra. The 
loca.tion and type of change indicated the disturbance in 
the underlying stochastic process generating the series. 
Literature cited above has teen mostly theoretical but 
Theil and Wage (36) formulated a stochastic model under- 
lying the procedure of ada'etivc forecasting of an economic 
time-series. Ilerlovc and Wag'3 suggested that given the. same 
underlying model adaptive forecasting was optimal in a 
nrudi \d.der sense. Gouts et. al.(l62) carried forward the 
general idea of Hcrlove and Wage for predicting those tjrpe 
of non-stationary series which can be reduced to a stationary 
series by finite linear transformation. Farmer (I 63 ) has 
presented the application of the theory non-stationary 
time-series prediction to electricity demand estimation. 

The method utilised an adaptive time-series approach repre- 
senting the load cycle in terms of characteristic function. 

New 16^ )■- appears to be one of the first inves- 
tigators to add probabilitj’- dimension to the forecasting 
.of electricity demand. He decomposed the time series into 
its components viz, cyclic, trend and noise. Forecasts 
and the associated probability distribution for each compo- 
nent -were obtained separately. The I-bnte -Carlo gaming 
technique proposed by Hammersley and Handscombe (165) 
was then used to canbine the probability distribution of 
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components to get the final forecast. He claimed that 
the use of the Konte -Carlo gaming technique offers certain 
advantages. Tlie advantages accrue due to the fact that 
the teclinique does not rcoiiire complex convolution scheme 
to combine composite distributions, Lathom et, al. , (166) 
have suggested a procedure for the integration of various 
probabilistic conponeiit forecasts into the total foreca.st. 
The integration vas obtained through a subjective distri- 
bution method. This procedure is not verir sound in the 
sense that different forecasters v;ould. obtain different 
forecast distributions from the same data. 

Stanton and Gupta et. al(l67,l68,l69»170)have 
applied stochastic time series models to long range and 
short range forecasting of demand for electricity. 

In a dissertation Gupta (171) has applied the . 
stochastic theory of prediction to forecast the weekly, 
montlily and annual pealc demands in an utility system. 

Some state estimation type of modelling for load 
forecasting were used by Toyoda and associates (172,173, 17^ j ' 
175>17^f 177)»Jiie sequential filter techrdque suggested by 
Schweeppe and Wildes (178)and Larson et. al.(179) were 
used by Toyoda et# al. for this short terra forecasting. 

The sequential filtering techniques for forecasting were 
proposed by Kalman and Bucy (180) , Sorensen (l8l) , and 
Aoki'(l82) , Sage and Husa, and Mehra (183,184,185). 
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2.3 other Kethoos Of Forecasting-; 

The Demand Electricitv 

The previous sections presented a revie\r of the 

literature on forecasting b 7 statistical, econometric, and 

< 

stochastic time series methods. In this section we shall 
present a few methods, other than tlie methods revie\fed in 
the pre'vrious sections , \;hich have been utili?:sd for fore- 
casting of electricity,^ deiaand. 

Hooke ( 186 ) proposed a method for forecasting 
the demand for electricity hy class of service for a time 
interval of three to five years. A secular trend was 
dra^m from past data to forecast conditions to. fiye 

years hence. An estimate was made of expected deviations 
from trend in the near term future. By summing up the 
estimate of sales in each category and then adding losses, 
a forecast was made. Expected peak demands were determined 
hy a process of correlating peal^s vjith consumptions and 
system output, 

Goddard (18?) determined future electricity 
energy requirements by estimating the sales in each 
category of demand namely residential, conmercial and 
industrial. For each class the demand was estimated hy 
relating the cx>nsus 5 )tion to some economic and demographic 
variables . 
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Schiller (188) studied the demand pattern of 
different class of consumers by regression analysis 
tecliniques. In Electricite' de' France demand forecasting 
has been the subject of severcCL studies (189,190,191 , 19^) 
These studies have emphasised that growth in demand for - 
electricity is closely tied up vri.th technical progress, 
and economic expansion. The probability distribution of 
future consumption has been expressed as a function of the 
rate of general economic development as forecasted in the 
process of national planning. Ifodels representing succe- 
ssive demands increments by a Markov process have given 
satisfactory estir.iates of global demand. 

Arnoff and Chambers (193) have advocated the 
determination of peaJi: povrer demand from estimated KWH 
sales, since the error in estimating annual sales is 
normally less than error in estimating peak demand. 

The effect of marginal cost pricing for large 
customers, on peak demand has been studied by EDF (19^, 
195) • 'Studies indicated that there vra,s succe- 

ssive improvement in annual load factor due to adoption 
of marginal cost pricing principle. 

The study by NCAER (I 96 ) on determination of 
the future requirements of electricity for India has 
applied the 'method of end uses* techniques. The future 



demand has hr on obtained hy projecting: the outputs in 
various electricity consuriin;;:, sectors and then deter- 
mining the total requiTCrente of electricity from the 
knov/ledgc of specific electricity requireiients per unit 
of output. Similar techniques have been applied by the 
energy survey conni'ctce of the Government of Inc’ia 
to forecast future electricitj' requirements, 

A mathematical fomula requiring the use of 
population forecast only has been developed by Scheer 
( 197 ) of U,S, The formula has been developed from the 
thesis that for eveiy hundred fold increase in per capita 
generation the rate of gro\rtli in generation will be re- 
duced by half. This formula has been used for forecasting 
generation requirements in utilities , 

Mukhopadhyah and Tripathi (I 98 ) have advocated 
the use of flow graph analysis techniques to load fore- 
mens ting. Taking economic prosperity’- as a starting point 
a flow graph may be drawn relating growth of demand for 
po\;er due to the urge for economic prosperity. 

Methods for handling -weather sensitive demand 
has been discussed by several authors including 
Hieneman ( 57 ) f Thompson (199) 1 Davey et. al. (200) 
Ashbury- (201) , Clair and Siniirefclier (202) , Cbrpening 
et. al. (203) . The peak dec^ds and energy are separa- 
ted to weather sensitive and non-weather sensitive 
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component. Weather models arc determined by correlating 
demand with appropriate ^^eather variables. Wcatlier models 
are determined from data sampled more freouently. 

The literature survey presented in tins diapter 
is not claimed to be exhaustive, forecasting techniques 
are in an area of gradual develojonent and more and more 
research is being done on this area. However the more 
advanced techniques of forecasting to a large extent make 
use of a combination of simple approaches with modem 
developments of prediction theory. 

A surve3'- of the presented literature reveals that 
although considerable amount of research has been done in 
the area of forecasting the demand for electricity on the 
international scene, very little work seems to have been 
done in this area in the Indian context. This has motivated 
us to investigate and formalise strategies for long range 
probabilistic forecasting of annual peak demand and energy 
for India. In Chapter III, the methodologies selected for 
projecting demand for electricity in the present study 
haafebeen presented. 
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METHODOLOGIES OF ELECTRICITY DlJiA’ID FORECASTING 
TIHIR EVALUATIOII AND SELECTION 

3.1 i-ietl:odologies of forecasting 
EU^CTRICITY DEIAI2D 

The main types of teciiniqiics used for pro;jectinf; 

electricity demand can be di ’voided into the f ollo\\dng groups ; 

1 . Direct extro.polation of rates of change in requirements 
on the basis of observed characteristics of a past 
trend. 

2. Extrapolations fron observed relationsJiips based on some 
deduced functional relationships and causality between 
variables . 

3* Application of statistical and stochastic time ~ series 
analysis and theories of prediction. 

Methods of comprehensive plans - Extrapolation of 
specific energy requirements per unit of product to 
planned production. 

5, Input ^ output techniques . 

6. Comparative international analysis - Evaluation of 
future development from comparative study of its 
hcmologues in different countries. 
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In the following paragraphs we shall disciiss in 
brief the main features of the aforementioned methods and 
present an evaluation of these methods . 

3»1«1 ^rect Extrapolation On The Basis 
Of Observed Qiaracteristics Of A 
Past Trend 

TMs procedure is widely used for short term 
projections (up to one or two years aliead) of electricit 3 >^ 
demand. The basic assumption in all these methods of fore- 
casting by trend extrapolation is ^influences tliose were at 
work in the past vTill continue to exert a similarly evolving 
total effect in future. Two distinct variations of the 
method may be envisaged. They are: 

1 . Extrapolations maj’’ be carried out fran a rate of 
change deduced frcri past statistical data. 

2. Extrapolation by application of a mathematical 
formula vdiich fits development in tlie paist and is 
assumed to characterise the likely changes in 
future . 

To the same general class of mathenatical extra- 
polation techniques belongs, an attempt to fit an 'S» 
curve derived frcm expressions similar to a logistic trend, 
which accelerates to a point of inflexion, beyond which 
growth declines slowly and then becomes asymptotic* In 
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general straightline, parabolic, ex.)onential and modified 
exponential curves are applied. Out of these the one wliich 
appropriately reflects the gro\^}i is selected. 

Extrapolations from rates of change in past time 
series can be obtained from matrix tables which sum up all 
rates over a specified period. Guch methods normally 
assume the retention of a fair degree of stability of demand, 
Extrapolation can also be carried out using graphical 
methods. Graphical methods include assembling a series of 
long term trends as established from consecutive years in 
the past and plotting them all together from the same point 
as origin. The resulting cluster of graphs are stimmarised 
by plotting the median curve and upper and lower quartiles . 
This method is normally used as a starting point for other 
extrapolation techniques . 

3.1 *2 Extrapolations From Functional 
Rel ati onships 

This mode of projection inspite of its many 
variants offer a distinct and widely used extrapolation 
method. Electricity requirements have significantly high 
correlations -vidth certain econoroic variables. The deduced 
functional relationships between demand for electilcity 
and these variables provide the basis for extrapolation 
of demand. 
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St'uc33'-*of tlie many types of associative formulae 
shows that econonic indicators of the following types are 
most frequently incorporated in the relationshjips . 

1 • Demographic data of - 

(a) Total population. 

(b) Total labour force. 

(c) Employment in industry. 

2. Econonic data sudi as 

(a) Gross national product or national income, 

(b) Industrial production. 

(c) Fixed capital formation, 

(d) Per-capita Income. 

(e) Steel consumption. 

(f) Agricultural output, 

(g) Transportation indices. 

(h) Indices of coinmercial activity. 

The multiple regression technique which expresses 
linear and non-linear functional relationships is generally" 
employed to lirOi the demand for 'electricity with the above 
variables. The coefficients of such relatior^hips repre- 
sent the effect of unit change in one variable on the 
variable to be forecasted. 

Even \diile restricting the approach to linear 


regression equations based on method of least squares, 



many variations are possible. The values of variables may 
be specific (or absolute) or they may be represented by 
indices. Further, they can be expressed in linear coor- 
dinates or in logarthirdc form. It is evident tliat the 
periods for which functional, relationships are estimated 
should not be too long, if structural changes in parameters 
are lilrel^’- to occur, 

3,1,3 Forecasting By Stodiastic Time-Series 
Analysis .Ind Theories Of Prediction 

Any set of observations arranged chronologically 
constitute a time-series. The basic idea of stochastic 
theory of time-series analysis is to regard tlie time series 
as an observation made on an ensemble of random variables. 
Assioming that the peak power and energy demand are stochas- 
tic and stationary in nature the principles of time-series 
analysis can be applied to forecast electricity demand, 
Non-stationarity of data can be handled by converting 
them to a stationary series by means of appropriate 
transformations. The auto-correlation and partial auto- 
correlation function are used as tools ,to identify the 
nature of persistence. Cyclicity in the data, if 
any, are revealed and identified by applications of 
spectral analysis techniques. Having identified the 
trend, auto-regressive, moving average, and cyclic compo- 
nents In the time-series of past data and the order of 



the stochastic process, the parameters of the model are 
estimated by refilled statistical iterative tecliniques . 
After validation of the model to be fitted by means of 
diagnostic checking, the stochastic theories of predic- 
tion can be applied to forecast future values of demand 
and obtain tlieir range for a given level of significance. 

Comprehensive Economic Programmes And 
Aggregative Forecasts 

This method of projection seeks to cover the 
entire field of electricity supply and demand, where as 
the three methods discussed above are mainly adapted to 
determination of macro-level forecasts. Energy require- 
ments for a given economic system should be determined 
comprehensively so that the requirements of different 
consuming sectors are mutually consistent and in line with 
expected trends in conversion efficiency and substitution. 
Such estimates allow for effect of teclinical progress and 
social developments, which initiate structural changes in 
the growth of industries , 

Centrally directed plans and projections which 
adopt a ccanp3rchensive approach in respect of electricity 
requirements include the one-year, five-year and long 
term or perspective plans. Beqiiirements for electricity 
in the principal industries and consuming sectors are 



50 


worked out by opplyins to planned levels of output, the 
specific electricity input per unit of output. Separate 
sectoral estimates are aggregated and they are harmonised 
\ath loccl and regional estimates. An organised study of 
tlie nrobable interplaj' of factors making for substitution 
replacement is essential for the comprehensive method of 
forecasting. 

The metl'xod of approach in preparation of compre- 
hensive forecasts relics on the follovang three distinct 
type of inputs . 

1. Assembly of local and regional extrapolations by 
supplji" undertakings based on detailed special knowledge 
of requirements. 

2. Careful studies by individual consuming sectors 
(market studies in residential sector, production 
forecasts in industrial sector, etc.). 

3. Macro studies linking electricity re quire i^ients mill 
indices of economic gro^rth. 

The comprehensive forecasts for electricity 
demand is based on the reconcilation of the information 
obtained frean the above mentioned inputs. 
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3 • 1 • 5 Forecastlnr: By Input-Oatput Tedinicues 

Tlie input-output model is represented by a system 
of simultaneous equations \;'iich lirJ: output of all sectors 
and industries to the inputs used in tliose sectors. It is 
utilised to anal 3 'se and predict all the productive inter- 
industiy transactions (by industry'- categories) and others 
that go into determination of gross national product. The 
input-output methodology provides a means by v/hich economic 
forecasters can convert their estimates of aggregate final 
demand to estimates of total required output in every sector 
including electricity and amount of resource inputs that 
are consistent i/ith GITP projections . 

3.1.6 Pro .lections By Comparative Analysis Of 
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income and electricity consumption per capita for a wide 

* 

range of countries. This information can be utilised for 
the forecasting of electricity demand. 


3,2 AN EVALUATION OF THE lETHCDOLOGIES OF 
ELECTRICITY DEI-AIID FOBECASTI?]G 

In the previous section we have presented six 
different tecliniques of forecasting electricity demand. 
These techniques do not necessarily exclude one another 
completely. They all reduce to some logically consistent 
mode of extrapolation offered by observed regularities of 
past behaviour. The inherent assumptions are the conti- 
nuation of the mode of evolution of factors generating the 
demand for electricity. In each tedinique provisions have 
to be made for uncertainity and unanalysed interplay of 
various technical, social and economic factors by use of 
informed judgement. It has been observed that despite some 
common features tliese methods of projection necessarily 
differ vridely in their application and use. Each method 
has its distinct advantages, disadvantages and problems 
of application. 

The two most important assumptions in the trend 
extrapolation approach are that the growth process conti- 
nues to remain same tliroughout the historical period and 
that ihe factors responsible for grovrth continue to evolve 
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in the same manner. These assimiptions distinctl 3 ’- appear 
to he implausible . There can he no assurance that forces 
which influenced past growth iriLll continue unaha.ted and no 
ne\; factors v/ill he introduced. Ragnar Frisch (20if) an 
eminent econometrician calls tj:e trend extrapolation method, 
the most primitive technique of forecasting. It needs to 
he pointed out that the trend extrapolation method does 
not generallj’’ take into consider'''tion major structural and 
investment changes in the economy. Cases may exist when 
during past years there might have heen unusually liigii 
rates of growth due to a low datum of electricity consump- 
tion. Extrapolation from such high rates, which usually 
are non-linear, i.dll resiflt in unusually high estimates of 
demand, \/hich raa 3 ^ not at all he realised. Hence results of 
trend method need to he modified and revised hy exercising 
a great deal of judgement in anticipating future conditions. 
The application is limited to short range forecasting. 
Despite these limitations trend extrapolation technique 
certainly is simple and easy-to-apply . It has the advantage 
and utility in providing a convenient first approximation. 

The difficulties •i.rhich underlie the econometric 
methods are numerous. The structuring and f emulation of 
a comjxLete and detailed simultaneous equations ecceometric 
model representing the realistic linear and nonlinear and 
coB^slex technical, economic and behavioral relationships 



is a difficult task. There always exists a possihility 
that the assumed form of functional relations liips do not 
represent the forecasting situation realistically. In 
addition there are problems of identification, specifica- 
tion, choice of appropriate independent causal variables 
and errors in data and their measurer.' ents . The inherent 
and fundamental assunpeions in the metliodology, of non 
existence of multicollinearity in independent variables and 
autocorrelation in errors, wliich are difficult to realise, 
create problems of estimation. Even though long series of 
data are advantageous and desirable for purposes of estima- 
tion, the distinct possibilities of structural changes over 
tine may not be ruled out. Parameters and coefficients may 
turn out to be unstable and non-stationary, requiring fre- 
quent revisions. The forecasts obtained by the application 
of this method can only be relied upon if the structural 
coefficients are stationary, proper regressors are chooseii 
and reliable estimates of future values of the stochastic 
regressors exist. Since the econoEietric models are essen- 
tially static in nature, forecasts need fairly frequent 
amendments, Hotmtlis tending above problems, econometric 
methods are emerging as excellent tools of fore ceis ting 
the demand for electricity, since they explicitly take 
into account the significant causal relations that exist 
between electricity demand and eccnomic variables. 
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nUie facility afforded by stochastic time-series 
analysis approach for choice of a forecast function 
appropriate to tlic particular probleri under consideration 
is the nost distinctive features of t3iis ncthod. Never- 
theless it is just this characteristic wliich has led in 
sore instances to criticism of the approach and difficul- 
ties in its practical application. The method provides a 
cla.ss of models and strategies by \daich particular models 
can be chooson from the various classes wliich reflect the 
properties of individual tine scries . The eventual 
forecast function is dictated to a large extent by the 
data used. The models capture tiie totality of information 
coming out of the data. Further, in contrast to tsrpical 
exponential smoothing techniques , this approach allows a 
good deal of freedom and judgement in the choice of an 
appropriate model. The method has been criticised on 
grounds of requirements of large amount of time for analysis 
and need of better expertise . The freedom of choice offered 
by the approach may lead to the selection of a poor model 
producing forecasts of low quality. The method does not 
guarantee numerically identical conclusions from the sane 
set of data by different analysts, A moderately long 
series of data are required for its satisfactory applica- 
tion, S'opliisticated and complex computational procedures 
are needed for model identification, estimation and 
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conter.iplated since there exists several potential area^ 
\dierc these models are likely to find applications. 

The method of comprehensive forecasts is an 
aggregating, com^'^onsating and cross-checking approach 
vdiich replaces the reliance on the relative assurance of 
an economic plrji. It uses overall economic forecasting 
combined with aggregation of a large number of sectoral 
and regional estimates of electricity needs. However, for 
all countries the long term development and structural 
economic changes pose problems which cannot be forecasted 
once and for all. Perspective plans are controlled by 
macro-level projections prepared witli variations allowed 
for alternate developments . The formulation of energy 
needs can only be arrived at in stages since they are . 
mainly derived from output plans. I&iderlying iiiis metho- 
dology are the problems of obtaining data and predicting 
the specific needs per uaiit of output in future. 

Projection by ccanparative analysis of hopologue 
trends is by far the least well understood and documented 
of projection techniques. The major problem in using 
this methodology arises due to the uncertainities invol- 
ved in forecastiiig the nature of long term economic 
develo^ent. The use of the comparative metiiod to study 
demand relationships in different countries and to make 
projections based on these relationships, aids in 
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obtaining improve-'" forecasts hy reducing the apparent 
disjDersion betv/een tendencies in aifferent countries and 
by all owing projections to be made more realistical.ly . 

This however is only tnie if the data clioosen are linked in 
a related development process. Such a development process 
can be discerned in the world phenomenon of characteristic 
growth rates for electric energy, which might be explained ' 
in terms of clianging rates of capital investment and evol- 
ving levels of labour productivity and personal consumption. 
The value of t’lis methodology is implied in the emphasis 
given to international comparability in various surveys of 
electricity demand. Characteristics varying with structure 
of national resources and tlie resulting pattern of indus- 
trial structure in different countries are less likely to 
yield useful resiilts except -idiere their structural features 
can be matched. The difficulties that underlie the applica- 
tion of this methodology is in procuring, organising and 
co-oitlinating the international data. 

A comparative/relative evaluation of the metho- 
dologies presented above lead us to the conclusion that 
none of the methodologies are entirely satisfactory and 
fool proof. The usefulness of any particular method 
rests on tiie nature of problem in question, availability 
of data and considerations of time and resources. 
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3.3 SELSCTION OF ^ETHODOLOGIES FOE 
TilE PEESEIJT STUDY 

The selection of one or more methodologies for 
forecasting is influenced by various factors. Amongst 
several considerations the availability of suitable data 
for the past, and their tine span are major factors, in- 
fluencing the choice of a methodology^. An ideal projection 
technique shoiild take into accouiit the influences of many 
interrelated factors of rates of growth, of changes in 
structure of activities and technical progress; since 
projections allowing for these factors can be checked for 
consistency. Further, due to tlie uncertainities involved 
in forecasting, forecasts obtained by the application of 
any methodology should be specified for a given probability 
level by defining upper and lovrer limits of forecast. An 
additional factor that should be given due attention while 
selecting a methodology for a given forecasting situation 
is the adaptability of the methodology to the nature of 
forecasts , 

A relative evaluation of the methodologies indi- 
cate that lujne of the methods of projection are entirely 
satisfactory for a given application. However, after 
careful consideration and evaluation of the relative 
merits and demerits of tl:a six methodologies presented 
we have selected the following two methodologies for 
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the pvirpose of the present study. They are 

(1) Econometric me thodolog:/ 

(2) Stochastic time -series analysis and forecasting 
methodology. 

In tlie follovn.ng paragraphs vre provide reasons for the 
selection of these methodologies. 

3.3*1 Reasons For Selecting The Econometric 

Iletliod For Electricity Demand Forecasting 

It is evident that the demand for electricity has 
a significantly’- lii{^ correlation \ri.th t}ie levels of econo- 
mic aiid industrial activity as well as with certain 
demographic variables . I-fecroeconomic theory pro-vides us 
the overall frenev/ork for forecasting an econoroic variable 
(lil:e elec trie i-ty consumption), by providing a sound 
hypothesis which explains economic behaviour, and based on 
■v-diich causal variables can be chosen.. It is the econome- 
tric method which explicitly considers the functional 
relationship bet-ween causal economic variables and the 
demand for elec tricity, and also provides a procedure for 
quantitative estimation of these relationships. Further 
it needs to be pointed out that the implications of any 
policy, which is normally required for planning and devi- 
sing policies, can only be determiited by formulation of 
a model which recognises the effects and interactions of 
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econondc variables on demand for electricity. The applica- 
tion of the method, of extrapolating the demand for 
electricity, from quantitative econometric relationships 
are rare in the Indian context and projections made expli?*,. 
citly as prohaLilistic state) ^ents referring to a particular 
poiiik in time are iiostly unavailable. It is these afore- 
mentioned reasons, along idth the availability of reasonably 
good quality data, viiich has motivated us to use this 
technique for electricity demand forecasting. 

3.3.2 Reasons For Selecting The Stoche.stic 
Time-Series Analysis And Forecasting 
Technique 

The following considerations have provided us a 
motivation to use this methodology for electricity demand 
forecasting. 

Tl:ie theory of stochastic time -series analysis 
proposes a general class of models oat of idiich particular 
models can be choosen according to the properties of the 
individual time-series data of peak and energy demands. 

The eventual, forecast function is to a large extent dic- 
tated by the data, a property referred to as "letting the 
data Speak for itself". The method captures the totality 
of information that can be derived out of a time series. 

As pointed out earlier, the present method offers a good 


62 


deal of freedom and choice to exercise judgement in selec- 
tion of appropriate models , V/ith the availability of more 
data models can be modified and revised through a more 
thorovigh analysis. 

Stochastic methods involve the development of 
probabilistic models idiich have as their output the evolving 
electricity dema.nd, Tliis nethocblos'- offers us models 
possessing naximui:i simplicity and minimum number of para- 
meters concordant \-rith representational adequa.cy. The 
obtaining of such models may be important as thejr nay provide 
us information about nature of the system generating the 
time-series. Further once an appropriate model is built 
optimal forecasts can be prepared rapidly from past and 
current data. 

In the present dissertation the above mentioned 
two methodologies have been applied to the problem of 
electricity demand forecasting for India. The details 
regarding the econcaaetric and stochastic time-series ana- 
lysis and forecasting metliodologies have been presented in 
Chapter I? and Chapter V respectively. 


CHAPTER IV 


FOrJBCASTIUG OF ELECTRICITY DEI JO® 

BY ECOHOIIETRIC ^ETITODOLOGY 

This chapter is devoted to the development of 
models for forecasting the future recuirements of electri- 
city by the applicotion of econoeietric methodology. We 
propose macrolevjl econometric models to analyse and 
estimate the demand function for electricity. The models 
have been applied to forecast future requirements of 
electricity in India for the post sample period (from 1976 
to 2000 A.D.), In section ¥.1 we present the macro econo- 
metric model for forecasting future peak demands and energy 
requirements. Section 4,2 deals with the formulation and 
estimation of models for analysis of industrial demand for 
electricity. Results obtained by the application of the 
methodologies and their analysis are presented in Chapter VI, 

4.1 EXOGENOUS VARIAEIES AFFECTING THE DEMA!® 

FOR ELECTRICITY 

Electricity being a factor of production has 
a derived demand. By principles of economic theory x^re 
postulate that th.. demand for electricity in an economy 
depends positively on the measures of economic activity, 
which use electricity. E<-onomic and demographic variables 


which are expected to have a significantly high correlation 
\jith demand for electricity are gross national product 
(national income), real per capita income, price of elec- 
tricity and price of their substitutes, industrial output, 
agricultural production, transportation activities etc, 
Electricity is utilised in all sectors of the economy and 
coal, fuel oil, gas and other forms of energy resources are 
■used as its substitutes , 

The choice of macro causal variables (which are 
to be used as independent regressors) effecting electricity 
demand can be confined to gross national product. Industrial 
output, poptilation, per capita i ncome, price of electricity 
(relative price) because variables like agricultural pro- 
duction, transportation activities are indirectly included 
in gross national product values. 

Econometric Methodolog y - 
Its Tools And Problems 

Before presenting the actual models, we discuss 
in brief the econometric methodology - its tools and 
problems , 

Three distinguished economists, Sanu^lso^, 

Koopmans and Stone (20^) have ably described econometrics 
as a quanti-tative analysis of actual economic behaviour 
based on concurrent development of theory and observation 
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related "by appropriate methods of inference. Econometric 
studies involve systematic collection of information about 
actual phenomenon and construction of empirical correlates 
to theoretical concepts . Tliis activitj* dra\;s mainly on 
mathematical statistics and economic theory. Hie objective 
of an econometric model is to obtain quantitative state- 
ments that explains the behaviour of variables already 
observed anc’ forecast their future values. Econometric 
models are essentially stochastic wiiose distinguisliing 
feature is the presence of a disturbance or error term. 

The models are generally in form of single or a 
system of linear and non-linear equalities and inequalities 
with numerically specified coefficients . An ideal econome- 
tric model should have properties of relevance, simplicity, 
theoretical plausibility, explanatory ability, consistency 
and accuracy of coefficients and forecasting ability. 

The econometric model building and its validation 
is generally accomplished by 

1) stating the problem clearly, 

2) choosing an appropriate model or maintained 
hypothesis to be used for the problem, 

3) observing relevant data, and 

if) using st'tistical infeience techniques based on 
the maintained hypothesis and drawing inferences 
about the problem from data. 
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The difficulties underlying the above mode of 
attack arc to obtain accurate knov^lcdge about the appro- 
priate model. The main source of maintained hypothesis is 
economic theory v:hich at best specifies the conceptual 
variables to appear in an equation, the algebraic sign and 
value of certain partial derivatives and the functional 
form, but tljls is not possible i;i all sitviaticns. In such 
an eventuality it is recommended to choose a functional 
form of equation i;hich is reasonable Cii gro'ands that arc 
in part theoretical and empirical. One d^oicc is to brj'' 
out several different theoretically plausible functional 
forms in a sort of experimental fashion and select the one 
amongst them vMch fits the data veil. However this method 
is beset, v/ith certain problems. These problems stem frcci 
statistical inference and description of temporary or 
accidental features of data rather than enduring systematic 
features , 

^.1.2 D evelo p ment Of Econometric 
F orecasting ^^:^el For 
Electricity Demand 

Wc specify the following implicit functional 
relationships between electricity demand and the chosen 
independent varia'^les, 

^ = 




( 4 . 1 ) 
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PD^. = (GNP^, I^, P^, Pe^, Psj., t) (4.2) 


W}icrc, 


E 


t 






Annual energy consumption (demand) In year t 
in million KWH. 

Peal: demand in year t in IIW. 

Index number of gross national product in 
year t. 

Index number of industrial output in year t. 

Population in year t (IlUions). 

Price index of electricitj' in year t. 
Relative price index of substitutes of 
electricity in year t. 


In the above relationships it is assumed that all money 
oriented indices are deflated by a suitable price index 
to neutralise effects of inflation in jthe variables. 
This would resiilt in indices at constant prices , Demand 
in this study is defined to be identical to consumption 
or satisfied (effective) demand. 


The demand functions specified by equation 
(4.1 ) and (4.2 ) are implicit. Tlie various econometric 
models studied are given below. A relationship linear 
in parameters and vaidablcs has been asstimed. Four 
variations of the linear model are proposed. The first 
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model is of the sijnple linear forr-, Model 2 incorporates 
first differe/ices of variables. The third model is of the 
double log-linear fonri. Model ^ incorporates first diffe- 
rences of log transformed variables. Models proposed for 
peal: poi^er and energy demand are identical. The models 
are as follOA.^ : 

Model 1 

^t = f^o ^2 ^t ^3 ^t ^®t 

+ ^5 t + 


®t = 

Po 

•f 

'^1 

GNP^ + 

1^2 h* 1 

A 

-3 

* ''’I ^■ 





PSt + 

'*6 ♦ <='t 



Model 2 








AE^ = 




A-GNP^ 


+ 

P 3 




Pi 

ZiiPe^ + 

+ 

» If 

1-6 * + t 

APD^ = 

/3o 

+ 

^1 

^NP^ 

*** 

+ 




t 

K 

APe^ + APs^ 

t 

Pe t + g ^ 


Ifodel 3 

Log = a-Q + log GNP^ + a2 log log 

* % log ^6 

= + a!| log GNP^ + a^ log + a^ log 

+ log Pe^ + a^ log t + 


(^.3) 


CifA) 


Log PD^ 
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Model 3 + 


Alog = 


^Icg PD^ = 


+ b^ Alog GNP^ + b^ Alog + b^ Alog P^ 

+ bj^ AJog Pe^ + b^ Alog Ps^ + b^ t + 

(^. 6 ) 

b^!^ + b^ Alog GNP^ + b^ Alog + b^ Alog P^ 

+ \ Alog Pe^ + b^ Alog Ps^ + b^ t + uj. 


where , 

represents first differences of variables, 

/^ , a, b, etc. are the coefficients of models to be 
estimated, and 

^ , 11, etc, are the error terns. 

The models proposed above include all the independent 
variables. Various other models can be formed lath alter- 
native combinations of variables. The tine variable is 
included as a surrogate for all tine-trended variables like 
rate of electrification, and technological progress, etc. 

Ho published data on electricity price and prices for 
substitutes of electricity are available. However, data 
on relative price of electricity vdth respect to fuel 
price index are available, Tlierefore, the variables 
Pe^, and Ps^ are replaced by the relative price index, 
since it is construed tha,t this ratio reflects price of 
electidcity, and other fuels. Further as most of the 
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variaMcs arc significantly time -trended due attention has to 
be given to the problem of mnlticollinearity in indepen- 
dent variables • The problems of statistical estimation in 
presence of imilticollincarity and autoc or related error terms 
have been dis classed b-^' Johnston (60) , Theil (65) , 

Tintncr (69) and others. 

^.1«3 Data Used For The Problem And 
The Estimation Procedur e 

The time series data for the period 1928 to 1975 
on electricity consiamption, peali demand, GUP, popialation, 
industrial output, relative price of electricity have been 
collected from various sources. The sources of the data 
arc listed in Appendix. Except population all other data 
for the independent variables are expressed in terms of 
indices with base of 1961 - 1962 equal to hundred. The 
data on electricity consumption arc expressed in million 
KVCI units and peah demand in Ilogavatts, 

The ordinary least squares technique has been 
used to estimate the coefficients of the equations . A 
computer program has been written in Fortran TV based on 
the linear least squares algoritlm. 


age number 71 has been omitted by mistake in the sequence of page numbers. 
Please read page number 72 in continuation of the matter in page number 70. 
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4 , 1 .4 Forecastlnr The Future Values Of 

Electricity Energy And Peal: Demand 
Eased On Estimated Functional Relations 

To obtain reliable forecasts of future require- 
nents of electricity energy and peak demand on the basis 
of an estimated relationship \ic need to have reliable 
estimates of future values of independent variables . 
Reliable estimation of future vodluas of C-NP, population, 
industrial output is a ciifficult task, as these variables 
depend on a variety of economic, technical and social 
factors vhich are highly stochastic in nature. Ilethods 
for prediction of economic variables have been discussed 
by ICLein (61) , Theil (63,64) , Buttler et. al. (76) 

and various other authors. Tne application of tliese tech- 
niques for prediction of future economic and demographic 
variables are beyond the scope of this work, 

V7e have resolved the problerc of estimating the 
future values of relevant exogenous variables in a rela- 
tively simpler manner. It is assui-ied that the past rates 
of grov/th in these variables \n.ll be maintained in future. 
The future values of the independent values are extra- 
polated from these past grov?th rates. Another choice is 
to use the planned and anticipated gro^rth rates for these 
variables for extrapolating their future values. Substi- 
tution of these future values in the estimated relation- 
ship provides us the point forecasts , The upper and 
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lover ranges for these forecasts liave Leen obtained for 
a given probability level. 

The vslues of forecasts c" tsined for India for 
the period 1976 - 2000 ai-e presented in Chapter VI, 

k.2 Ai^ALYSIS OF IITDI}SIBI J. DI]:ih.~ID 
FOR ELECTBICITI 

In this section ve present econometric models 
for analysing the demand for electricity in tlie industrial 
sector. The approach adopted involves construction of 
models incorporating hypothesis based on accepted economic 
principles, Tlie quantitative estimates of various struc- 
tural coefficients are obtained using the historical data, 
Tiie objective of this anal 3 rsis is to obtain a clear under- 
standing of the basic factors influencing past developments 
of industrial electricity demand, Tlie proposed models may 
be used for forecasting electricity demand for the indus- 
trial sector. 

The analysis has been carried out for the 

following groups of industries viz. Iron and Steel, 

Textiles, Pape]p, Vehicles, Engineering, Mining and 

Quarrying, Nonferrous Metals, Chemicals and Food in 

Indian context. The various models proposed, are based 

/ 

cai postulates, idiich describe relationships between rele- 
vant exogenous variables and demand relying on macro- 
econoiBlc principles. 
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4,2.1 Foimialotion Of Hypothesis ^oout Doteniiinants 
Of !berand For Electricity In Tlie Industrial 
Sectors 

The approach adopted in bi’’ ' study is to treat 
electricity as an input enterinf into the nroduction func- 
tion, Using this philosophy ue develop specific hypothesis 
about the form of relationsiiips involved and particifLar 
variables whidi appear to be relevant. \-le will describe 
three basic models each derived from a differing primary 
hypothesis. 

In the industrial sector the maiii features of 
electricity consumption in the past has been that consump- 
tion of electricity grew faster thsin the output in every 
industry. The following three explanations can be put 
forth for this phenomenon, 

1 . Relative price movements may induce substitution of 
electricity for other fuel inpiits and possibly labour, 

2. Technological changes may be such that they lead to 
innovation based on use of electricity. Further 
technological changes may also induce substitution of 
electricity against other fuels and labour. This 
substitution is attributed to the technical advance- 
ments in the usage of electricity rather than a 
favourable price movement for electricity. 
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3. The normal explanation given hy economic theory’- for 
non-proportional variation in input and output is the 
existence of varying returns to scale. Electric power 
does not fit into this f rar.iewor]-: since its relations Mp 
with output is indirect thrcu.gh tl-e meditu.i of plant 
and machinery. The d /iiand for electricity is a derived 
demand and hence it is effected, by the demand for out- 
put of all electricity us.inf commodities. In the long 
run therefore the variation of electricity demand with 
variation in output may depend on the follo\ 7 ing two 
factors , 

i) Factors ^diich determine how cayital stock varies 
as output varies , 

ii) Factors which determine the electricity using 

characteristics of increments in capital stock. 
Thus when we observe that electricity demand grows 
faster than industrial output, this would suggest the 
follo\ri.ng: 

(a) Greater proportional increase in capital stock 
than in output might have led to faster growth . 
of electricity demand. 

(b) Tlie electric intensity of the additional capital 
stock reeded for increased output might be above 
the average for the existing stock. 
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Thus ■when relative fiJll price and tedinological change do 
not appear to have influenced the grovrbh of electricity 
demand, faster gro\/tli of electricity demaiid as compared to 
output may be due to more than proportionate increase in 
capital stock as compared to output. The absence of ade- 
quate data iniiibits the explicit inclusion of electricity 
intensity and capital stock relationships in the develop- 
ment of forecasting models . 

The three models proposed in the folio-wing sec- 
tion are based on liypotheses, which use the first t-wo 
explanation. u given in this section for the higher rate 
of consumption of electricity as compared to the industrial 
output • 


^•2,2 Models For Industrial Electrici 'tv’^ 
Demand 

Model 1 


The development of this model is based on the 
straightforward application of the tneo-’y of demands for 
inputs. We assume that firms have a production function 
of the Cobb - Douglas (206) type > defined on inputs such 
as labour, capital, rai^ mate rial *5 and various forms of 
energy. Let the output of an industry be Q. Then 


Q = 


<5^0 






• ♦ ♦ 


(^. 7 ) 
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where 

= relevant inputs, j = 1, 2, k 

^ = parameters for j = 1, 2, ... k 

It is assumed that firms v.'ish to mnirnisc the total cost 
of production for any output. T.et the total cost of 
production he C, Then, 


C = p^ x^ •{' P 2 ^2 + . .. + Pj, X,. (^.8) 

where p.j, P 2 , Pj^ are the prices of inputs x^ , X 2 , 
... Xjj^ respectively. 

Minimising (^.8 ) subject to ( ^*7 ) gives tlie first 

order conditions for a constraint cost minimum 


P1 - X0C1 Xo ^ 


^'2 *^3 


0(1 ■^2 


- 1 


P2 “ o<^ X2 


o<'3 


^3 


• ♦ Xi- ~ 0 

x \ 

^k 


^1^2 

Pk “ k ^o ^2 






s. - ■' 




f cy ^ 
•1 '2 


c/. 


Q -cXq xi 


= 0 


= 0 




\diere ^»is the Lagrange multiplier associated -vd-th 
(eq. 4 , 7 ). We now have (K + 1) equations in (K +1) 
unknowns which can be solved for obtaining the values 
of ^ Xgj Xjj,. 


• • ♦ 


Let us assume, that electricity is the K input. Then 
the demand for eloctilcity, is given "by the expression 

^'1 4 h [c+1 

P 2 ••• Pp Q (4- .10) 

idle re 

= parameters of p^ •••Pi, ,1 = 1, 2, ... k and 

/5 .. parameters represent combinations of c<!, . , j = 0, 1, ,. .1 

Equation ( 4-.10 ) shov;s that demand for electricity is an 
exponential function of K input prices and the output. 

The demand model given by Sq . ( 4- , 1 0 ) forms the bas is for 
the set of models given by equations ( 4-,11 ) to ( 4-.14- ) 
for obtaining electricity demand. 

In deciding on the appropriate variables to be 
included in the equations for demand, a reconciliation is 
necessary betireen the needs of the hypotheses and the 
limits imposed by statistical estimation theory. Equations 
vith the most desirable economic preperties and explanatory 
ability may well have undesirable statisticaJ. properties, 
particularly because of multicollinearity between indepen- 
dent variables, Wlien the risk of bias in parameters are 
felt to be particularly acute it is preferrable to take 
theoretically less satisfactory conbinations of variables . 

With the above considerations we propose the 
foUoidLng group of four equations based on model given 
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by equation (4,10 ), The linear form of equation ( 4,10) 
is obtained by ta]:inc logarithms of both sides. The 
price variables chosen are : 

(1) price of electricity relative to fuel price index, 

(2) price of electricity relative to vage rates. 

Each of these variables along \7ith output have been included 
in seperate equations to mininise the effect of imilticolli'- 
nearitj' betvrcen independent variables on bias and loss of 
significance of parr.ncters estimated. 

Log + /?-, Log + 6 ^ ... (4.11) 

Log Log Log ( ^)^ + (4.12) 

Log ^ Log + 7^*2 Log (g)^ + (4,13) 

Log Log + ft 2 ^ t + (4,14) 

where 

= Electricity consumption in KVJH 
= Index of production 
P = Price of electricity (Es.) 

F = Price index of all other fuels 

W = Average Wage rates (Es.) 

t = Time -(Years) 

= The disturbance term 
values are the parameters. 
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Ibdel 2 

I'bdel 1 lays emphasis on relative price and out^^iut 
as explanatoiy variables , Model 2, though complementary to 
Ibdel 1, lays emphasis on changes in fuel technology. 
Generally in studies relating to inclusion of the technolo- 
gical progress as a relevant independent variable, the 
time trend is tised as a surrogate for technological progress. 
For the present study this surrogate variable is undesirable 
both from statistical arid economical point of view. Since 
most variables used in the analysis are strongly time - 
trended inclusion of time as an independent variable will 
lead to multi coll inearity and consequent statistical problem 
of estimation. Further time as a variable does not explain 
anything of the process by which technological cliange taJ:es 
place. Instead it is felt that in many industries the most 
important change due to advancement in fuel technology 
would be reflected in the declining use of coal. Hence it 
is our contention that coal and cohc consumption for the 
industry can be included as a surrogate for technological 
improvement. 

On the other hand some of the substitution 
miglit have been due to relative price changes, Shat is 
why Model 1 and Model 2 are complementary. If Model 1 
suggests that relative fuel price had little effect 
on electricity demand, the coal variable will reflect 


ai 

influence of change in fuel teclinoloey. In case electricity 
der.'.and appears to have been influenced by relative j)rice 
changes, the coal \’ariablo \7ill represent total effect of 
price and teclinological changes. Lvonthougli inclusion of 
the fuel price relative and coal variable together may 
separate out the two effects, the problem ol statistical 
distortion due to intercorrelation bet’ween these two 
variables mahes tn^is approach unattractive. 

The follo^.dng groups of equations ( 4,1 5 ) to 
( 4,19 ) arc bs-sed on the above considerations, i.e,, coal 
consumption is an effective surrogate for tec}inological 
change, Purtlier to tost the effect of lags one period 
lagged coal consumption is included in equations (4,15, 4.16) 
The other major variable is the industrial output. An 
attempt has been made to study the effect of inclusion of 
a time trend in equation ( 4.16 ). Empl oyment/ output 
ratios are included in equation ( 4,18 ) to isolate the 
labour intensity aspect of change in technology. Invest- 
ment/output ratio has been incorporated in equation ( 4,19) 
to isolate the capital intensity aspect of technological 
progress . 

\ = *=^0 "^1 ^2 “^t ^t,l‘ ^t (4.15) 

2t = «^o *^2 ^t "^3 ^t-1 ^ 
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~ o 1 ^t '*’ 2 '’’ ‘^ 3 ^ t (^>-*10) 

^*^2 + niy^ ^"b (^i"»”19^ 


Vfiiere 


"t-1 


Q 

I 


OfO 

^t 


t 


Electricity cons^’nec’ in coal equivalent units 

Coal consunec" in period t 

Coal consiuned in period t-1 

I-Iuniter of employees 

Index of industrial production 

Gross fixed capital formation 

values are parameters of the equations 

values are disturbances 

Time. 


Model 3 


This model is constructed on the hj’pothesis that 
there exists a one to one relationship between changes in 
output and electricity consumption. Tlie deviation from 
these changes are assumed to be induced by dianges in 
relative prices, changes in labour intensity aiK?. changes 
in capital intensity. 

It is assumed tlaat demand for electricity is 
directly related to output by the following relationships 


• <k 


f • • • 


n 
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t = 1, 2 


(4.20) 


oC^ = r ( ... ) 


(4.21) 


]fJliere 

= Demand for electricity in I'.VJrl in time t 
= Proportionality factor betv^ei. outp-ut and 
electricity demand in period t 
= Output of the industry at time t 
^1f ^2t. ••• ^mt = Values of relevant independent 

variables such as price of different fuels , 
labour intensity, and capital intensity. 


From equation ( 4.20 ) ve get, 

(2), (4.22) 

To study the effects of relative price of different fuels, 
labour intensity and capital intensity on consu3|iption of 
electricity per unit of output, ve propose the follov/ing 
equations. Equation ( 4,23 ) studies the effect of labour 

r 

intensity by including the ^ variable . The influence of 
capital intensity on electricity is determined by equation 
( 4,24 ) by incoi^orating the ^ variable. The impact of 
relative price of different fuels represented by tdie 

p 

variable ^ is obtained from equation ( 4,25 )- Effect of 
price of electricity relative to wage rates on electricity 
d®nand is studied by equation ( 4.26 ) incorporating the 

p 

^ variable . 




Log ( I 
Log ( ^ )^. 
Log ( I 
Log ( Q 


^°S ( f )t + -_t 

% + Los ( 5 )t * 4t 

a© + Log ( f )t ■*■ 

a© + a^ Log ( y )^ + *1^^ 


(^. 23 ) 

(4.24) 

(4.25) 

(4.26) 


■v^ere, 


M 

I 

Q 

P 

F 

W 



= Kimil'er of employees in the industry 
= Gross fixed capital stock formation 
= Output of the industry 
= Price of electricity 
= Fuel price index 

= Wage rates of the workers in the industry 
etc . are the parameters of the equation 

are the disturbance temi'-i . 


4.2.3 The Use Of Time Lags 

The equations (4,11 ) to ( 4.14 ) and (4,1 5 ) 

to ( 4,19 ) presented above have been reformulated to take 
into account tire possibility that tlie equilibrium level 
of the dependent variable will not be attained until some- 
time after the changes in independent variables liave taken 
place. The following pattern of adjustment of dependent 
to independent variable is assumed. From equation (4.12) 


of riiodel 1, we can express the equilibriur; demand as 

= fCQ^jCp)^) (^r*2V) 

It should be noted that tills level of demand is only rea- 
died in the lone run equilibrium. Assumine a flow adjust- 
ment cype of equilibrium we have 


’I-i 


D* 1 

( ^ ) 

"^t-l 




(if .28) 


v/here X is the parameter of adjustment. 


actual demo.nd 



is given by 

0->) 

/ 


Therefore the 


(if.29) 


Tailing logarithm of both sides of equation ( 4.29 ) and 
substituting for from equation ( 4,27 ) we have 

Log = A log [ f { (|)^ )J + (1 - A) Log 

(>*.30) 

Assuming that equation ( 4,27 ) is exponential tlie long-run 
elasticities can be obtained by simply dividing through 
by an estimate of X from equation (4,30 ). It can 

be proved that the weights in the lag distribution are 
geometrically declining. For proof reader is referred 
to IIyock*s work (207 


h , 2 Data Used For T!ie Study And 
T3ie Estii^c~t~ion Procedure 
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For Iron and Steel, Textiles, Paper, Vehicles, 
Engineering Kinirig and Quarrying, Ilonferrcus Ifetals, 
Qiemicals, Food industries the relevant da.ta on the various 
variables Tor th.e period 196I - 71 has been obtained fron 
various sources. The sources of data are listed in 
Appendix A . The ordinary least squares technique has 
been used for estinaticn of the structurrl. coefficients. 

The results of this analysis, their interpreta- 
tions, and conclusions drawn on the basis of the results 
are presented in Chapter VI, 


CHAPEER V 


STOaiASTIC TI^E-SERIES ANALYSIS AI® 

FORECASTING OF ELECTRICITY DEMAl® 

In tills chapter we shall present the laethodolofiy 
for 'building, identifying, fitting and checking time-series 
models for electricity demand forecasting. The methodology 
discussed is appropriate for discrete ( sample d-dat a) sys- 
tems . The proposed time-series models are used to forecast 
the future values of peak power and energy demand for 
India for the period 1976 to 2000 A.D. The motivations for 
applying stochastic time-series analysis for electricity 
demand forecasting has already been presented in Chapter 
III. 

The process generating the time -series of electri- 
city demand is a dynamic S 3 rstem exhibiting changes in 
characteristics over time. The time -series of electricity 
demands (peak power and energy) are sets of observations 
generated sequentially in time and they constitute a 
statistical phenomenon that evolves in time according to 
probabilistic laws. For such a time series prediction of 
the future variations of demand can only be made based on 
probabilistic models . Hence it is contended that the 
theory of stochastic time series analysis and prediction 
can be applied to time -series of electricity demand. 
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For an extensive exposition of the theory and 
analysis of stochastic process references are made to 
Doob (145), Kannan (135) » Wiener (l4), VJold (136), 

Qnenouille (1^9), Parzen (14-6), Bartlett (14-3), Box and 
Jenkins (4-5), Jenlcin and V/atts (209), V/hittle (138), 
Robinson (211), Anderson (14-7), Grenanden et. al (212), 

Rao (213), Wald (210), Bailey (1^-4) ,• Bartholomew (150) and 
othe rs , 

In the folloi^ang paragraph a brief statement of 
the forecasting problems i/hen viewed from time -series 
analysis point of vie’.; is presented. 



Let ^t-1> ^t-2 *** ^t-n chronolo 

gical observations on present and past electricity demand 
at time t, t-1, t-2, ... t-n respectively, at discrete 
equispaced intervals of time. Derating (1) as the 
forecast for the time period t + 1, made at origin t, for 


lead tiiie 1 and as the actual demand for period 

t+l, our objective is to obtain a forecast function vdiich 


is such that the mean square of deviations (Z^^^ - (1)) 

between the actual and forecasted values is as small as 


possible for each lead time 1. In addition to calculating 
the best forecasts we specify their accuracy so that risks 
associated with decisions based on these forecasts can be 


evaluated. The accuracy of the forecasts are expressed 


by calculating the probability limits on either side of 
the forecasts for specified levels of confidence. 
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5.1 ASSUJIPTIOI® 

The following ass\i 3 nptions are made for formuilating 
time -series models for electricity demand forecasting. 

1. The time-series of observations Z^^- ... 

constitute a stationary- stochastic process. This 
implies that the generating mechanism of the process 

is independent of time. The parameters and probability^ 
distributions of components of the tiKie series are 
stable and tirae invariant. This implies that no struc- 
tural clianges in the parameters and process has tal'en 
place over the time span of observations, 

2. If the observed process is non -stationary, then the 
non-stationarity is of the homogeneous category, and 
it is assumed that it can be converted to a stationary 
series by a finite linear transformation, 

3. The joint probability distribution associated \7ith m 

observations Z. , Z. , Z. ... Z. made sequentially 
^1 ^2 ^3 

at any set of equis paced intervals of time t^, t 2 •••'^32 
are the same as that associated with m observations 

2(t, + k) > + k) ••• + k) 

(t^ + k), (t 2 + k) ... (tjjj + k) for all valiaes of 

m and k. 


89 


5. The stochastic process is covariance stationary. TMs 
implies second order stationarity. Mathematically the 
second order stationarity can be expressed as 


I, for all t ^ T E = a constant 

II, covariance of Z^ and ^t+k (i;here k is the lag) is 
a function of 3: only and is independent of t. 


Cov 



Z 



I I 

E |_(Z^ -A.) -^)J 

® G^t. Gtj 

V k 


where is referred to as autocovariance at lag k 

and is a function of k only. 


6, The stochastic process genera-ting the time series can 
be represented as a linear aggregation of random shociis 
and has a finite number of parameters , 


7. The linear stochastic process is invertible. 


5,2 STEPS IIT BUILDING TBE-SEEIES MODELS 
FOR ELECTRICITY DEMA!® 

We shall follow the following tliree basic 
iterative steps that are involved in building time -series 
models as proposed by Carlson et, al (21^) and Box and 
Jenkins (h$), Tlie iterative approach is represented 
diagiummatically in figure 5»1 33he steps are: 


t 
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1 . Identifier' t ion 

Using the data and any otlier infornation concern- 
ing the process, rougli methods of identifying a sub-class 
of t]ie general class of models are developed. The anal3"- 
tical tods of correlation analysis and spcctreil anaIl3'’Eis 
are used to identify the nature rind order of tiie process, 
degree of differencing and form of the models which are 
tents. tively entertained for parsimonious (21 5) represen- 
tation of tiie stochastic process. In addition the identi- 
fication process also yields rough preliminary estimates 
of parameters , 

2, Estimation 

The tentatively entertained models (obtained at 
the identification stage) are fitted to the data and its 
pararaeters are estimated. Rough estimates of parameters 
obtained during identification stage are used as starting 
values in more refined iterative statistical methods for 
parame ter e s tima tion , 

3 * Validation 

The model to he finally adopted is validated hy 
diagnostic diedcs. These checks are applied with the 

J 

; ' objective of uncovering possible lade of fit and diag- 
li; nosing the cause. Residuals of the fitted model are 
I ! analysed to test model adequacy. If no lack of fit is 


xnaicated, the riodel is finallj'- selected. In bne eveix oi 
an inadecuate model the iterative cycle of identification 
estimation and checkins is repeated till a suitable repre 
sentation is found. 


I Postulate general 
: class of models 



Identify the models to I 
be tentatively enter- 
! tained. ! 


Estimate parameters for 
tentatively entertained 
! model by refined methods 


j Validate the models by 
I diagnostic checking. 


'the modelNs 
adequate 


Use model for 
forecasting 


Fig, • Stages in iterative approach to 

time-series model building . 
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We no\f describe iii brief the procedures for identificatiori, 
estimation and validation of the time-series models for 
electricity demand forecasting. The procedure folloued is 
essentially the same ;;hich is suggested by Box and 
Jenkins (45) . 


5,3 IBIWTIFICATIOIJ 07 THE TIIE-SE7tIES liODLLS 


Identification methods are procedures applied to 
the set of data to indicate the kind of representational 
model \/orthy of further investigation. 

The electricity demand time series can generally 
be considered to consist of a random component, and non- 
random components that include trend, cyclicity and 
persistence. The random components refer to pure random 
component of the process. Trend refers to tiie long term 
behaviour of the tine series and is a function of time alone 
Persistence refers to linicage or relationships existing 
between values at a given time i;ith earlier vaD.ues and mscf 
be due to internal or external dependence, Cji^clicity refers 
to the periodic components that repeats themselves at defi- 
nite periods. The final representation of the time-ser5.es 
in terms of seveial components is given by 



^trend ^cyde '** ^ar ^^t-1» 

^ma ^ Ct-1 > €t-2’ • * * 


•** 

€t-l^ * Ct 
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idiere and represent the autoregressive and moving- 
average components, representing internal and external 
correlation respectively. ^ is the nurc rrndoin component 
at time t, 

Refering to lox and Icrlcins (45), Anderson (216), 
Nelson (217)} the autoregressive process oi order (p) is 
denoted by AR (p) is given by 

0 (B) 

The moving average process of order q is denoted by 
IIA (q) and is gi-veii b 3 '- -the relation 

6 (B) = a^ 

An autoregressive integrated moving average process is 

denoted by ARIIIA (p, d, q) and is expressed as 
d 

0 (B) = 9 (B) a^ 

The specific a,im at the identification stage is 
to identify the trend, persistence and periodic components 
of the time series and obtain some idea ot the value of 
p, d, q needed for the general linear ARIiiA model, Tlie 
identification stage also yields initial estimates of 
parameters. Tentative models identified at this stage 
provides a starting point for application of more effi- 
cient methods. 


9 ^ 


5.^ STAGES II-I TI-IE IDE IWIFI CATION PROCEDURE 

To identify an appropriate sut-class of the 
general ARIMA model (4-5). 

t (B) = 9o ® 

The following steps are follo^Ksd. 

5 .4- . 1 Identification of Trend 

A graphical plot of the data indicates the exis- 
tence of a trend, and their characteristic. A suitable 
functional relationship betvTeen and t; wliich is generall 3 '- 
linear, exponential or poljTioiiiial is assumed. Ttie coeffi- 
cients of the functional relationship are obtained by 
multiple regression analysis techniques (218). Statistical 
significance of the coefficients at 95^ confidence level 
are tested by *t’ values. In case the trend is found to 
be significant, the trend component is subtracted from the 
data to obtain a trend free series. 

5.4- , 2 Identification Of Degree Of Differencing 

The trend free data may be non-stationary. This 
series is then ddfferenced as mans'- times as necessary to 
produce stationarity there by reducing the process under 
study to a stationary ARMA. process (4-5) • 
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The order of the resultant stationary' autoregross- 
ive-novinc average process (APMA.) has to he identified no\7. 
The principal tecliniqucs of identifying the degree of 
differencing, and order of the process is by/ the application 
of certain analytic.al tools of identification and estimation, 
to be described next. 

5.^.3 Analytical Procedures For Iden t ific etion 
And -dstimation Of Tirie-Series Ilodelr 

The periodic and persistence conjDonent of the feje 
series are identified by correlation and spectral analysis. 
The basic analytical tools for time -series are the auto- 
correlation function, partial autocorrelation function, a:id 
spectral density;- function. 

5 Correlation Analysis Of Time -Series 

For the present study the autocorrelations of tZic 
time series for different lag values have been calculated 
by the methods suggested by- Box and Jenlvins (45) j Anderson 
(216) and Jenliins and V/at-ts (209). The pouer spectra 
of the time series have been calculated by the methods 
suggested by Jerfleins and V/at-ts (209). To have a crude 
check whether the autocorrela'bion is effectively zero 
beyond a certain lag k, a formula given by Bartlett (218) 
has been used. The partial autocorrelation of Hie time 
series have been obtained by solving the Xule-Walker 
equations (219j 220), 
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^.5 IDS.ITIFICATiai 0? T.X /jaiA PHOCESG 

Having estinatecl ine valuer of tho PACF, and ACF 
for specified lag values, the PAC? a.KL aCI are plotted 
grapMcally as a function of the lags, deferring to the 
procedure suggested in (4^?) for identification of the degree 
of differencing for the tr-^nd-froe scries, use is made of 
the -Dro'iorties of stationar-' tine scries that if none of 
the roots of 'die stationary moder are close to the bounde.i^r 
of the unit circle the ACF quiclcly dies out. Tendency for 
the ACF not to die out quiclily indicates non-stationarity. 

The degree of differencing that is necessary to achieve 
statio2iarity lias been acliicved rhon the ACF of the difference 
series dies out fairly rapidly. As suggested in (45) the 
time series is differenced at the most tvo times and the 
first !I/4 (maximum 15 to 20) estimated autocorrelations a.re 
inspected. 

Having decided the value of d the general 
appearance of the ACF, PACF, and the spectrum of the 
appropriately differenced series are studied to provide 
dues about the presence of periodicities, and values of 
p, and q, (order of AH process and MA process respectively). 
In doing so the cliaracteristic behaviour of tlie theoretical 
ACF, PACF, and spectrum for AR, lA, and AR>iA process 
described in (209, 45, 216, 217) are utilised. 
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The riToclels that h?ve been tentatively’- identified 
for the electricity (fenano time-series based on the infor- 
mation provided by the estimated ACF and PACF are presented 
in Chapter VI. 


5.6 IIIITI/iL SSTIIATES OF PARAI^EronS 


Tile refined iterative techniques for best estima- 
tion of para;-.x?tcrs require initial estimates as starting 
poin-ts. These initial estimates of parameters are deter- 
mined from the estimates of ACF and PACF obtained at the 
identification stage. For obtaining initial estimates of 
A .Pi, parameters the Yule-VJalker (219j 220) equations are 
solved. For sol-ving the Yule-Walker equations -tlie theore- 
tical autocorrelations are replaced by estimated auto- 
correlations, 

.deferring to Box and Jenlcins (1+5) ^Te \\rrite 
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r ai« the samixLe autocorrelatiors . 


0 


The solution of parameters J25's in terms of sample auto- 

correla.tions are given by 

/s -1 

0 = K . r 

P 

The initial estimates of moving average parameters are 
obtained, by solving the follov/ing set of nonlinear q 
sin-ultaneous eouations (45). 


®1 ® k +1 + *'* 


+ V • Qr 

0 “K C 


, 2 2 

(1 + e^*- + O 2 + . 




k = 1, 2, 


I’or a mixed AH!-iA. process, Box and Jenlrins (45) have sugges 
tod a computational algorithm i.'hich can be used for obtain 
ing initial estimates of parameters of IJRl'A process. 


The results for the initial estimation of para- 
meters are provided in Chapter VI, 


5.7 MAXILUM LIIOELIHOOD ESTI1IA.TES Of PATiAMHlTERS 

The identification process having led to a tenta 
tive formulation of the model, idle effi''ient estimates of 
parameters arc obtained by a refined iterative technique 
referred to as the maximum likelihood technique. Details 
of the procedure are presented in Box and Jenkins (45). 

An exhaustive exposition of the likelihood principle has 
been presented by varioi.:^ authors such as Fisher (221), 


Barnard (222), Tirnbaun (223). In the follovri.ng paragraphs 
MB give a brief description of the likelihood function. 


Let us associate an N dimensional random variable 

to the soi.-.plo of II observations Z-, , ... 2^ lAiose lino\m 

1^2 n 

probability distribution p ( Z depends on unlmovrn 

parameters , 

where, 


Z = vector of observp.tionc 

vector of parameters ju, 0 and of size P + q + 1 

For a fixed , before the data are available, p (Z \XJ 
wi3J. associate a density with each outcome of Z, After 
the data are available , ve are led to conternlate the various 


possible values of 




giving rise to tlie set of observa.- 


tions 2, The function for the purpose L ^ called 

the likelihood function. The likelihood principle states 
that all the data has to tell us about pareneters are con- 
tainedi in the likelihood function, Likeliliood function 
is generally unimodal for moderately large and large sample 
of data, and can be approximated by a quadratic function 
near the ma:ximum. Values of parameters maximising the 
likelihood function are called maximum likelihood estim8,tes 


(MLE), For further details on the maximum likelihood 
metliod, the reader is referred to Ba.o (213) » Whittle (2!^), 
Durbin (225) i Bartlett (1^3)* Barnard et, al (226), Kendall 
and Stewart (227), Mann and Wald ( 228 ), Hannan (135)* 
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5 *7 • 1 LiI;eliIiood Function For AlillJl Process 

Let ^-d+2 *** “o’ '^1 *** 

trend and cycle free tiKO scries vritli II = n + d observa- 
tions. Assuniin^ that the atovc series is generated hj 
an ALIIIA. ('p, d, c) process, •«>r generate a series ^ 2 , 

, , , Vn o- n = II - d observations, 

d 

^7here = /'V'H 

Fitting ARIi'IA (p, d, q) to the original series (trend and 
cycle free) is equivilent to fitting ARl-A (p, q) to the 
\i series , We can express the residuals a^ as 

H ~ '^'^t “ -^1 ''^t-l ” ^2 '^t-a • * • “ -^p ^t-p ®1 ^'t-l * * • 

%. • t = 1, 2, ... n 

Iloruially "vri-th d/ 0, the mean of \r^ series is zero, 

__ n 

Otheriase . w=< ^ w ) / n can be substituted for the 

t=r 

mean of w series, and w^'s have to be replaced by (w^ - u) . 

It is assumed tliat the a^’s are normally distri- 
buted. The 3 oint probability distribution of ••• 
is given by 

, ^ ^ / 2 tfa^ > 

p (a^ ... a^) ® 

Referring to (^5) unconditior^ likelihood for the 
general ABIMA model \;ith IT = n + d is given 
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1 (0, e, <) 


= f (0, ©) 



The imconditional STij:i of squares function B (0, 0) is 
given by 


S (0, 0) 



\Aiere 

^ ! 0, 0, w I = E \ a^j 0, 6, v 

L ' t ’ 

For moderately large n, f (0, 0) is insignificant and the 
log likelihood function is dominated by S (0, e) . Hence 
contours of unconditional suiii of squares function in 
parameter space i f6, 0) are very nearly contours of the 
likelihood function. Parameters estimated by minimising 


sun of squares function provide close approximation to 
maximum likelihood estimates. 


To calculate the unconditional sum of squares 
the a^*s are computed recurs ivel 3 ’-. A prelimineiy back calcu- 

r ■ 

lation provides values for I 3 = 1> 2s... 

to start the forward recursion. procedure described 

in Box and Jenldns (4^) supplies ts the unconditional 
sum of squares to any desired degree of approximation. 


Referring to C^5) the backwaid model 


0 (F) w 


t = 0 (F) 





generates the back forecasts 
In practice because of l^ie stationary 


1D2 


character of the AE operator 
some point t = - Q^with 0 
essentially eqtirl to zero. 


, estimates of "beyond 

of moderdte size becomes 


With sufficient approximation ue • rite 

-1 

(E) e (E) a^ 

''i>o 



In general the following eouations are used to generate 
bad: forecasts, and a^ values. 

0 (F) ~ ® f®t^ generates bad: forecasts 

0 (L) ~ ® 


(B) 


W 


generates the a4. values 


5,8 mXlWK LIKELIHOOD EXTIIIATES OF 

PARAl^lETERS BI NOhUIJEAE ESTIMATION 
ALGORITBM 


Maximum likeliiiood estimates are closely approxL- 
mated by the least sqiiares estimates which mal:e 


s i0, e) = 


n 


t=- 


- v2 


0, 0, V : a minimum 


In practice t^s is replaced by a manageable finite- stje 


n 
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Consider?ble sinplification occurs in the minimisation of 
the sum of squares function mth respect to j if eadi 
f^ ( /3 ) , (t = 1, ... n) is a linear function of the para- 
meters . Box and Jenliins (45) have showed that the 
linearity status of /B is different for AR, and I-IA process. 


Ho\:evcr tlie parameters are estimated in lineari 
sing the model for a^. 

Let Bo ~ ^ ^ 1 , 0 * &,o **’ ^ ) represent the guessed 

sot of parameters. Expanding | Q-t i Taylor's series 
we have 


k 



where 





ITow if X denotes tlie- (n + Q) x R matrix x^^, the (n + Q) 
equations may he expressed as 


r%]= ^ - w 

The adjustment ( jB-': which minimises s (^) = s i0, e) 

’a'j may now be obtained by linear least square, 
i.e. by regressing the ^ 'to the x*s. The adjus- 

ted values are siibstituted as new guesses and the process 
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pure random scries (viiite noise). Testing goodness of 
fit, by scrutiny of the residuals are described in 
Ansco:ibe and Tulcey ( 231 ) and the methods of normal plotting 
by Daniel ( 232 ). Box and Jenkins (*+5) suggest a procedure 
of checking the models by over fitting. 

For the purpose of the present study the model 
adequacies arc exarr.ined by tliree tests. Tl^ey are 

1 • Correlation anal 3 ''Sis of residuals 

2 . tost (chi -Square tests) 

3. Spectral anal 3 ’'Sis of residuals. 

V/e describe them in brief in the follo\'d-ng x^ara- 

graphs. 


5«9*1 Correlation Analysis 

Referring to Fox and Jenkins (45) the residuals 
are expressed as 

= 9^ (B) t (B) 

It is possible to prove that if the model is adequate, 
a^ become close to white noise as length of senes in- 
crease. The autocorrelations of are calculated by 
the method described in (45). Anderson has proved that 
with the knovledge of true parameter values J2S, and 0 the 
estimated autocorrelation rj^ (a^.) for t = 1, 2, ...n 
and k = 1, 2, ••• k are uncorrelated and distributed 
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approximately normolly with mean zero, and variance 1/n. 

Hence if the calculated autocorrelation lie within 95/^ 
confidence limits the model is considered to he adequate. 

5.9*2 7l^ Test Of Goodness Of Fit 

Durhin (233) has pointed out that it is dan, -porous 
to assess the statistical significance of apparent dcscre- 
pancies of estimated autocorrelations , fra.: their theore- 
tical values based on the standard errors which are actua3.l3' 
appropriate for theoretical autocorrelations of the residuals . 
Box and Pierce (2^) subsequently derived the large sample 
vare.ance and covariances of all the autocorrelations of 
for any APil4A. T)rocess, It was proved by then that use of 
1/n^ as standard error for rj^ v/ould under estimate the 

statistical significance of apparent departures from zero 
of autocorrelations at low lags but could be employed for 
moderately largo and high lags and n should be treated as 
an upper bound for the standard errors of r^^ ^^‘s) . 

To circumvent the above problem, instead of 
considering the individual autocorrelations the first t\jenty 
or so autocorrelations of a.j.*s are taken as a whole to 
indicate model adequacy. I>ue to a result by Box and 
Pierce ( 23 V), it is possibLe to show that if the fitted 
model is appropriate 
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K 

« = " Z 

]c=1 




(irficrc ) are the first K autocorrelation of 

residuals of the ARIIIA (p, d, q) process , is approxiriiately 

■y 2 

distributed as /\^Q-_p_q) vherc n = II - d is the number of 

v’s fitted to the hiodel and k -- p - q = degree of freedom. 

Test of hypothesis of model adequacy is made by comparing 

the calculated Q value vTith theV ^ ... ^ values wliidi are 

/V. critical 

obtained from a talle of percentage points of^X-^. If 

Q critic al given level of significance, the 

model is assmned to bo adequate. 


5.9.3 Snectral Analysis Of The Residual Series 

Jenliins and Vlatts (209) have proved that for a 
pure random series in the frequency range 0 f ^f^, "the 
spectral density function is constant and is equal to the 
variance of the process , 


The number of degrees of freedom for the spectra is given 



■vdiere 




n 

f 



and n 
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The sample spectral estimates arc distributed about the 

population spectrum according to ( distribution. 

From tables of ( /'r'> distribution the confidence 

limits for a specified level of confidence (1 can 

be obtained ^other\n.se fren distribution, tables 

T" obtained. Dividing these 

values by "1 gives the aunropriatc li;.iits for X \ . 

\ ^ 2 ' 

As Yf) 3- A- A — distribution, the series can be 

considered pinre rando.m at 100x (1 percentage confi- 


dence level if 

71^/2. T- 


r* 


^TTfl 




2^ 


'T 


where 


G ’(f) = spectral density of the pure random series, 


Gj^ (f) = the spectral density of a^ series. 

If the spectra of tiic scries lies \;ithin the above 

100 (1 - Aj) % confidence limits, then series is 

considered to be random at the given confidence level. 


Results of diagnostic diecks for testing model 
adequacy have been presented in Qiapter VI. 
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5.10 FORECASTING TI!E TUTUEE VALUES OF PEAK 
POVJER AND ENERGY DElIA!® BASED Oi: Tlttl 
THE SERIES hODELS ADOPTED 

Having selected the models of the time series of 
electricity demands the forecasts of future values of 
electricity demand, the variants of theses forecasts and 
their probability limits for specified levels of confidence 
are calculated by the procedure described by Eox and 
Jenlcins (45) . 


CHAFEER VI 


RESULTS, DISCUSSIOI'B AMD COMCLUSIONS 

In C3iapter IV econometric models -were developed 
for forecasting the future demand for peak power and energy, 
and for analysis of electricity demand in the industrial 
sector. Stochastic time-series forecasting models were 
formulated in (Chapter V for determination of forecasts of 
future demand for electricity. This chapter is devoted to 
the presentation of the numerical results obtained by the 
application of the above models to actual case studies. In 
Sections 6.1 and 6,2 results for mcro-level forecasts of 
peak demand and energy for India for the period 1976 - 2000 
are presented. Section 6,3 deals with the results of the 
analysis carried out for industrial electricity demand. 

The industries considered are Steel, Textiles, Paper, 
Chemicals, Non-ferrous metals. Food, Vehicles and Mining 
and Quarrying. Secticai 6.4 presents the numerical results 
of forecasts of peak povjer and energy demand for India 
for the time horizon of 1976 to 2000 A.D, by the applica- 
tion of stochastic time -series models. 
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6,1 CQCFi'Icmiii ESTIl'iA.TES OP PCC:iOlCSG?RIC MODELS 
fceec/^timc peai: po\jeh aid eierui rniiDs 

The functional relationships presentee! in 
Chapter ^ (equations ^*5} !^.6 ) for * 

energy and peak deioand i^cre estin.ated by orc3inai^' least 
square technique. In the first stage all the explanatory 
variables were incorporated in the equations to 1*® es ui- 
mated. Tlie simple correlation matrix indicated significsait 
collinearity bet'ween all the indenendent varialoles, indi- 


catin^i the existence of imiLticollinearity . Because of 
this, it was inferred tliat most of the parameter estimates 
were hi^y biased, and insignificant. As edl the explana- 
tory variables are significantly time -trended, inclusion 
of time as an independent variable vm.th equations was ruled 
out, to avoid the problem of multicollinearity. Various 
combination of variables for each of the simple » double 
log, first differenced, and first differences of logged_ 
variables form of equations \rere tried, satis- 

factory' fora of equations on tlie basis of botli economic 
theory and statistical inference in its estimated fora for 
forecasting energy and peak demand were found to be ox ohe 


following. 


Model 1 

E. = -6183.96 + 200.99 GITP. + p . = 0.9777 DW = 1.87o 

^ (7.07) (^.97) ^ 
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Model 2 


\ = -21336.23 + 399.22 I. + 

^ (10.99) (26.22) ^ ^ 

2 


H = 0.9372 

DW = 2.235 

I'bdel 3 



E = -79OOI+.7 + 254.58 P. + f:.. 

(12.67) (16.7Q) ^ 

2 

F. = 0.8580 

LW = 1.579 

I-iOdel 4 




_2 

Log E. = 3.3M+ + 1.301 Log Gl’IP. +r-. R = 0.9786 DW = 1.866 

(25.97^) (45.90^) ^ 


?fodel ^ 

_2 

Log E. = -24.1+72 + 5.61+0 Log P. + ^. R = O.9965 DW = 2.02 

^ (83.857) (115.09) t v^t 


Model 6 

Log E. = -1 .709 + 2.397 Log I. + A . 

^ (5.69) (36.255) ^ 



0.9662 


DW = 2.65 


where 

E^ = Electric energy denanded in period t, 

GNP^ = Gross national product (indices) in period t, 
DW = Burhin - Watson statistic. 


Models Adopted For PeaJc Demand 


Model 1 

PD. = -638.681 + 31.299 S = 0.9878 DW = 1*932 

^ (6.361) (60.923) ^ 
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Model 2 

2 

PB. = -3130.17 + 63.2^9 If + R = 0.9773 BW = 2.189 

^ (17.180) (4lf.^62) ^ 

I-bdel 3 

-2 

PD. = -12209.00 + 40.225 P4- + £+ = 0.8931 = 1.965 

^ (14.57) (19.633) ^ ^ 

Model 4 

_2 

Log PB. = 2.642 + 1.102 Log G!IP. +(£V R = 0.9586 BW= 2.054 

^ (17.23) (32.68) ^ 

Model 5 

_2 

Log PD., = -21 .09 + ^.807 Log Pt + Pt = 0.9885 IM = 2.73 
* (46.29) (62.81) 

Model 6 

2 

Log PB+ = -1.857 + 2.079 Bog + et ^ = 0.9933 Bl\- = 1.793 
^ (16.216) (82.39) 

where 

pD =: Peak demand j^in MW) in period t. 

W 

For all the above equations the figures in the parentnesis 
under the coefficients are their corresponding *t* values. 

is the coefficient of deternrinatian adjusted for degrees 
of freedom. BW is the Burhin-Watson statistic. 

The coefficient of determination for each 
equation is significantly high; there hy indicating that 
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these structural equations are satisfactory and are adequate 
representation of electricity demand in the past. Gross 
national product, population, and industrial output are 
appropriate explanatory’- variatles for representing past 
demand for electricity. Further the »t* values for the 
coefficients estimated are quite high indicating that the 
explanatory variables are significant in explaining the 
causality hetyeen the dependent and independent variables. 

The values of the D.W, statistic -signify that the hypothesis of 
no au-tocorrelation in error tenns cannot be rejected at the 
95^ confidence level. 

The coefficients for the linear functional form 
of equations can be interpreted in the follov/ing manner, 

A unit increase in the iMex of GHP using 1961-62 as the 
base y’-ear vriLll raise the demand for electrici-ty by 201 
million KWH, From the relationship linlcing ind-us trial 
output and electricity energy demand v:e find that a rise in 
industrial output leading to a unit rise iii the index of 
industrial output ■'.dll generate demand for an additional 
399.92 million KWH of electricity energy. Further, from 
the estimated relationship for population, the coefficient 
suggest that a million increase in population will lead 
to an increase of demand for enei^ by 25^.58 million I'iWH. 
Ccanparing the R val-ues of the simple linear relationships 
per energy’- demand it is observed that population explains 
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about 85«B^ of tlio total variance of demand for electricity 
energy vhere as the corresponding figures for CNP and 
industrial output are 97.77^ ".nd 93 .72^^ respectively. Tliis 
leads us to suggest that amongst the throe explanatoi^’’ 
vaitLables the GMP is a more appropriate causal variable, 
explaining the largest variance in electricity consumption 
followed by industrial outiDut and populatio;i in that order. 

From tlj.e results vre observe that electricity 
requirements for a unit increase in industrial output is 
almost double that of the requirements for an unit rise in 
GNP, Tills result is as expected, V/e provide the following 
explanation for this phenomenon. The available data on 
electricity consumption by various sectors of the econon^y 
indicate that, industrial sector accounts for about 
of the total consumption. Hence a unit rise in index of 
output may lead to such heavy requirements. Further the 
gross national product of country comprises of output from 
all sectors of tiie econonry, A31 the sectors are not such 
disproportionately heavy users of electricity as the 
industrial sector, compared to their output. For example 
the transportation sector mostly relies upon coal and oil 
as its major source of energy. Despite the gromng needs 
of the agricultural sector for energisation of tube wells 
and mechanisation, this sector still has a minor share 
in the total consunption of electricity. In addition the 
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contribution of industrial output to GUP is only about 60;.J 
rieaninc thereby that a unit increase in GUP has received 
about 0.6 of its value fror. t’-e industrial sector. From 
these considerations it is clear tliat a unit rise in GUP 

require lesser aniount of electricity energy as compared 
to the requirements for a unit rise in industrial output. 

The stTUctural coeffi events of the log linear 
form of equation represent the ei.asticities of demand 
corresponding to tjie variable under discussion. The income 
elasticity of demand for energy as v/cll as the population 
elasticity ano. industrial output elasticity are greater 
than unity as well as for population and industrial output. 
One percent rise in the GUP gives rise to 1.35^ rise in the 
demand for electricity energy. Similarly l/fo rise in indus- 
trial output and Increase in population gives rise 
2,'^% and increase in demand for electricity energy 

respectively. From the log linear form of equations also 
we see that the requirements of electricity energy for 
rise in industrial output is almost double that of tlie 
requirements for one percent rise in GNP. As in the 
sinjiLe linear form of equations (equations ^.3 and h.h of 
chapter h) tiie variables have maintained tteir relative 
positions of explanatory ability. 

The results of the models for peal: demand can 
be interpreted in an exactly similar manner. Here also 
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the results iiic'.icate that GfJP is the r^st apuropi^atc 
explanator;' variable exelaininr the Irrj^est variance in 
electricity consuiiption follo\red by iiiclustrial output and 
ponulation. A unit diange in in.'e:: of CUP, industrial 
output and population gives rise to aiv ac'ditional peak 
demand of 31 •29, 63.2^, and ^0.22 ilJ respectively. Simi- 
larly one percent increase in GnP, industrial output and 
population give rise to a 1.10^, 2.08,^;, h,o0% increase in 
peal-; pov/er demand respectively. 

The kno\7ledgo of the elasticity coefficients are 
useful for policy decisions. For example, v^ien plans are 
foimulated vith an objective to increase the GI'IP by a 
certain amount, or the industrial output by a certain 
amount, siii table investments have to be made in the power 
sector to meet this additional demand as a consequence of 
increase in output. Similarly appropriate policy measures 
have to be taJ;en in terms of determining tlie optimal mix, 
location, and time phasing of these additional capacity 
that has to be added to meet the demand. 

6.2 FORECASTS OF FUTURE DStAISDS FCE PEAK 
POVffiR ARD ENERGY. FOR INDIA 

The forecasts of peal; pov^er and energy demand 
have been obtained for India for iiie yeais 1976 to 2000, 
on the basis of the estimated relationships. The future 
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values of GiiP, populate on, and Indus t‘>'; "'1 output have been 
estimated by applyin^^. planned and anticipated growth rates , 
The plaraicd and aziticipatcd growth rates have been talven 
from a publication of the Planning Commission, Goveriraent 
of India (235 )• Substituting these values of the indepen- 
dent variables in the estimated demand equation \fe obtain 
point value of predictions. Interval 2 )redictions for 95f 
level of confidence, (for t’ne individual values of fore- 
casts) have also been presented, IIuiTierical results have 
been presented for both peal: power and energy demand for 
twentyfivo years in future for the various groirth rates, 
Grovrth rates considered are - 5%, and 7% for GNP; 8^?, 

9^ and 10^ for industrial output; 1,5^» 2.0^ and 2 , 25 ^ 
for pop\3lation, Bie anticipated and planned gro’.rth rates 
according to the Planning Commission (235 ) for GilP, 
industrial output and population are 5^> and 1,5^ 
respectively. 

Tables '1'. to 3 present results for energy 
demand and Tables ^ to 6 present results for peal: 
demand forecasts obtained from the simple linear foim 
of equations (equations ^.3 and 4,4), Table: >2 to 
'Sj present results for energy demand and Tables 10 
to At', present results for peal: demand forecasts ob- 
tained from log linear form of equations (equations 
4.5 and 4.6), 
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Forecasts of Energy Demand From Gross ifctionol Product (GNP) Pyojpc-Vto'rvs 
(Log linear form of equations) 

Level of Confidence - 9^^ All forecast values are log transformed values (;o base 
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obtained from population projections are likely to be. less 
reliable. The forecast obtained from GUP and industrial 
output do not differ videly. However, the GNP based 
forecasts are lower in magnitude as compared to industrial 
production based forecasts, Tliis can be attributed to the 
fact that GNP includes output from all sectors of the 
economy and all these sectors are not heavj'- consumers of 
electricity in comparison consumption in the industrial 
sector, A five percent growth in GNP does not correspond 
to eight percent growth in industrial output. 

The 95!^ probability limit forecasts provide us 
quantitative information on risks associated mth plans 
based on these forecasts, Tliis is one of the chief advan- 
tages of adding probability dimensions to forecasting. 

The additional requirement of electricity energy over that 
of 1974 - 75. for the teiminal year (1979 - 80) of fifth 
plan of India is expected to be approximately 25300 Million 
KWH and for the terminal year of 6th plan about 60225 
million ICV/H, The additional peak power demand based on % 
rate of growth of GNP for terminal year of fifth plan and 
sixth plan ai*e 4704, and 10195 M.W, respectively. This 
implies that these additional amounts of capacities have 
to be installed for meeting peak demands 
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6.3 KCSULTS for ECONOJIBTRIC AFALYSIS op 

FOR ELECTRICIiT III THE INDUSTRIAL SECTOR* 

In this section the results of the Gconomotric 
analysis of cer.iand for electricity in the industrial sector 
liave been presented. The analyses has hcen carried out 
for the follovang industries, naroly Food, Textiles, Iron 
and Steel, Clienical, hon-ferro'iis netals. Engineering, 
Vehicles, Paper, Mining and Quarrj’-ing, 

The ordinaiF^ least squares techniques has been 
used for estimating the econoLictric equations. The results 
of the analysis provide us quantitative estimates of elas- 
ticities and lag terms , Further, the results provide us 
\at3:i information for carrying out the following kinds of 
comparisons . 

1. Direct comparison of goodness of fit of particular foiras 
of equations for a particular industry group. To the 
extent that these equations can be considered as alter- 
natives , we can draw conclusions on the relative merits 
of each . 

2. Inter-industry comparison of size and significance of 
parameters and goodness of fit of equations. This 
type of analysis could be follox^ed up by a cross 
section analysis by relating the inter-industry 
differences to economic factors . An examination can 
be made of as to ho'-j differences underlying economic 



en^ironiTiGnt in wliich industries o'^ercite condition the 
responsiveness of their electricity demand to c} 7 an::es 
in certain causal variables over time. 

Besides these two types of coj^parisons , the 
effects of inclusion of particifLir varirhles on the general 
fit of equations and on the size and significance of other 
paroj’eters have 'boe:-! studied. 

The results obtained are interpreted using tlie 
above stated basis of comparisons. The interpretation of 
the results are presented in a later section, 

6*3.1 Evaluation Of Goodness Of Fit Of Equations 

The procedure that has been followed in comparing 

the goodness of fit of equations is to calculate for each 

2 

equation (over all nine indust ly groups) the mean R 
(1-fultiple Correlation Coefficient) and the mean percentage 
of standard error (S.E) of equations together with the 
corresponding standard deviation of tliese means. The equa- 
tions are ranked on the basis of (a) highest R , (b) lowest 
standard error of equation Cc) lo^/est standard deviation 
of R^, (d) lowest standard deviation of the standard error 
of equation. 

The degrees of freedom vary between equaticxis. 
Hence it is possible for R^ and S.E of equationsto give 
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conflicting results. The extent of arnhiguity can be esti- 
mated. by acte rmiuDJig tne ranli correlation coefficient 
betv^een ranlting on the basis of (a) and (b). It is founrA 
that this coefficient is 0.925 and therefore there is little 
significant ambiguity . Comparison of : 7 ith its S.D.j 

and S »E , of eouatior \lth ?.ts S.l',, indicates the extent 
to \diicli equations perfor;' conci.stently wei.l or poorly over 
all industries. In fact li: io found that those equations 
which perform best perform consistently best having the 

smallest standard de\riaticns. Tne ranh correlation coeffi- 

2 

cient between E and its S.D.; S.E, of equation and its 
S.Di will give an idea bet-ween the relation between rela- 
tive consistency and relative goodness of fit. Again we 
find reasonably high correlations. Let p denote the ranli 
correlation coefficient. The follo\7ing results have been 
obtained for 9 • 

^ (R^, S .E of equation) = 0.925 

f (R^, S.D. of (R^)) = 0.9^-6 

j>(S.E.; S.D. of (S.E)) = 0.729 

An examination of Table 1 3 indicates that the equation 
ranked highest on R^ is also ranlied highest on S,E, 
of equation and it also tends to have lo^rest S.D, imply- 
ing greater consistency ov.-r:all industiO" groups,. 

■ Although there is some smtehing in position as t3ie 
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TABLE 13 

— oE goodness of fit of Gouatioi'is for* 
industrial electrl city”~doniand. 


Ho. of { 
the ( 

Equa- 5 
tion .J 

Mean y Rank 
•n2 5 liO. 

^ 5 

- 1-- 

5” iTean 
iS.E of 
5laua- 
.gtlonC^: 

1 Rani: 

0 Ho. 

)5 

1 S.D. 

0 of 

y p2 

5 

r Rani: 
5 He. 

5 

5 

I S.D. 

0 of 

5 S.E. 

0 

f Rani: 

0 No. 

0 

0 

4.12 

0.981 

1 

5.12 


2 

o.c 5 ^ 


2.01 

2 

4.13 

0.979 

2 

4.53 


1 

0.C48 

1 

1.79 

1 

4.16 

0.975 

3 

5.79 


4 

0.068 

3 

2.46 

‘ 5 

4.11 

0.974 

4 

5«46 


3 

0.079 

4 

2.23 

4 

4.18 

0.970 

5 

6.07 


5 

0.0G3 

a 

2.67 

7 

4.15 

0.968 

6 

6.24 


6 

0.094 

6 

3*16 

10 

4.17 

0.963 

7 

6.43 


7 

0.099 

7 

2.58 

6 

4.19 

0.953 

8 

6.96 


9 

0.083 

5 

2.13 

3 

4.25 

0.949 

9 

7.34 


11 

0.113 

8 

2.84 

8 

4.23 

0.944 

10 

7*01 


10 

0.120 

9 

2.99 

9 

4.24 

0.938 

11 

7.34 


11 

0.134 

10 

2.84 

8 

4.26 

0.902 

12 

6,78 


8 

0.148 

11 

2.46 

5 


Explanation of symbols : 

2 

R - Coefficient of determination of equation, 
S .E - Standard error. 

S,D, - Standard deviation. 
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criteria are changed, yet results "bet'/een grouns of indus- 
tries are more stable, £nuations(^,12 c: 4,l3)are marginally 
tetter tnan er[Uations(4-,‘l^c.o(^,^3)on al3. the foup criteria. 
Hov'ever, the differences are very small. In %lev 7 of the 
fact that L-ie record group is suoject to iioasurement error, 
the results of the record group are good. The inclusion 
of a time variaolo in equation(^,16)may possibly account for 
goodness of fit there , ■’jhile in the rest of this groun the 
presence of the coal variable has a considerable effect on 
goodness of fit. The t’lrd group of equations correspond- 
ing to Model 3 are clearly consistently inferior on the 
criteria selected, 

'* Hoirever the narrovmess of this basis of comparison 
and smallness of differences, especially between group of 
equations corresponding to I'bdel 1 and Model 2 do not permit 
a firm conclusion \diether one set of equations (and there- 
fore the hypothesis involved in that model) is superior to 
another, for purposes of forecasting and analysis. Similarly 
lid thin the groups of equations differences are even smaller 
and comparison even less conclusive. 

6 . 3.2 Explanatory Power Of Variables 

The explanatoiy power of an equation included 
is determined by explanatory power of the variables and 
the form of relationship that it is assumed to take. 
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2 

Though R is an important criterion Tor toctinj tlie good- 
ness oj. rit ojid in assesing tlie rel^.uive usei'dlness o_’ 
any set of equations, in tliis an-13'sir t'lis sIjOuIcI not be 
the only criterion ul'dch shorJLd be used. ? e couse of in- 
clusion of a lagrcd dependent variable it is possible for 
an equation to have hifli R^ but ’'it]i none of the meaning- 
ful explanatory variables significant. This may be either 
due to lack of exnlanator}" pover of variables concerned or 
due to presence of multicollinoarity , which is then a 
wrongly specified equation. Hence we adopt tiro further 
bases of compo.rison. These correspond to (1) the gejieral 
level of significance of variables in an equation, and 
(2) tlie degree to which multicollinearity appears to be 
present. 

Clearly it would require a large amount of space 
to present the whole set of parameter estiicates for each 
variable in eadi equation in each industry group. There- 
fore, we present a summarised picture of it in Table "iPt. 
Tliis table gives for each variable in eadi equation tiie 
number of industry groups in wiiich the parameter estimated 
was significant at the 0.95 level. From t!ie table we get 
an idea of the relative explanatory power of particular 
variable. From Table 1^it is clear that apart from 
lagged dependent variable, the variables which are the 
most common determinants of changes in electricity 
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of e^mlamtor^r T30wef» of variables for industrial 
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6.3.3 Inter-Indu::trv Elastlcitv 
Coefficicnts Cor,m:irison 

For inter-industry coranaris ons , ve select froii 
the large nu’-’her of vnrj'ing estrhiates of coof.fici6nts of 
a particular variable, the one unich is considei’od best 
for a particular industry* group on the criterion given 
above. 

Table 1-5 gives soi.e selected resiiLts from the 
first tv70 groii'.'S of equations. It can be seen from table 
9 that tlie industries fall broadly into tvo categories. 

Those industries for \7hich significant index of production 
elasticities are greater than unity. The industries wliic’i 
belong to tiiis categor:/ arc Food, Nonferrous metals, 
Textiles, Iron and Steel, Vehicles and Mining and Quarrying, 

None of these industries have significant fuel 
price elasticities . The price wage relative is Important 
in three classes namely Textiles, Nonferrous metals. Iron 
and Steel. The remaining industries do not sliow ary res- 
ponsiveness of electricity demand to changes in price 
relative to wage rates . 

It is noticeable that food, and nonferrous metal 
industries have non— si gnificant coal variables • This 
suggests absence of developments affecting the pattern of 
fuel use such as have oc cured in other industries. In 



TABLE 15 

Estimated elasticitj^- coefficients of industrial electricity demand. 
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consiuiiption arc the index of prcxiuction, and Coal and 
Colvc consur.iption • ijic oinic variable "ii equation is 

si£piificnnt in surprising:!}’- foi; cases and tliis may be due 
to int-?rcor relation uith and C^, both oi vhich are 
strong].y time trended . Anart from perhaps the (P/F) and 
(^/Q) variable in cl:e Lodcl 3 equations, there are no 
other variables i.lt]i fairly ccneral explanatory pover. 

The follo\;ing conclusions can be dra^m frOij the result 
presented in the table. 

(a) As between groups of equation there is lit-Qe to 
choose between first group (Model 1) and second 
group (Model 2). Each group contains one very good 
equation, the remainder being fairly good or indi- 
fferent, It therefore api^ears that use of linear 
form and measurement of electricity in coal equiva- 
lent tons (C.E.T) makes less difference than what 

WG had opined before. Both these groups are markedly 
superior to the third, partly because this third 
group relies on gcncra.l!}'- less successful variables 
and also because of the form assumed, 

(b) In comparing equation between groups there is in eadx 
group one eq-uation which is markedly better than 
others. In first group this is equation (^#12)* In 
second group equation(*t-,i^J*^s -tlie highest general 
level of significance. 



the absence of any other significant determinants it would 
seem that in these two industries capital stock effect of 
groirth in industrial output accounts for the In.gher magni- 
tude of electricity consurntion rather than substitution. 

Those industries whicli have significs.nt index 
of production elasticities less than unity. The indus- 
tries ’./liich I' el one to this group are Chemicals, Engineering 
and Baper. 

The fuel nrice relative is not significant in 
case of Ghenical, Bp.per and Engineering industries. For 
chemicals alone there is a high and significant elasticity 
of the (JV'Q) variable suggesting a labour substitution 
effect, which in the presence of a non-significant (P/W) 
variable is probably the result of technological develop- 
ment. For Engineering and fhper industries the significance 
of coal elasticity is difficult to interpret. This may be 
in part due to substitution as a result of relative price 
movements. For these industries, the results do not show 
any discrimination between price and technological subs- 
titution effects. 

Finally we consider the results of the consump- 
tion ratio form in Table A6 . The general lack of signi- 
ficance of variables in this case reflects the inappro- 
priateness of the assumption that elasticity of ^ectricity 
demand with respect to output is unity. the other 



TABIX 16 


E.os'iil'ts o f* elasticity coemcionts obtained fron 
Model - ^ Bruat i ons 


Industry’’ 

Group 

1 ^D/Q,P/F 

|®D/I5,P/U p 

D/Q,1VQ 

1 Vq,!/? 

Food 

-1 .917 

-1 .456 

-2.206 

1.094 

Chemicals 

-0.587 

- 0 . 113 * 

- 0 . 071 * 

0.01^!-* 

lion ferrous 
Metals 

-0.887* 

-1 .044* 

0 . 298 * 

-0.062* 

Iron c": Steel 

-1 .5ih1* 

- 0 . 690 * 

0 . 179 * 

0,616* 

Engineering 

-0.i^16 

- 0.319 

- 0 . 149 * 

- 0 . 181 * 

Veliicles 

-1 .V 9 O 

- 1 .389 

0,468* 

G.OOo* 

Textiles 

-1 .425 

- 1 . 090 * 

0.462* 

0 . 195 * 

Paper 

-0.401* 

-0.029* 

0 . 227 * 

0.126* 

Mining & 
jQuarrying 

- 1 . 059 * 

- 0 . 518 * 

-0.61 2* 

1 . 035 * 


* Indicates that these coefficients are insignificant 
at 0.9U level of significance. 

E-n/i^ •n/i:. - Elasticity of D/Q 'jith respect to P/E. 

D/Q,P/F 

Similar interpretations are to be made for other elasticities . 


vdiere in Table 15 the index of production elasticity is 
rel'.tively close to unity, the results on significance 
of variables in the first and third models, tend to coin- 
cide oarticularly for relative price variables. Never- 
theless the poor results of the model 3 vhich is based on 
consumption ra.tio suggest that the first tuo models 
constitute a better approach to electricity demnd analysis. 

6.3«^ Interpretation Of Results 

The main points that we ^dll discuss are the 
inter-industrj’’ difference in results and the results for 
the relative price variables. However, before discussing 
the resijlts we consider the various statistical reasons 
for cautious interpretation of the results. 

The use of distributed lags may give rise to 
autocorrelation with consequent bias in estimates of st^- 
dard errors and therefore it may result in mistaicen con- 
clusion on the significance of parameter estimates. The 
D-W (I>urbin - VJatson) statistic has indicated that auto- 
correlation is present. 

Further althou^ sane attempt was made in the 
formulation of equaticais to minimise the risk of inulti- 
coiTLinearity between independent variables, it evident 
from the simple correlation matrix of all variables that 
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TABLE 17 


The effect of correlation between Ya -ria'hTpg 
levels of sigitiflcarifS ■ 



Iron Cl Steel ~0,^6 - 0,769 7.6^ 2.9 -0.223 4,75 


Mning & 

Quariying 0.6?8 0.753 6.97 1.94 O .075 5.03 


Chemicals 0.858 -O .965 4.93 1.05 -1.823 3.83 


Non ferrous 

metals -0.795 -0.897 8 .I 9 2.07 -0.102 6.12 


Paper -0,820 -O. 9 O 6 3.97 1.758 -0.086 2.20 


The symbols P/F and PA have been explained in Table No. 14 








inrw 


The level of significance of variables and. la.bour 
intensity of jjidus tries arc tabiaated in Table - 1^'. 

There are clear differences beWcen labour- 
intensive and capital intensive industries in respect of 
the significance of output and coa!3.-consuniptiQn vario.bles . 
The output variable is significant in all the industries 
concerned, where as the coa3. va.riable is not significant 
in all industries. This mean that the relationship betireen 
output and electricity deiiand may to a large extent be 
deterjained by the capital -stock characteristics. The 
explanation of the pattern of significance of coal consump- 
tion variable lies in the teclmologies of the industries 
concerned. 

P 

The results indicate that ( g ^ variable is 
not significant in mining and quarying, vehicles, paper, 
food, and chemicals. This suggests tliat the possibility 
of substitution of electricity for labour in these indus- 
tries are necessarily circumscribed by the fact that 
electricity may often be embodied in expensive machinery. 
Although relative change^ in electricity price and wages 
ma; ?- to some extent lead to i-xjchanisation, it still may 
represent a proportionately small change in the relative 
costs of production. However some industries like textile, 
engineering, non-ferrous metals, aixl iron and steel do 
appear to have some price— rra.ge responsiveness to electricity 
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table 18 

Levels of significance and laLour intensity 


Industry Group 

"h- 

■J - ^ 

^ P/F 

Variabl 

>rT 

es 

I/Q 

fii/Q 


lining and 






— )L , 

Quarrying 

-f 

- 

- 

- 

- 

4 

Vehicles 


- 

- 

- 

- 


Textiles 


- 

+ 

- 

- 

4 

Engineering 


- 

+ 

- 

•f 

- 

Paper 


- 

- 

+ 

- 

4 

Non ferrous metals 


- 

+ 

- 

- 

- 

Iron Cc Steel 

+ 

- 

- 

4" 

4 


Food 


- 

- 

+ 

4 

- 

Qiemicals 


- 

- 

- 

4 

- 


+ Denotes that tliese variables are significant, 

- Denotes that these variables are non-significant. 

Hote ; The industries are ranked in order of increasing 
capital intensity. 

Symbols Q, P, F, W, I, M, C have been explained in Table 14, 
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demand. These conclusions are of course subject to 
reservations stemming from the resu 3 .ts tabulated in 
Tabic 17 . 

The fuel price relative is almost insignificant 
in dl the industry groups. This is expected. From the 
available data on the value of electricity costs per ■ 
thousand rupees of output, there does not appear to be 
any importance of electricity as a cost factor .and hence 
there is no 1-12 s pons iveness of demand to chaises in price. 
For industries vhere (P/F) variable is not signific.ant 
there are two other possible explanations . Die first is 
that there is a strong conrplimentarity between fuels and 
capital equipment. It is not possible to substitute among 
fuels 1-ri.thout changing type of plant and raachinery used, 

3 0 that other things being equal changes in relative fuel 
costs has little effect on relative costs of producing from 
different equipment. These relative costs can change in a 
vay that it may not talce advantage of favourable relative 
price movements. The second explanation is a statistical 
one. The fuel price relative is derived free general price 
index for fuels used by manufacturing industries and so is 
same for all industries, Althou^ this index is a good 
indicator of the general level of price paid, the actual 
price paid b3’’ firms vary according to area and size of 
firms and also between industries. There is a degree of 



approximation involved vhich nay imply that results do not 
measure actual response to price changes, Tliis second 
explanation raises the question of the reliability of 
significant estimates obtained. Can we infact infer a 
causal relationship or do the significant elasticities 
simply'" represent correlated time trends of tiro variables • 

To test this hypothesis equation (4- .14 ) vas estimated 
■with (P/F) replaced by time trend t. The resiilts for 
this equation for all industries shovred that n^, and index 
of production elasticities I'/ere roughly same as in equaticai 
(4- .IS) • simple correlation betv/een price and time "v/as 

0.94-1. Thus it is difficult to refute that the price 
variable has acted as a time trend. 

Finally as noted above there appear to exist a 
relation between capital intensity and significance of 
index of production. It is also noticable that capital 
intensive industries have non significant price effects, 
and have index of production elasticity greater than unity. 
This tends to confirm fairly closely to our hypothesis 
suggested in section 4-, 2.1 of Chapter IV, However, there 
is a qualification "which must be made on statistical 
grounds , The output for each greup of industry is a base 
weighted index for broad group of produces, I'jhile ^ectrl-" 
city consumption is sin 5 )ly unwei^ted total, dearly 
therefore it may be possible for equal propoiticsiate 
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growtli in electricity consumption and output to take place 
in each sub-industiy while the system of v;cighing for index 
of production could lead to a different relationship in 
aggregate, 

CONCLUSION’S 

The main conclusion of this analj’^sis is that the 
relative price changes are not important deterniinants of 
growth in industrial electricity consumption. The cliief 
deterninants are gro\7th in industrial output and changes in 
technology taken at a face value, the results for the 
relative price variable suggest that price elasticity of 
demand is highly insignificant. This, if valid, would seer; 
to have a relevance for current developments in the energy 
economy. It suggests that considerable changes in relative 
price may have to be necessary to offset even partially tlie 
gro\;th of industrial electricity demand and that dther form 
of energy a 3 ?e higlily unlikely to have any significant effect 
on electricity’s share in the industrial energy ccoisiimption. 

On the other hand the results na?'" have led us 
to an incorrect interpretation. The possibility thatthe 
price variable mostly acted as a time trend lias been 
demonstrated so that where the variable is significant 
ve cannot necessarily infer that there was a significant 
price effect. At the sane time we also cannot infer that 
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where the variable is not significant, there would not in 
fact he a signific^t price elasticity. This arises because 
the price variable has been derived from wliolesale price 
index of fuels used by manufacturing industries. It does 
not necessarily reflect the actual price paid by a particu- 
lar set of industrial consumers, of a given size and region. 
Thus the results for this variable are no more than first 
approxiriationo . Given that on technological grounds there 
is no reason to expect the absence of an important overlap 
of areas of substitution between electricity and other 
fuels, the statistical limitations imposed on the analysis 
go some \Ta.j in explaining what must be to the theoretical 
economist the rather surprisingly unimportance of price. 

However, the discussions and results have also 
brought out the important part that must be played by the 
associated equipment in any business decision on the clioice 
of fuels. Unlike raw material inputs fuels are demanded 
for services the}'" perform within the production activity. 
Electricity is a part of a group of complimentarj'' goods 
composed of the various types of electricity consuming 
capital. The indices by which the relative price have been 
measured, because they are based on electricity and other 
fuels alone may not be an adequate indicator of price 
relatives of electricity’s whole group of compLiments , 

It is the price elasticity of these compliments which is 
more relevant to -theoretical expectation of raticaial 


behaviour* 



6,4- RESULTS FOR STOCHASTIC TIl-lE SERIES ANALYSIS 
AI3L' FORECASTING OF ELECTRICITY DE}'ANI) 


In this section ve prc’sent the results for 
stochastic time series anal^’sis of electricit;/ demanc’ and 
prediction of the future values of the tire series. Results 
are given in the order of the procedure adopted for identi- 
f 3 ^ing, estimating, validating the model, Tlie future values 
of electricity demand arc’ forecasted hy using the selected 
model . 

5 efore presenting the results ue discuss about 
the data used for the study. 

The data on tirje series of peak power demand and 

energy demand have been collected from various sources. 

\ 

The sources of data are presented in Appendix A, Both 
the time -series data of pealc demand and energy demand are 
discrete in nature, VJe denote the time -series of electri- 
city demand os series A and peak poorer demand as series B. 
Series A is a time series of fortyseven annual observations 
of electricity-. 6o-^umption (in million II'JII units) and 
series E is a series of fortyseven observations of annual 
peal: pou'er demand (in ^IW) for the whole of India. 
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6.4.1 I dentification Of Tr end 

A srapliical plot of series A and T indicates the 
existence of a distinct trend. The trend appears to 1)0 
an exponential function of tiii':c. 


The trend comprnents of series A and T are assuned 
to he of the folio*. 'ing forms : 


Model I (Polynomial trend cor.roonent) 





(5.1) 


Model II (Linear trend component) 


The linear trend for series A and E are assumed 
to he of the following forms. 


E 


P D, 


T 


^o ® 


a^ t 


( 6 . 2 ) 


1=0 ® 


h^ t 


where 



Trend component of series A at time t 
Trend component of series E- at time t 


a,b^, j3 


etc. are the parameters 
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Takinc lociarithjn of koth sides of equations (6.1) and 
(6,2) u'e have 


Loe- 


Log P L 


Log E. 


~ C<p * c?Ci ^ 

= los + fi^t* 

= Log + a. t 


(6.3) 


Log P 


= Log + b^ t 


The trend component is a rionotonicallj^ increasing function 
oi time, without any fluctuation around a nean brend line. 
This leads us to tenta.tively suggest that cyclic conponents 
may be absent, Tlie presence of cyclic conponents have 
to be tested "by correlation and spectral analjrsis. 


The trend conponents are estimated by multiple 
regression analysis. The ordinary least squares teclinique 
have been used forthis purpose. The estimated equations are 
presented below. 


Trend models estimated for serie s A, (Energ3’‘ Demand) 


Ilodel 1A (Polynomial Trend) 


= 2.977 + 0.030 t •*> 0.0003 = 0.996^ DW=1.7^ 

(173.82) (18.7^3) (9.033) 

^fodel 2 A (Linear Trend) 


2.0591 + 0.0452 t 

(154.021) (67.217) 


1^=0.9899 30^=1.963 



Trend modeJ ._ s__e_stini?ted for serlr<^ r (peal: demand) 

Mode l lB _( rol^mominl Trc-nd*) 


1?7 


P. ^ = 2.519 + 0.016 t + O.OO.Cif 
^ (101,78) (6.C50) (9.555) 


H = 0.9897 BVJ = 1.766 


Model 2B (Linear Tr°nd) 

P = 2.339 + 0.033 t 

(84.002) (37.947) 


_2 

H =0.9690 DW=1.59G 


where 

_2 

R = Coefficient of determination of the equation. 
DV/ = Durbin ~ V/atson statistic. 


The values in the parenthesis below the coefficients are their 
corresponding ’t' values. 

_2 

For Models 1A, 2A, 1E, 2B we observe that R values are quite 
hij;^ and the 't' values of the coefficients are significant. 
This leads us to suggest that the trend coriponent is signi- 
ficant for both series A and B. 


We have chosen the linear -trerd for both scries 

_2 

A and B, since the difference betireen E values for poly- 
nomial and linear Hbrendi are . negligible . 

The linear trend compor^nts are then subtracted 

from the original time-se'f’ies to obtain a trend-free 

% 

series , 



6.4.2 Hogults Of Identif ication Of Of 

— , P srsistGncc And Cycli c ity 
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As siuiiini_, that the resi0.t3.nt ahci. trend— free seriES is 
an MaVA process of order (p, d, q) the next step is to 
identify the values of J), d, q. 

The principal tools that have been used for identi- 
fjang (p, d, q) are the autocorielation, function partial 
autocorrelation function and po^-er spectrua of the time 
series , 

The autocorrelation and partial autocorrelation 
functions have been calci0.ated for various degrees of 
differencing. The computational algoritlv' described in 
Box and Jenleins (45 ), has been used for this purpose. 

To identify the degree of differencing (d) to 
produce a stationary ARMA. process use is made of the pro- 
perties of ACB and PACT of stationary series (45). The 
ACF for the non-differenced (d = 0) scries A ar^i B does 
not die out quiciay indicating a non-stationary process. 

Hence these series are differenced once and ACF and PACF 
for different lag values are computed. The first ten 
autocorrelations and partial autocorrelations for the 
undifferenced and first differenced series are presented 
in Table 19. Tests of significance (45) for the ACF 
and PACF (whether they are effectively diffei^nt froia 
zero beyond a certain lag) have been conducted. 
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To identify the existence of periodic components, 
if any tlie snoothed spectra for degree of differencing 
zero, and one, for various lag values have been obtained 
for Series A aiKi E by the nethod described by Jenkins and 
Watts (209). The computel rav and shiooth spectra for 


various lag values are pr:;sented in Table 19.' 


From the characteristics of spectra that has been 
tabulated and recalling the properties of periodic compo- 
nents (209), we are led to suggest that no periodic compo- 
nents are present in both series A and B. 


Based on the information provided by the ACF 
and PACF and spectra of series A and B for degree of 
differencing zero and one, and recalling the characteris- 
tics of ACF and PACF of various AEI-IA processes (45j 209, 
217) the fo?lo\ri.ng orders of the ARII-iA process have been 
identified. Since identification methods are rough pro- 
cedures applied to a set of data to indicate the kind of 
representational model worthy of further investigation ten 
models have been tentatively selected for series A, and 
ten models liave been chosen for series B, They are 
presented below. 


ABIMA models' tentatively aeleetod for Series A and B 


1 . ( 1 , 0 , 0 ) 
5» (2, 0, 2) 
( 2 , 1 , 1 ) 


2 . ( 1 , 0 , 1 ) 
6. (1, 0, 1) 
10 . ( 2 , 1 , 2 )^ 


3. (2, 0, 0) 
7 . ( 2 , 1 , 0 ) 


h. ( 2 , 0 , 1 ) 
8 . ( 1 , 1 , 1 ) 
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Rc s ulL * b ^ pjT initirii. o sting *01,0 n o f rojiG t g 

Tlie initial estinatcs of the parpji;eters oi the 
procesces arc roouirod to he used as starting points for 
laore refined iterative es tination procedures . Tlie initial 
estimates are obtained fr 'ii the estirates of ACF and PACF 
at the idenbifiCi. tion str;,i, l.ie coj'’.putational algorithm 
proviclec' in Eon and Jciiliins (45) have beer, used to obtain 
the initial estim^^tes. A surema-y of the .'nitial estima- 
tes of paraii'.cters for the tentatively selected models are 
presented in Table 20. 

Results of m aximum lll'iolihood es timation of -pa rameters 

Tlie identification process having led to a 
tentative formialation of various models, the efficient 
maximum likelihood estimates of parameters are obtained 
by using the computational algoritlim provided in (45) • 
Marquardt’s (230) non-linear es/lmation routine has been 
used for tliis purpose. The parameters and their variances 
estimated by the laaximura likelihood principle for all the 
tentatively selected models are presented in Table 20* 

Results of validation of models 

The models are validated by moans of diagnostic 
checks applied to the residuals of model. For details of 
the procedure reference is Eia.de to Box and Jenkins ( 45 ) • 
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.^o oi'lL 3 '’ prc.ierib the r-f cults cf diagnostic checks 
for the .:odel which is finall^/ sriccted on the basis of 
these tests, ihc I'looels hiat a.T’o fjnall” selected for 
series h anv.. E arc iatntjr>Gl /J'.hA process, of order (2, 1,0) 
with different va3.uos of nara, t '-. ors . The result, of diag- 
nostic checks are present d in able .11, 


Proiii the auto c jrrelr/'ion IV.ncti''n of residuals 
it is observed that their value' are statr.-tica'ly not 
different fron zero. For a confidence level of 95;.' all the 
smooth spectra values lie beti'cnn thr upper and lower 
limits of and respectively, leading xis to 

suggest tJi’t at the specified level of confidence tlie 
residuals are pure random. Further the chi -square test 
shov? that the calculated value cf is less than the 

critical value at the specified level of significance. 


Eased on the results of the above tests for 
validation of model, we hai’e selected the 'iBII-iA. models 
of order (2, 1, 0) for both series A and E. 

The models that have been finally selected for 
representing the time series of electricity demand are 
presented below. 


For Series A (Energy Demand) 


(1 - 0,3316 B - 0.2?f6l B^) (1 - B) 
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or 

= 1.3316 z^_^— 0»Q855 Z^_2 - 0.2461 H 

•^2 _ 2 , 

= 0.883 X 10 ’’ 

For Scries D (Peak Demand ) 

(1 - 0.3617 B - 0.3631 E^) (1 - E) 
or 

Z^ = t-,3617 Z^_^ + 0.0014 Z^_2 “ 0.3631 Z^_^ + 

= 0.1847 X 10“^ 

The forecast values of peak demand and energy are obtained 
by using the demand models selected. Forecasts of annual 
peak and energy demand for India for the period 1976-2000 
are presented in Table 22. 


TABLE 19 

Autocorrelation Function of Series 
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TABLE 20 (conbinuGd) 
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S,S - standard Error 


TABLE 21 
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Autocorrelatj-On Function and Smooth Spectrum of 



Residuals 

of Models Finally Selected. 


Lac 

— 

-^FOLb _. .A 

1 SERIES B 

8 ATT? ' (prnootn Suectrafl 

' ACF 

iSniootli SDectra 

1 

-0.0186 

0.1065x10**^ 

-0.1080 

0.lC62x10“^ 

2 

-0.0834 

0.1079x10“^ 

-0.1759 

0.060^-x10”^ 

3 

-0.0183 

0.0791x10“^ 

0.1610 

0 .0989x10“^ 

4 

•0.3398 

0.0707x10”^ 

0.2341 

0.2843x10"’^ 

5 

-0.1390 

0.1 280x1 O**^ 

0.1494 

0.3685x10"^ 

6 

-0.0046 

0.1659x10“^ 

-0.1681 

0.3331x10"^ 

7 

—0 . 09^'4 

0.148 1x10"^ 

-0.0435 

0.2715x10~^ 

8 

-0.1051 

0.0781x10“^ 

-0.0092 

0.2^!-64x10‘‘3 

9 

-0 .0209 

0.0516x10”^ 

-0.1027 

0.lC6i-!-x10"^ 

10 

0.0486 

0.11 29x1 0“3 

0.2044 

0.l364x10’^ 
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TABLE 22 

Forecasts of Pea Ji Power and Energy DemajiC’ by 
Time -Series Analysis 

confidence level . 

|Peal c Pove r t s jj E nerg^ T l^^a nd ^Fo^ casts 

YEAR H Ife an rUpper " TLo\;ef "T’ Lean 5' Upper 5 Lower 

5 Value 5 control 5 control 5 Value 0 control 5 control 


5 

5 limit 

0 limit 

A . 

5 limit 

5 limit 

1976 

10160 

11780 

8831 

^550 

27630 

22390 

1977 

10990 

12670 

8710 

26980 

29680 

22960 

1978 

12220 

13710 

8790 

29560 

31340 

23600 

1979 

13650 

14860 

8979 

32660 

34200 

24380 

1980 

15350 

15140 

8642 

35990 

36730 

25410 

1981 

17300 

17580 

8732 

39540 

39620 

26360 

1982 

191^0 

19590 

8872 

42460 

43550 

27540 

1983 

20890 

22230 

9057 

45600 

47860 

28840 

1984 

23330 

25230 

9311 

48980 

52600 

30270 

1985 

^890 

28640 

959^ 

52600 

57680 

31770 

1986 

27100 

32510 

9931 

56490 

63240 

338IO 

1987 

29650 

36990 

10300 

50670 

69180 

37 180 

1988 

32280 

41780 

10720 

65310 

75680 

39540 

1989 

35^00 

47420 

11190 

67010 

82990 

39170 

1990 

395^0 

53700 

11720 

75500 

90360 

41400 

1991 

42170 

60670 

12250 

80910 

99080 

43750 

1992 

46230 

68710 

12880 

86900 

108900 

46340 

1993 

50510 

77800 

13550 

93510 

117200 

48980 

199^ 

54990 

87700 

14290 

100500 

127900 

51 880 

1995 

5998c 

99080 

15070 

107900 

139300 

54950 

1996 

65610 

101700 

115920 

115900 

151200 

58^0 

1997 

71610 125900 

16830 

1^500 

165200 

61 800 

1998 

71860 129100 

17820 

133700 

179900 

65460 

1999 

85510 159200 

18800 

143500 

195900 

104700 

2000 

93330 

179100 

20000 

154200 

212800 

116900 
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6.5 Discussions OF RESULTS OBTAIIffiD 
FOR FOREaiSTS m USING 
THE -SERIES ANALISIS 

; ARIiA models oI order (2, 1, 0) \d.th differinr 
pararr.oter values have been used for forecasting peak power 
and encr£3^ demand. The models identified for the processes 
generating these eenand series turn out to be non— stationary’’, 
ilon-stationarities of the demand serj.es are confirmed by*' the 
existence of a ti^end component i/hich is an increasing func- 
tion of time. V/e observe large differences in the magni’tudes 
of forecasts obtained by the econometric methodology and tjbne- 
series analysis methodology. It is very difficult to provide 
specific reasons for this \dde differences, A possible 
reason might be tlio.t the parameters of the ARIK/l models have 
undergone structural changes over the tine span of observa- 
tions , The tjLme-series analysis methodology suggested in 
this dj.ssertation 1 not pro’vide accurate forecasts in t!« 
e-ventuality of time variant parameters. Further from the 
length of the record available it is difficult to draw any 
concO-Usions regarding the nonstationarity of parameters. 
Adequate leng’tli of records are essential to es’timate a re- 
liable model of "tlie process generating the tine-series.. 

In our S’tudy the length of the time •‘Series is not adequate 
for a realistic and accurate representation of the process. 

The possibiLlity of errors of measurement of data is not 


ruled out 
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In the nature of things it is difficult to realise 
the assumptions, based on which projections have been made. 
However we have souglit on the basis of data available some 
independent estimates of likely levels 6f deriand. Generally 
it has been found tliat estimates resulting from those based 
on macro economic variables provide reliable forecasts . It 
is expected that a.s more and more data becomes available 
ccmparison of actual values mth fore cats will enable a 
meaningful revision of forecasts. We are reasonably satis- 
fied that .unless there is a development of a kind of an 

uneapected nature forecasts may be treated as adequate and 

♦ 

realistic. To have quantitative Informations on the risks 
involved ifith plans based on these forecasts we have added 
the probabilistic dimension to forecasting. It needs to be 
pointed out that it is not possible to have the same measures 
of confidence on the long term forecasts beyond 1990. The 
experience in the next few years will provide a clear indi- 
cation of the direction in which forecasts have to be revised. 

The demand for electricity depends on many factors 
such as future pattern of economic growth, demai^ for the 
output of electricity using coinnodities . In a country like 
India which has a central planning system, lin diPSiilaifej'Pp 
the plans factors such as regional development, inter- 
sectional relationships, imports, exports ard. social 
fare are given due ■vjeightage and final targets are deter- 
mined based on these considerations, lowever, It slwjuld 
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recognised that while plans represent desired levels of 
econoiiiic development in accordance vdth desired social 
objectives their execution depends on a variety of factors. 
Hence forecasts based on planned and anticinated growth 
rates may not at all be rerO.ised due to improper implemen- 
tations , 

Forecasting of elactricity demand, in its truer 
meaning and form involves substantially more than the 
application of ratios, formulae and historical relationships, 
it is very much a venture into the unknown, a journey into 
the land of 'if’, 'probably' and 'perhaps'. Any study of 
the future requires exercise of judgement. However, judge- 
ment is a much abused term. At its worst it is an invitation 
to whims and conjecture rather tlian an adjunct to sound 
rigorous thinldLng. Much of what may pass for judgement is 
in truth relies heavily on a sour^ methodological frame \rark* 
Forecasting being both a science ar^ art, scientific tools 
and techniques are as essential as judgement aixi insight. 

Eml ngufi 

After studying the material presented in previous 
sections the reader is clernLy jiistified in saying "this 
was all somewhat informative, but T^t are tl» 3?eal demand 
projections?" The ansi^er to tee implied criticism in the 
querry must be that we cannot offer definite d®BMd 
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functions , that electricity demand functions are evolutionary 
variables v/hich must be re-exeimined and changed when impor- 
tant determining factors arc changed, and that the deter- 
mining factors are made up of an ill defined inberixLay of 
physical, economical, social and political forces. Elec- 
tricity demands are time dependent functions that reflect 
the status of our technology, our economic v;ell -being, our 
social system, our leisure time habits and existing political 
realities , They vn.ll change vriLth tine . 



CHAPTEE VII 


PLAM'HIJG FOR GAPACITI EXPAI'EIO'.I OF ELECTKCC 

POVJER SYSTEM - A LITERATORE SURVEY 

The literature on planning investments in electric 
po%'?er systcns to satisfy future demancis covers a fairly 
"broad spectrurn and is both qualitative and quantitative. 

In this chapter ■vre present a brief review of the literatxire 
in the area of planning for capacity expansion witli special 
emphasis on electric power systems. 

A principal difficulty in determining niimial cost 
optimal strategy for capacity expansion is that the invest- 
ment costs for capacity increment typically are subject to 
significant economies of scale, lia thematic ally tliis means 
that the investment cost functions for capacity expansion 
are non-convex. Thus the determination of a minii'isl. cost 
expansion strategy is Ittely to entail the inlntalsatioii of 
non-oonvex function and we are faced with the problaii of 
distinguishing the true glotal minimum from pos stole local 

minima. 

Even the static location models, in uhidi demands 
are fixed for a given point of time and it is desired to 
find plant sizes and locations to mlnloise total invest 
ment, operating, and distributing costs for meeting this 
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demand suffer from the problem. Most of ths models of 
traditiana3. location theory are of tMs form. A survey 
of the literature of such models is given by Bos (237) 

For the case 'vriLth demand of uniform intensity over an 
infinite plane surface and a cost function composed of 
a fixed charge plus cost proportional to the size of ex- 
pansion, he shows that average cost per unit of dci-and 
may be transformed to a convex function \diich can be solved 
directly. If mar3:et boundaries are finite and demand are 
concentrated however this sinplifi cation is ruled out. 

Baumol and V/olfe (238) developed a simple mathematical 
programming model for warehouse capacity problem vdth a 
finite number of demand points and locations which as they 
pointed out would find only local optima. Kuehn and 
Hamburger (239) demonstrated that such local optima might 
be poor approximation to tlie global optimum and developed 
a simple heuristic method to find reasonably good solutions 
for location problem. Other heuristic methods have been 
developed by Cooper (2^0,241 ), Feldman, Lehrer and Hay (2^2), 


The static location problem has been formulated 
as an integer programming problem with zero-one integer 
variables corresponding to decision to not construct 
or construct expansions at varioxis locatiois. For sswU 
problems solutions have been found by the isethc^. of 

# 
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complete enumeration by Vietoriss and Manne ( ^3). 
Approximate solution to larger problems have been obtained 
by the "one-point move" algorithm due to Manne (244). 
Efiiomson anc. Ray (24-5) suggested a pror.jising approach 
for solving this specialised problem by uso of braiich and 
bouna integer prograniming methods. Tlic case of uncertain! ty 
in deniano. for static location problem has been studied by 
Gregory'- ( 246) . 

The problem of capacily expansion over time for 
a single location has been studied by Qienery ( 247), 

Me Do\rell (^8 ), Manne ( 249), Coleman and Yorh ( 2 ^), 

Manne ’s work also considers a special case of probabilis- 
tic demand growth. Veinott and Wagner (251) have observed 
the mathematical equivalence of capacity expansion problem 
vri-th otlier \- 7 ell knoTO problems like determination of econo- 
mic lot size inventory decision and equipment re-.'lacement 
policies . The single location capacity models is closely 
related to the inventory model of Vflsitin (252) and Hadley 
& VBiitin (253). A dynamic programming apnroach for single 
location problem -with finite time horizon and arbitrary 
increasing demands is given in Manne and Yeinott (254)* 

The general problem of planning capacity expan- 
sion for several locations with demand growing over time 
has been addressed by Ghosh (255) in case of Cement 
industry in India. He did not incltrf^e the ccusplicacies 
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of econonies -of -scale in investment cost, and solved the 
problen by linear prograiimiins, transportation model. 
Application of linear prograi-jiiing to capacit'”' pla,nning in 
various industries have been sunmarised by Ward (256), 
ihe additional element of econor.ies— of— scale in investment 
costs for expansion is included in the models developed by 
Manne (257). ^''or the case of tuo producing locatiois, 

liannc uses a simple two-phase cycle -model in which the 
time interval betvjeen expension is constant and identical 
for both locations. For the case of several loc-ations he 
developed a heuristic method ^diich assumes a constant cycle 
time interval for expansion at individual locations. 
However, the model permits the length of the cycle to vary 
among locations, provided that each such interval is an 
integral divisor of some longer period called major cycle. 
The problem is forimflated as an integer programming model. 
The solution method utilizes heuristics and does not gua- 
rantee optimal solution, ffenne (257) applied tliis method 
to three industries in India. The industries considered 
were Cement ■'.ri.tli ten locations, Caustic Soda ano. i-ertilizers 
with fifteen locations each. 

Another integer programming formulation coupled 
with the use of a branch and bound method has been used 
in investment planning study for tire Brazilian Steel 
industry by Kendrick (258). Kendrick- s model is n finite 
horizon model and is quite comprehensive. 
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Tl c dynarjiic mid.tilocation capacity planninr prdblcr.. 
is also closely related to the multiproduct nu.1 ci -facility 
inventor^' uod^cl of Zangwill (259 )y with capacity irscrcipcnts 
at diffeicnt locations corresponding to orders for d.ifferent 
products or the so.jC product at different facilities. 

Sreedliaran and Woin (260) have proposed r conti- 
nuous tine model v/ith several types of plants. Tlioir 
emphasis is on the optimal timing of several candidate 
sequences of plant installations. Erlerhottcr (261) con- 
siders a model in which the type of plants are dcd.c mined, 
by their locations . Ec seeks the optimal amount and tiining 
of capacity expansions over a finite horizon. In his paper 
the system operating cost is the cost of transporting the 
product of plants to the points of demand, lie also dis- 
cusses a stationary planning model irith an infinite horizon, 

• In planning the expansion of electric power S7,rstems 
there are two popular typos of models : simulation and 
programming. Galloway et, al. (262 ) have discussed a 
model to assess the stochastic variations in the available 
capacity. Giguet is one of tlie earliest investigators 
to deal \jlth the problem of capacity expansion of elec- 
tric power systems. The main feature of tlic model is 
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to study the relative profitability of a p^rticul.-.r plant 
in i-<ilation to a reference thermal plant. All plants 
;;hich arc con-jidcred to be potential plants arc raiilccd on 
the basis of their relative profitability. The liLants 
which have the highest profitability arc selected until 
the total demand is satisfied. Tliis method is satisfactorjr 
as long as various plants do not interact among themselves 
and there are no econcxnics -of -scale involved, rnc cut-off 
method of Giguct has been an acceptable mctliodology till 

195 ^. 

Some investigation into tlie area of plaining the 
capacity expansions of pokier systems have dealt with the 
classical project by project engineering, economic and 
political trade-off between hydel, csteam and nuclear power 
generation, Jacoby (263) gives a good description of the 
methods presently used in practice. Host of the methods 
have either used linear programming or simulation. 

The French Federal Power Comniission has clone 
considerable work in using linear prograijning in numerous 
studies of investment planning in electric poorer industry. 
Masse and Bessier (26^), Lfesse and Gibrat (26^) are some 
of the important contributors. 

The model developed by Dantzig (266) is in the 
same lines as Kasse and Gibrat* s model except tJiat it 
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includes luc tijic dimension. 


In tills inodel elenenc of 


power transmission are not considered. 


Here recently the Llectricite' dc, Trance- (TDI ) 
published the invcstiient >85 riocel (267) conrictinr of 
159 variables and 53 constraints. The co2T\-ex quadratic 
cost function and the linear constraints assumed in the 
model uses the Wolfe's reduced gradient optimisation 
routine. This routine makes use of the fact that the cons- 
traints are linear? and hence employe partially the sinplex 
procedure of linear programming. 


Lack (268) formulated and solved a linear djmamic 
model of the diversity problem in which the problem of 
determining the trade-off bet;;een adding more gener''ting 
stations at one of the two inter-connected load centres or 
adding more transmission capacity was considered, Chen 
et, al (269), and Sautter (270) formulated and solved the^saiiie 
d.iversity problem by using linear prograiiiming Models. 


Sequential probabilistic linear progranriing was 
employed by Ilanne to calculate the optimal electricity 
plant mix decision during the decade 1980* s given the un- 
certainities of the date of availability'’ of breeder 
reactor. The model allowed for the possibility that 
future uranium resources sca.rcities might lead to an 
increase in electricity prices and hence a reductlcm in 
projected demand. 
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liarasirnlian’s (271) lineal*- prograniminc model 
considered the transmission of pov/er bet;/cen generating 
anci load centres along with the capacity planning prcblen 
of determination of size and location. Capacities and 
energy capacities during ench season of nev/ plants were 
t alien as the variables of the model. Die objective func- 
tion \T3.s sum of cc3t of total energy for each plant wiiich 
included fixed '.nd running costs. Constraints of the model 
were upper bounds on capacities and energy balance of the 
entire system. In tlij.s model transmission of power mbs 
talien to incur on^.y fixed cost of transmission equipments 
and transmission losses were ignored. 

In a dissertation by Gosai (272) the L.P. model 
has been used to determine location, type and capacities 
of power plants to be added to an existing system. Ihe 
concept of product rhjc and integrated grid system was 
incorporated in this model. The objective function for 
his model minimisco. the annual amortised cost, annual 
generation cost and annual transmission cost. The ccxis- 
tralnts tool' into account the power losses in the system, 

Tikaria (273 ) has developed a nonlinear progra— 
mning model of planning for capacily expansion in a 
grid. The objective was to minijnise the sum oi trans- 
mission loss cost, annual amortised of plants and 

running cost. The model formulated has a cpadisatic 
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olDjecti\T 2 function linear constraints, vliich has been 

solved by using sequential unconstrained minimisation 
algorithm. 

In England Dale (274-) suggested t’nc use of dynamic 
prograr.Tuing for the problem of do termination of size and 
time of installation of ■oo\rer plants, neglecting the loca- 
tion aspect of the problem. 

Integer line'll’ programming formulation has been 
used by Qliada and Ilauai (2'/5), to select sites, tinje and 
capacities of nci^ plants •‘..i order to meet the specified 
demand \d.th minimum cost. A step^/isc cost capacity’- cliarac- 
tcristic has been considered in this fomid-ation uMch is 
more realistic than a linear characteristic, irliich has been 
assumed in all linear programming models . But pover losses 
in the systcfxi and upper and lo\;e:’ bounds on capacities of 
plants have not been considered. 

Various boohs and articles have been piiblished on 
the system simulation approach to the vholc problem of 
planning size and location and time pliasing of future 
series of pov/er plants . Host of the publications present 
fairly complex models but use no optimisation. Sensitivity 
analysis is used in most cases to obtsiin the various good 
feasible solutions, Marglin (276), Helson (277^ 

Bary (278) present very ccciprchou&ive descri-ption of such 
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aspects or the problem as pricinc, cost to serve, inter- 
connection l'ct\;ccn systems etc., but leil to give r 
cor.putr.tiona3.ly efficient approach, that uill include all 
t}ic cconoi.dc aspects of the problem.. 

Cazalett (279) presented an extension of the 
Everett's generEilised lagrangc multiplier approach to 
iincons trained minimisation problem. Ilovcvcr only the 
strategic plaiit capacity sub-problem vas analysed. Hic 
tacticaD- variables i.e., those related to operation by 
variaas nlants in the sj'-stem were appro'ar,ntcd by using 
a h 3 rpotlie tical operating policj'- that loads plants in order 
of efficiency. A fairly complex simulation oriented 
to^■/ards Nuclear plant feasibility analysis was published 
by S.R.I. (280) in conounction with Lexican Federal Power 
CoiiUTission 'using Cazallet's concepts'. 

Jacoby (263) presented a ver\^ comprehensive 
simulation model which analytically included many of the 
ideas of Baldwin (281,282 )and Bary (278). Although 
the approach included operation, reliability, environ- 
mental, social value and evaluation submodels it did 
not include any optimisation procedure to determine the 
best schedule of plant installation. 



The sur\'iey of the reported literature reveals tl^at 
nathcmatical programming models have been extensively used 
for solving the problem of capacit3'' exnansion of electric 
poucr systc- . ilovfcver the models proposed, have trcl-acd 
t]:e proble;., of capacity’’ expansion and operational plaxining 
independently and scperatcl'"-, \dth canacities as coiitinuous 
variables, v.'e r=>roly find any evidence of a methodology 
uliich obtains an integrated solution of tlie overall capa- 
city expansion problem incorporating the capacity ox'‘)ansion 
aspects and operational planning aspects simultaneously. 
Further most models proposed, treat the capacities as conti- 
nuous variables, where is in a real life situation capacity 
increments take place onl3^ in discrete steps, due to stan- 
datd in force, availability of specific sizes of plants 
in the marl-et, and manufacturing considerations of power 
plant equipment. A capacity expansion program wliich does 
not consider the feasibilities of ox'oa.ision r.t a pa,rticular 
location of a given type and size of plant in a given period, 
is uhlikel3’' to provide a pragmatic solution to the capacit3’' 
expansion problem. If it were possible to draw up a com- 
plete list of all feasible power plant developnents in a 
region bhe problem would reduce to choosing the best combi- 
nation and sequence of power plants to meet the projected 
power demand over the planning horizon. But such a list 
of feasible development multiplies as soon as several 
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locations, sizes and construction periods and type of 
projects, (such as hydel, thermal and nuclear plants) are 
talien into consideration. Because of the problem of dimen- 
sionality most studies up to nou have used a screening type 
of studies to obtain a few possible power plant sequence, 
but have failed to obtain the best solution. 

In this dissertation general models that include 
the main aspects of electrical power system capacity expan- 
sion problem are proposed. The approach is neither purely 
simulation nor analytical oriented. Model simplification 
ha.s been made in the interest of obtaining numerical results 
bub such realism as transmission lines, operating policies 
are included. Gradual improving expansion policies are 
obtained and each time their feasibility with respect to 
satisfying both capacity and energy are tested. 

The complexities of power system planning is due 
in part to the technical interaction between generating 
plants and economic impact of tneir interactions. The eco 
nomic importance of fuel costs require the analyst to study 
the-cbordinated operation of the individual plants as well 
as capacity expansion of the system. The large number of 
possitle combinations of individual plants with a realistic 
planning period has been the major obstacle in us 
optimisation techniq'ues for planning the capacity expa^is* 
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of electric poorer system. Fodel deconposition and iterative 
concepts prove very useful m partly resolving this combi- 
natorial problem. 

The purpose of system modelling is to construct 
a. mathematical model such that response of both the real 
system and matnematical model to the same input are almost 
id.entical • But a mathematical model cannot includ.e a!}.l t}!io 
detailed aspects of a large scale electricrl power system in a 
developing region. Thus the aim m.ust be to represent as 
realistically as possible on3.y those aspects of the problem 
I’-hich are important. The theory and applications of the 
mathematical model of a coEiplex system arc still in an early 
stage of development, A complex realistic model with many 
locaD.ly good solutions ideally requires an efficient optimi- 
sation routine that vdll find the globally optimum solution. 
However, this is not alviays feasible and the problem is again 
one of determining the best trade-off between efficiency and 
cost of optimisation techniques. 

With the aforementioned considerations an attempt 
is made in this dissertation to suggest a methodology for 
determining the minimum cost expansion of capacities of an 
electric power system in a region incorporating the capa- 
city expansion aspects eb" well as operational pi&Sjiiiiig 
aspects in an integrated manner. The suggested approa^ 
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systematically searches for an optimal solution. among 
feasible alternatives. The iterative use of capacity 
expansion planning models together with that of the opera- 
tional planning model is is \^at distinguishes our metho- 
dology from other mathematical models. In Chapter VIII 
v;e present the mathematical formulation of the capacity 
expansion problem and the operational planning problem and 
their solution methodologies. 



CHAPTER VIII 


PROEIJSK FOEl-IOLATION AND 
SOLUTION lETIICDOLOGY 

This cliapter is devoted to the problem of 
planning for capacity expansion of an electric power systcu. 
Section 8.1 gives a statement of the overall capacity ex- 
pansion problem and the need for a decomposition method 
to solve this problem. The overall capacity expansion 
problem is treated as a combination of tvzo sub-problens, 
viz. capacity expansion sub-problem and operational planning 
sub-problem. Section 8.2 presents the mathematical model 
of the capacity expansion sub-problem. In section 8,3 
\iQ provide the solution methodology for the capacity expan- 
sion sub-problem. The formulation of the operational 
planning sub-problem and its solution methodology is pre- 
sented in section 8,4. 

8,1 PROELEM STATEMENT 

The purpose of this study is to analyse using an 
integrated approach the cemplex problem of electricity 
supjxLyj loads, and econondc constraints in planning 1^© 
best location, size (capacity), lyp®* and time of insta- 
llation of a series of power ^ants in a region to laeet 
the fnture growing demands for electricity. 
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In the analysis presented here, the overall capacity 
planninj^ prohleci is decomposed into tvo parts . They are 
referred to as tlie capacity expansion sub-problem and the 
operational planning sub-problem. 

^ • Ce-Pd.eity; Expansion Su b— probLem - Tliis problem involves 
the deteriaination of a minimui'! cost (discounted) secuence 
of future plant instaJ-lations such that the total system 
which might consist of a nix of hydel, theimEa and 
nuclear plants has at all times a total capacity suffi- 
ciently large to meet the peak load demand plus reserve 
requirements. This problem is referred to as the sjy^^- 
ter<ic problem. 

2, Operational Planning Sub-problem - Tlais problem invol- 
ves the determination of the discounted Eiinimum post of 
operation' of stich. a. sequence of powfer plant -ingtallafelons. 
The sub^probleai is inferred to as the tactical problem. 

From the operation planning sub -probLem the pro- 
duction and other system characteristics arc calculated. 

The parameter values obtained from the solution of the 
operational planning sub-probLem are then compared \ri.th 
parameter values assumed in capacity expansion sub-problem. 

If there is a difference in assunqption of -aVbrage opera- 
tion in (1) and the calculated operation results in 
(2) the assumed parameter values, used in t}» capacity 
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c:q)ansion sub-problems are altered. This feed back is 
repeated until the expansion sequence in (1) lias no further 
inprovcnient. The same approach has to be carried out for 
different demand schedules. Ey the methodology of decom- 
posing the prolfLem into two parts, the primary purpose is 
to find a quasi -optimal solution to the overall large scale 
problem. 

The strategic problem corresponds to the capacity- 
eripansion cub-nroblem and is the portion of the procedure 
that emphasises the capacity installations (MW) aspects, 
i.e., determination of quantities of different types of 
equipments to be constructed. 

■ Tlie tactical or operational planning sub-problcn 
emphasises the operational aspect of the plants chosen bj?" 
the strategic calculations. In particular the operational 
planning problem determines the minimum variable cost 
operating policy to satisfy the projected energy demand, 
(thB9' providing the feasibility of the current expansion 
policy) . 

The feed-back between the two problems gradually 
improves the f\iel cost estimates so that a true minimum of 
fixed plus variable costs are ultimately obtained for the 
overall capacity expansion problem. If the objective is 
the minimisation of fixed costs orily, rather than fixed 
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plus variable costs, only the capacity eiipansion problem 
need be solved, .since investment cost for each alteri'^e.tivc 

t 

plant is hnovni. Howeverj the variable cost of each plant 
are not knoi'm as they depend on the alternatives chosen b 3 ’’ 
the capacity expansion sub-proTaLem . The variable costs 
■;/ould be knovm if these plants operated in isolation fro:;i 
each other, i.e,, vdthout any systems co-ordination. 


The above methodology is based on the observation 
that tentatively fixing the values of the plant energy 
variables renders the problem simpler, wliich is solvable 
in this case by an integer programming algorithm. These 
plant energy variables (which are tentatively fiasd) are 
referred to as the complicating variables. Benders (283) 
was the first to formalise a decomposition procedure for 
solving such problems , By fixing the complicating variables 
given problem \ra.s reduced to a linear programming problem 
(parametrised by the complicating variables). 


The general formulation of Bender’s semi-linear 


problem is given by 
0 

IJaxiiolse | x + f (Y) • 
X 0 L- — J 

■y X 


( 8 . 1 ) 


such that 


Ax + i‘(y)-b>0 


X, y arc vectors 
A and F arc matrices 
f is a scalar function, 

Tlic code proposed by Bender ( 283 for findin,- tlic optimal 
value of the complicating variables, involves two steps. 
The steps are: 


(1) llanipiilate the above formulation of the problem by 
projection of space of T variables (I-iefer (28^)) 
to yield 


Maximise... 

y 



such that 

y £ I pi V 

where v (y) = Supremum | x + f (y) i 

U- — ' 

X 


( 8 . 2 ) 


such that 

A X + F (y) - b ^ 0 


(8.3) 


and 


V 



: A X + F (y) 


b 


0 for some 
X ^ X 


So the projected problem equivalent to (8.1) is that 


an optimal solution y of (8*2) readily yields an 

ijt 

optimal solution X,y of (8.1). 


190 


(2) Since v, and V are only knoim implicitly ‘solve (8,2) 

• * 

by a cutting plane method tlib.t builds up a tangential 

^ • • 

approximation to v and V, Linear prograi'.aning 
duality theory is employed to derive the farily or 
cuts characterising their representation. 

The electrical pov/er pla 'nine, problem in tliis 
dissertation can be formulated in an integrated form as 

flax P- (x) - f (y) t 

— 1 

s.t P (x, y) ^ 0 

where 

X = vector of 0, or 1 variables that represents 
decision to build or not to build plants and 
transmission lines, 

y = vector of annual plant energy productions. 

The vector of constraints P (x, y) 0 includes satisfac- 
tion of capacity (IM) constraints, i.e,, P^ (x) ^ 0 and 
satisfying the energy constraints i,e., P 2 (x, y) 0. 

By using the non-linear convex duality theory 
Geoff rion (285 ) has recently extended Bender’s approacli 
to a broader class of problem in viriich tlB jjarametrised 
sub— problem need not be linear. A master program is 
developed that would direct the sequence to final dhoices 
of the complicating valuables , Tho computational algori- 
thm consists of iteratively solving {8*3) a relaxed 


m 
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master problem that ignores some of the original inequality 
constraints. The solution to the rela 3 ^d master progrcu^s is 
such that (8,3) is feasible. 

The effec civeness of the above algorithm is limited 
ho^rover to a special class of problems such as variable 
factor prograriTiing problem (286 ) where the solution of the 
relaxed master problem do not too often lie outside V. 

Here a direct non-linear progranming technique 
can be used for the operational planning sub-problem clue 
to the availability of solution algorithm for capacity 
expansion sub-problem and operational planning sub-problem. 
One of the advantage of this method is that by analysing 
the operation planning problem results at every iteration, 
one obtains useful insight and information. This can be 
used to modify the capacity expansion sub-problem, at the 
next iteration. 

8.2 PORMJLATION OF THE CAPACITY 
EXPAIISION SUB-PROELEII 

Before formulating the model we discuss in brief 
about the necessity of a coordinated system design in tne 
following paragraphs , 

Various types of power generating plants avai- 
lable today, have comparative advantages and disadvantages 
depending on several factora, Ti^refore a coordinated 
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system design and operation of the different pover facilities 
is. necessaiy to minimise the total cost of meeting electri- 
city power and energy demand. 

A planner is not only concerned \;ith the coordi- 
nation between various type of plants that are devised but 
also \7ith coo ruination on a multiregional basis , betvrcen 
different interconnected systems. The problem of coordina- 
ting an existing system in a region is different from that 
of planning the coordination of fiiture expanding systems. 

To properly solve the later problem it is necessary to find 
out the best future system which consists of the existing 
system plus sepuence of plant installations throughout the 
planning horizon T, and then to carry out a system operation 
study which will determine tl:^ best annual operating policy 
of the entire plant mix. The operating rules will be 
changing as the systems gradually expands over time T, This 
problem is called the operational planning problem, to bo 
dealt in section 8,4 However, the minimum cost sequence 
of plant installations vTill be knoim only if annual 
cost of various feasible sequences are known, aiMi this re— 
<iuires knoifledge of system load factor, and the best 
operating procedure. The protG-on is apparently circular* 
This is the reason why this study suggests the adoption 
of a tvro step methodology, presented in figure 8 .I* 


193 



STOP 


mo. 8.1 : Diagranmatic representation of the 

msthodology# 






19 ^ 1 - 

The capacity expansion sub-prolilcins uses the 
criteria of minimisation over T years j the discounted, fixed 
cost plus average varia'M.e costs* The variable costs arc 
mainly fuel cost and depend on average plant factors of 
the plants in the system. The decision variables in tliis 
problem include the capacity or size (IM) of different t;-oes 
of plants, to be built in a number of available sites over 
a planning horizon of T years. 

One of the distinguishing characteristics of the 
investment, operating and maintenance costs in the electri- 
city industry are that they are non-linear in nature. In 
addition, the problem of choosing minimum cost expansion is 
an integer valued problem, since only standard size equipments 
a.re available, for capacity expansions. Also, the problerj is 
a fairly large enumcrative effort as soon as a feu alterna- 
tive capacity increments for each typo of plant arc introduced. 

The general integer programming model differs from 
the linear programming formulation in that the later uses 
continuous variables . To use the implicit enumeration s cau- 
tion technique an even more restrictive constraint is re- 
quired, that is the decision variables must have only one 
of the two discrete values, zero or one. 

The general 0 — 1 integer prograwiing protleia 


is of the form 



3 = 1, 2 
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... - are the alterne.tive load centres consi- 

dered tliroughout the region, 

t = 1 , 2 . . . - arc the time period into which the total 

planning horizon T has been devided. 


8.2.2 ^fodel Foironlation 


The formiilation of the capacity expansion sub- 
problem is presented below. 

ACi A p- 

lanlmise ^ ^ 31-32, tlf 


such that 


^ijt = 0 1 

yio't ^ 0 ota 
Zi-t = 0 on 


t = 1, ... 

for i = 1, ... N 2 

j — 1 j ... 


(G.6) 


p ' ^ p ^ P P 

A, - "5” L.. + R.i • X,,. + ^ H. . V. .. A D 

30 X 3t 7 .^ X_ 13 I3t ^ Z- 

t=1 t=l i 3 i 


i3t • ^idt^ ^jt 
(8.7) 


( 8 . 8 ) 

for all t = 1, ... 1I.J 

31 = 1 , ... N 3 

32 = 1 ^ •##' 



^i,0‘lT32 " ^ transmission capacity of the trans- 

mission line alternative between load centres 
31 and 32. 


^0,31-32 


• t 

Plf31-32 - 








2 ... 


The MV/ transmission capability of transmission 
lines bet\/een centres 3I and 32 at the benirjning 
of planning horizon (T = 0). 

(1 + A) where A = percent heat and other losses 
th 

in the i transmission line between j 1 and 32 
during period t« 

Total base and peaking capacity respectively, 

'i'V 

of the 3 ^ region at start of planning period 
(t = 0). 

Total base and peaking capacity requirement 
respectively in the 3^ region at the start 
of t^^^ period (includes a certain reserve 
requirement') . 

Base load and peaking load plants respectively, 
retired in the 3^^ region at the start of 
the t^^ period. 

- E^fcents the alternat|ve plant types 
(size and capacity) available a|j|d alternative 
transmission line (size and cap^|&l 7 ^ available. 
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Definitions of Terminology 



0-1 decision variables that chooses the i^^^ 
type (representing size and type) during the t^^^ 
pe riod at the j ^ location , 



0-1 decision variables that chooses the exr^ort 
of i"*^^ type (representing type and amount) from 
the 3 ^^ region during period t. 


Z j. -0-1 decision variables that chooses the i 

i,0l-32, t 

transmission line alternative between load centres 
3 I and 32 during period t. 



- The swiie of discounted fixed plus estimated variable 
plant investment, operating and maintenance costs 
incurred in choosing the i^^ alter3;ra.tive plant in 
the 3 ^^ location during period t. 


V - The sum of discounted fixed plus variable 

U,3l-j2,t 

transmission facility, investment, operating anti 

. , . th 

maintenance costs incurred in choosing ine l 
transmission line alternative between load centres 
3 I and 32 during period t. 

- Plant capacity instaUation increment allomed Ky 
the X. . alteniatlve independent of time period t. 

- The import or export capacity of type and sine 1 
at 3 "*^^ location, during period t. 
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n 

Minimise ^ C. x. 

3 3 


(8.4) 


sulDject to constraints 
n 


\ a. . . X. ^ bj 
^ 13 3 ^ i 

3=1 


Xj- = 0 or 1 


1 

3 


1 » 2 , 

1 , 2 , 


.. m 

n 


In this section a zero - one integer prograEing model for 
. the capacity expansion sub-problem is presented for deter- 
mining the best strategy for' capacity expansion. 


8,2.1 Model Descrj-ption 

The capacity expansion problem will be formulated 
as an integer programming model. The decision variable is 
to choose an installation policy ( » ^ijt^ export- 

import policy expansion of electrical system 

in a region. By a policy ve imply a complete list of plants 
and trar^mission facilities to be installed and quantities 
of electricity to be exported or imported between difierent 
sub-regions over the planning period of the analysis. The 
planning period T is divided into n number of 5 yss-i 
plans during which a decision to build or not to build 
(^jt = 0 or 1 and = 0, or 1 ) and to import or 

not to import ( = 0 or 1 ) is »ade. 
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and ^ ^ 

l-jSdi-3-2, 


^i}31-j2,ti'*‘ ®o,j1-j2 
J 


(8.9) 


iiiQUation (8,5) is the objective fiinction v;h-ich 
minimises the sum of discounted fixed costs and total esti- 
mated variable costs and operational and maintenance costs 
lor all plants and transmis sion lines . Equation (8.6) is the 
zero-one restriction required for the use of enumeration as 
a method of solution. Equatiors (8.7) and (8,8) are tlie 
restrictions necessary to keep the total available peal:ing 
and base load capacity at the J region in every period t 

at least as large as the required peaking and base load 
P 

capacity which is obtained by projecting the 

P b 

present load requirements in the region. include 

certain reserve capacity which is necessary for meeting any 
expected plant failures. Because of transmission losses, 
^ijt actually a reduced import in equations?. (8,7) and 
(8,8), Eq-uation (8,9) restricts the capacity of the trans- 
mission lines between load centres. For two particular 
load centres j1 and j2 the total transmission capacity 
during period t must be greater than or equal to the 
maxiinuin of the capacities imported or exported bet%reen jjl 
and j2, increased by the amount of losses expected for 
the transmission lines. 
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To calculate the K8s and C*s of the otic ctivc 
function it is first necessary to determine the typo of 
plants and oransmission lines that can be built in a region. 
Kno\;ledge o^' various typo of power plants and transmission 
lines available today and in future are nccessar]’’ to i.alic 
estimation of those co&ts. In the present capacity exoansion 
model total variable costs for various typo of plants alter- 
natives are estimated by .-assuming approximate average 
plant factors j since the details of the co-ordinated 33^5 tern 
operation is not known until the sequence of installations 
are chosen. 

In general, for most regions there will be a largo 
number of alternative power plants (t3T)0S and size) and sites 
(i and j of the model) . The planning agencies usually carry 
out site screening studies to limit the number of choices 
on site and sizes and type as much as possible. Of course 
it may so happen that constraints also determine the parti- 
cular site and tyj>e of plant to be built in the next fe\r 
years. In order to solve the capacity expansion problcn 
formulated in this section it is necessaiy to carry out 
such screening studies based on several factors, 

8.3 SOLUTION l^THCDQLOGI FOR THE CAPACITI 
EXPANSION SUB-PEOBLEM 

The implicit enumeration zero-one inte^r variahi© 
algorithm proposed by Ealas { ) has been us<^ for solving 

the integer prograimaing model of capacity expansion 
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sub-prolDlcra. The algorithm uses heuristic tests suggested 
"by Holconhe (287), because of its high efficiencj’’ and lov; 
cost relative to other integer optimisation routines lihe 
the branch and bounds method by Land and Doig (288 ) , and 
Statistical sampling method by Bciter and Sherman (289). 

8. If FORl^LATIO:? OF TIIE OPEBATIOIIAL 
PLAIII'IIIJG SUE-PROELEII 

In section 8.2 the formxiLation of the capacity 
expansion sub-problem to meet future load schedules at 
minimum discounted costs v/ere presented. The capacity ex- 
pansion model determines the best sequence of future capacity 
installations in the region. In this section present tiic 
long-term operational planning model. 

The fuel costs are different due to the types of 
plant, economies -of -scale, technological improvement. The 
objective of the long tenn operational model is to meet 
energy requirement, over the planning horizon, at minirmni 
operation cost, i.e., to optimally, allocate load among 
available plants. The solution of the operational pla-ming 
model vill provide us the annua3. plant factors of the various 
pover plants in the region such that total discounted costs 

are minimum* 

The operational ilamiing model Is a non-linear 
prograaaing model. Ito decision enables are to choose 
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the yearly generation of energy (IMH) at each of the 3 
alternative plants in the region in each year over the 
planning horizon. 

Lefinition of Terminologies 


t — 1 j 2j • . • T “> 


i — 1 j 2 j ... •* 


3 — 1j 2j ••• Ng “ 



h P 

®13t 


> 

®i3t» ®i3t 


b 


E 





represents the time subscript which 
refers to the number of years , 
refers to type and amount of energy 
generated, imported or exported, . 
refers to location from v/hich energy 
is being generated or exported. 

Number of IIWH to be produced during t 
by i^^ (type and size) plant at 3 ^^^ loca- 
tion. This is a continuous variable, 
refers to base load energy and pealdng 
energy, 

the i^^ amount of energy imported (+ 
or exported (- (Tliese are cons- 

tants deteminod by the iiiiport-export 
capacities chosen in the capacity expan- 
sion problem,) 

refers to peak and base capacity. 

refers to net base peaking energy 
demand respectively in the 3 region 
during the t^ tiJie period. 
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%3t ^ 



m. - 4 . 

13 1 




fuel price for type of olant at the 
3* location during period t. 

The average not heat rate function of 
the i^^^ type of plant at 3 ^^ location 
during period t. This is a function 
of the generation historj" of the plant 



Variable operation and naintenanco costs 
for the i''^ type of plant at 3 ^^ location 
during period t. 

Penalty costs for adverse effects (such 
as pollution) of the i^^ type of plant, 
at the location during period t. 

The maxiraum total number of 1-iWH that the 
i^^^ type of plant at the 3 ^^ location 
is allowed to generate during the time 
period t. 


With the above terminologies, \/e specify the objective 
function as 


Minimise 



N. 


\jt ( Tiijt 5 »ijt • 


+ ^®ijt Pijt^ *^5^ i 


( 8 . 10 ) 



such that 




E 




i i 






( 0 . 11 ) 


for i = 1, ... 

3 = 1, ... .'I 2 

t = 1, ... T 

•^iot 4^idt fori = 1, ... 

3 = 1 , ... 112 ( 8 , 12 ) 

t = 1, ... T 

Equation (8.10) is the objective function which is the nini- 
misation of total discounted fuel costs, plant oper-^tion and 
maintenance costs plus penalty costs. 


In section (8,2) the general model of the integra- 
ted electric pov/er system ms presented using the ideas of 
Bender' s . In that formulation of the problem, the non- 
linear continuous function f (y) correspords to the opera- 
tional planning sub-problem here. The capacity expansion 
sub-problem described in section 8.3 is represented by the 
non-linear integer va.riable portion i.e., function C (x). 

The integrated Bender’s decomposition formulation is a 
non-linear mixed integer prograniBing probl^ idiere tlie 
complicating variables are identified with the integer 
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model. Farther into the future %;e project, tlie projections 
•become more and more uncertain. Hence adaptive modelling 
is recommended. Results of the model arc inrolemented for 
the iramectiate future (say five 3 ’‘ears for a planning ho risen 
of 'tliirty years) after which tlie model is modified or adapted 
to tahe into account changes that have taken place in the 
ii.TT.odiate and recent past. The problem is again solved and 
decision of the immediate future are implemented. Uncertai- 
nity of nuclear fuel costs calls for such an adaptive noddling. 

The primary component of the steam fuel costs are 
directly proportional to amount of energy generated and arc 
about ^ 0 % of the total costs of fuel. The scconda.ry costs of 
fuel are attributed to spinning reserves, starting and stopping 
of plants. The secondary fuel burnt is a function of the 
net heat rate and hours of use, and constitutes appro:Qj;iat,ely 
10^ of the total costs , 

The nuclear fuel costs comprises of cost components 
like value of the fuel burned, fabrication costs, shipping 
costs for both irradiated and new fuel, chemical reprocessing 
of irradiated fuel and such other costs. Pud turn up or 
irradiation level is given by MW days of heat generated per 
ton of uranium. 
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CHAPTER IX 


RCSULTS, DlSCuSCloriS, CONaUSlOJK AiiD SCOPE 
FOR FURTHER RESEARCH 

In tliio chaptsr ve shall prase nt the nuEierical 
resnlos obtained for a, case studj* by using the models 
developed in Chapter VIII for the capacity expansion 
planning sub-problei.i and operational planning sub-problem 
of an electric po^^rer system in a region. The sources of 
data used for this case study have been listed in 
Appendix A, 

The planning tine horizon for the capacity 
expansion sub-problem as well as the operational planning 
sub-problem has been assumed to be tliirty years (1971 - 
2000) . For the capacity ex33a,nsion sub-problem the planning 
horizon is divided into six time periods each of deration 
five years, i.e. 1971-76 corresponds to period 1, 1976 - Si 
corresponds to period - 2 and so on. For the operational 
planning sub-problem the planning horizon has been divided 
into thirty annual jilafis.. each of duration one year. The 
demand projections have been given for five year periods. 
For calculating the future requirements of energy in the ^ 
region a load fa.ctor of 70 ^ has been assumed* Fnecgy 
requirements for each period have been calcafl.at®d by *1*119 
following relationship. 
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Enerey reciuirernents = 0.7 x 8760 x project capacity dcnanc?.. 

Ease demand has teen obtained from the relations liips given 
below. 

Base demand (IIW) = 0.5x Baximain demand (IS/) 

Ease energ:^ demand = rase capacity x 876O (JIWII) 

Pea}: energy demand = Total energy de;..ajid - Fase energy 

demand . 

The life span for steam plants and hydel plants have been 
assiL'-ed to thirtyfive and forty years respectively. For 
calculating the total capacity demand for which planning 
has to be done, the reserve capacity is tcdien as fifteen 
percent of expected pea3: demand. 

t 

9.1 TI-IE MODEL FOR TiE SYSTEII USED 
FOR TEE CASE STUDY 

The entire region has been devided into turce 
zones A, B, C for tlie purposes of this study, baseci on the 
sirailarity of the characteristics of der.and in each zone. 

.Vs ejuilained earlier, toth the capacity expansion ani tho 
operational liannins models assicief. that the resion heing 
supplied hy tho eloctric..i>c«ervBystem under study has 
,^er, „ n ri centres betiffien vhlch long dista^xe transmission 

takes place. 
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The fixed cost of alternative 
projects not located at these 
load centres will be increased 
by the cost of transnission 
lines that v/ould have to be 
built fron the plant to supply 
electric! tj'- at load centres. 



This ijodel utilises the results of a screening; 
study nade by the plaiming agencies to cletemine the appro- 
priate alternative projects in the entire region. The rcsiats 
of this screening study supply the basic input data to the 
nodels proposed in this dissertation and hence the initial 
screening phase is vitally important. If the resulting 
expansion alternatives to be supplied to the ,.’.odel arc not 
good, the methodology suggested in this dissertation can oiii.:’ 
find the best ai^iong the alternative expansion policies which 
may not be optiraal. 


A selected sample of the results of such a screen- 
ing study have been given in Table - 23. for the entire 
region 97 alternativespowcr plants, ^ transmission line 
alternatives , and 44 importation - exportation alternatives 
have been chosen as feasible alternatives, to be supplied 
as input data to the model. 

The plant factors have been mainly used to calcu- 
late costs and do not represent initial operating policies* 
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table .23, 


A. sajnalj., .of alterrijat iro plants, transmissinn 

exportat i on alternatives f or 
oDtaincd from screenlnp; stud f ^ 

(i) Alto mat! VP Pngrv^ 


Zone 

of pro- 
5 ject (loca- 
5 tion) 

1 Plant typo 

2 and size 
\ 

i Probable i} S3Tnbol used 

1 installation 8 for this 
5 *D 6 riod 8 t.oT'naf-.i 

A 

Guyaiiiar - I 

100 MW 

Stcan Plant 

1971 - 75 

Xi 

A 

t 

Guyamas - I 

« 

25 MW 

Stean Plant 

1976 - 80 


A 

Baja 

30 MW 

Hydcl Plant 

1981 - 85, 
1986 - 90, 
1991 - 96 , 
1996 - 2000 


B 

• 

4 

1 

Obregon -I 

• 

• 

30 MW 

Gas Turbine 

• 

• 

1971 - 75 

• 


• 

• 

B 

• 

Obregon III 

• 

120 MW 

Gas Turbine 

• 

» 

1976 - 80, 
1981 - 85, 
1986 - 90, 
1991 - 95 j 
1996 - 2000 


C 

t 

» 

$ 

t 

Eacurata 

20 MW 

Hydel 

1971 - 75 , 
1976 - 80, 
1981 - 85 


c 

1 -iachis - IV 

200 MW 
Steam 

1986 - 90, 
1991 - 95 , 
1996 - 2000 

^ 5 ’ ^ 1 ’ 


(continued, on next pase) 



Z 'l t 


lABLL 23 (Continued) 

( ^ alter n atives frorii the screen ing 

s tudies 


Sjrmbol for 5 Import to 
each alter- 0 zone 


0 J^ount of 
0 import 


t Probable iS 
5 importation 0 Type 



Type 


Transmission line ; nro.iects between Zon es A and B 


3 - 5=^700 


795 ^1=^00* ^3=300 Z5=300 21^=300 ... z.,^=300 

^2=700 2i^=700 %=700 2^2=700 

T-rfl-nsmi .qsi nn line 'pro.iects bet\;ee n Zones E and_G 

795 2^3=300 3^5=300 233=300 

Z,,.=700 ... 


♦*S/C 220KV „ 

795 ACSR 2.J2-300 2^0-300 


S/ C Z \ ““700 ••• ••• ••• 

' •2168A0SR 

* Assraned. transmission capacity in MW 
S/C mean single circuit steel tovier 

+ ACSR refers to conductor size -i 

++ Z is the symbol (0 - I) variables used for tJaese alternatives. 
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Plant factors that liave been used for calculating variable 
costs are as follows : 

Hydel plants - 25 % 

Stean plants - 55 % 

Gas turbine plants - ^ 5 % 


The ca-pacity expansion r.Todc-1 for the region considered is 
as folio'll : 


Miniiiiisc 


where 



^ *1- U 

'it I ZT Z1 i ^t * \tl 

-J k=1 t=1 ^ 


\ 

j 



^t 

n^t 


= discounted cost of plant i constructed during 
period t (includes operation costs). 

= 0-1 decision variables for plant i built 

during period t. 

= 0-1 decision variables for transEiission lines 

k, built in period t. 

= present worth cost of transniission line k, 
constructed in period t. 


Demand constraints for Zone A 


U 


3 u U U D 

Rj. Ill + ^yO' ^11 ^ ’“l,1 




i=1 


i=1 


3 n 6 3c 0 * 

yit * Z z % "it " *1.0 - z 

t=1 i=1 


B 

It 


»®6 
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Demancl constraints fo-r? 7^r\f^ B 


6 

y 

1=1+ 


u ^ u u 


u 


u 


i=lK 


2,0 


^21 ^ 1,2 


^ U ^ ^ U 

Z 2 Z 

1=4 t=1 i=lf 

Demand Constraints for Zone C 


^ U U 

Z "2t^ 

.t=1 


6,2 


9 

ST 


U 


9 


u 


S . y. 
1 *^1 


11 


+ A 
i=7 


u 


30 "31 


u u 
4.1 > 4 


13 


9 


z 

1=7 


U 

^1 ^it 


6 


Z 

t=1 


where 


9 


Z 

1=7 


u u ^ u 

\ -^30 " Z ^3t 

t=i 


u 

^ ^6,3 


U 



u 


u 


^it 



base, and peal: 

various types of plant of type D retired. in 

zone n at beginning of tine period t 

generation capacity i^^ type of plant (of type U) 

import export capacity of alternative i (of ty-je U) 

0-1 decision variable for i^ toport export 
alternative in period t 

Installed generation capacity of type U in 
region '3 at tie staH of jaanning imrlzm 



caDacity deman:! of type u in region n 
during period j . 
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Transmission line co nst raints between Zone A and Zone B 
f — 4 „ 2 — . 




5^11 - Z ^ =^1-2 

u=base or :-1 i -1 — » 

peali 


6 2 


L j'l L 2 — - - 2_ Z *1 ht\^ "i-2 

usrbase or i t=1 i=1 — » 

peak 

Transmission line constraints between Zone B and Zone G 


^i1 " ^ ^i1 1 ^^2-3 


u=t>ase or 
peak 


4— U 6 4 -j 

^ nt- Z Z %t 4 ’^a-s 

u='base or i — : 

peak 


where 


t. = transnis S' '^n capacity of transmission line i 

T = ’bransmis sicn capacity Ije tween zone A and Zone B 

1-2 

at start of plannin" horizon. 

^ - transm'-ssion capacity hetween zone B and Zoi» 0 
at start of planning horizon 

L F. and L Fg ais tho loss factors of trawntosion linos 

tetween A and B and B anfl C wpootisrtj. 
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Acldltiong. . cons tralnt s for liml tatinn of the total 
instaile^capacity at the alternatilm hvdro sites 


6 

Z_ 

t=1 


^it 


4 : ■' 


for i = hydro alternatives 


Tills Tiiodel has got in total 199 variables and 5^ constraints 


9.2 OPERiTIOm PLAmMG MODET. 
FOR THE EEGM ’ 


To find the minimwa annual cost operation of the 
electricity system as it expands over the planning horizon- 
according to the schedule determined by the capacity expan- 
sion model, the t^dro planning sub-model, first calctilates 
the operation of hydro-electric plants that maximise their 
firm on-peak energy, ThJLs energy is subtracted from the 
projected annual energy demands. 

The operational planning model then calculates the 
annual energy generation (and thus the plant factor) of eadi 
pla.nt durir^ each period to get a minimum discounted cost 
of operation over the planning horizon. It is appropriate 
to emphasise that tlie purpose of the operational planning 
sub-problem is not to obtain a policy tliat will be recocroien- 
ded as the operating proceduia to te followed in fliture 
since only the expected variable (curating and :B^ntemiice 
costs plus fuel costs) costs per unit of energy are inclialed 
in the operational planning sub*'«)deil* 
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Having found the system expansion policies from 
the cap? city expansion sub-model and knowing the maxiiMur. 
firm on peak available energy from the hydro-electric plants 
in each zone, the operational planning model can be forrau- 

As observed from bhe general forriul8,tion of the 
operation planning model, the objective function consists 
of the sum of total discounted fuel costs plus plant opera- 
tion and maintenance costs. In particular the net heat rate 
function, II (number of BTU/ICWH) which determine the fuel 
costs a.re very difficult to identify without adequate data 
and more importantly mthout knovang future plant operating 
policy. Hence as per the data sunnlied, \]e a,ssumc 2,5 KWIv' 
unit of fuel burnt at beginning of policy. TJrLs malies the 
operational planning problem an uncoupled linear program, 
one for each year of the operational planning horizon. 

The variable cost and fixed cost associated \;ith 
each plants are calculated for a rate of discount of eight 
percent, Tliis data has been pronded by the planning autho- 
rities directly. 

The operational planning for tne „iiyst i1jergjM.p j2 
has been formulated as a series of linear parograniBing 
protlem^, one for each period (one year) of the planning 

horizon. 



The general form at 


xear t is given by 


Minimise 


such that 




^ r^ljt * ”'i; 


s: _ _ ^ 

1=1 3=1 


-3t 


"ijt 




J 




S' S' 


3 


e. .. 
I3t 


where 


p p ^ 3=1. -.Ho 

d. > E.. + y c 

ijt 2t - ^ ^ °ijt 

4 °13t 

1 = 1, ... N, 

0 ~ ) • • • 

^i3t = total fuel cost/\mit of energy 
t = 1 - generates the first year operational planning 
model 

t = 2 - generates the second year operational pla:ini‘i; 
model and so on. 


calcifLated from the equation 


^i;)t 


fuel costs ( = cost ner unit amo^t) 
(IvWH produced per unit amount) 


The operationaJ. planning studies for year 1971 has a coiW'* 
traint matrix of dimension 15 x 9* A snail programme was 
written to generate the cost »trioes, flae I.*P* pactoigeprosraiaca 
available at the Ccmputer Centre, I.I.T. Itejur has been used 
for sol vine this onerational plamtng pre^ae®* 



9*3 FQRimiATIOIl OF TIIC CAPACITY EXPAHSIOII 
MODEL TO BC USED FOR COMPUTATI Oil 


Prcm the data available it is possible to derive 
in a straight fonra.rd manner the integer programming model 
for the system. 

The objective function is formiolatcd as 

r 97 2N- . 

mnimise) 27 ^ ^ 

^ ^ 5=1 3 

where 

» 

x^, and are 0-1 variables, 

Cj^s and k^’s are the discounted total costs of insta- 

*L J 

nation, operation and maintenance of plants and trans- 
mission lines resjoectively. 

The satisfaction of base and peak demand constraints for 
Zone A at every period is given bye the following cons- 
traints im posed on tlie objective function: 

100 X 25 Xj + 25 \ * 50 X 5 + 100 1 (Base demand) 

45 X 3 x15\+ 5 I 1 +IOY 2 + 15 Y 3 •;^1 2 (Pea!: demand) 

m 

100 1- 25 Xg + 25 X^ + 100 + 25 X^o ^ ^15 

150 X25 + 300 X27 + 100 x^ -»• 25 + 50 \2 

+ 100 ^3 + 150 ^ ^29 CBise) 
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30 + 60 \ + 30 + 6o Xg + 90 X^ * 60 

+ 11 X.jg + 7 X-jr^ + ^ .j2o + 11 X22 + 7 X23 

+ 90 X2^ + 120 x^^ + 7 + 180 Xq ^ + 270 :%2 + -“^o 

+ 7 Xg^ + 30 -:. 5 q + ^20 Yi^q ^ 640 

(Pcc.k) 

Similar Gnuaiions can "be formulateo. for other 2c ■ *2. 


9 .If PIS gissi 0" _0F i-mERi qai . kesitlts obtaine: 


r'QR TPIE SELECIF'D G^r: STUDY 

The alternatives chosen for plants, transmission 
3.ines , and exportation - importation ty the capacity expan- 
sion 37odei at first iV- ration are shown in Tables 2 k, 2 ^ 25 . 

The resell ts for the first iteration of the opera- 
tional planning model is shown in Figures 9,1 to 9 » 3 « 
thirty year estimated, and calculated plant factors for chosen 
plants arc given in Table 27. 


The first obs'jrv. that can be made fron i’lG. 9«1 
is th.at d^ (energy from the first unit added to the sy.' torn) 
is zero until the end of 1971 - 75. Similarly \ (energy 
from first unit added to tl-® Cfcregon region, and fourth 
in sequence of insteO-iation) is not needed iintil the second 
year of period 1976 - 8g, and dg = 0 till tlw teginning 
of 1976 - 80 period. 
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TAIIiE ^ 

Secmences of installations of power pi ante . 
( Results fmn _caT)acitv expansion sub~pr~oiQ.em ) 


ZQHE 

5 Project chosen 
5(or existing) 

J J2Z1 

T Initial variable H Decision 5 Year 

5 0 & li and fuel 5 variables 5 of ins*- 
5 costs* IBLlls/KWHO 5 tallation 

i 

0 2 


1. G 0 

4 

(1 5 

A 

Guyamas 

96MW(S) 

]i-»2 

®1 


B 

Obregon 

2GhW(G) 

10.4 

®2 

- 

C 

Hochis 

if1^!VKS) 

4.1 

^3 

- 

C 

Culiacan 

l4lftKG) 

10.3 

% 

- 


A 

Guyanas 

100NVKS) 

4.2 

<^1 

1971 

A 

Hernosille 15MVJ(G) 

10.6 

d2 

1971 

B 

Nova j oa 

501W(S ) 

4.4 

^3 

1971 

B 

Obregon 

901'M( G) 

10.5 


1971 

C 

Mochis 

lOOlM(S) 

4,1 


1971 

A 

Hernosille 901^KG) 

11.1 

4 

1976 

B 

Nova j oa 

50MW(S) 

4.4 


1976 

B 

Obregon 

60KW(G) 

10.5 


1970 

C 

I’bchis 

lOOl'M(S) 

4.1 


1976 

C 

Hiites 

lOOMW(H) 

- 


1976 

A 

Caborca 

50MW(S) 

5.0 

^10 

I9GI 

A 

Hernosille 601'R'KG) 

11.1 

dll 

1901 

B 

Nova j oa 

50I&J(S) 

4.5 

*^12 

1931 

B 

Obregon 

120l'W(G) 

10.7 

^^13 

1981 

e 

Mochis 

lOOIM(S) 

4.1 

*^14 

1981 

c 

Culiacan 

60]iM(G) 

10,5 

^15 

1981 

c 

Culiacan 120I'&f(G) 

10.6 

*^16 

1981 

c 

Iioj)e2mateos79C2'®!fOE) 

- 

- 

1981 

A 

Caber ca 

751'M{S) 

5.0 


1986 

A 

He rnos ille 1 20MW ( G) 

11.2 

^18 

19^ 

B 

Novajoa 

50^M{S) 

4.1 

<taba.e cdistliiied) 

1986 
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TAEU2 2?+ (Continued) 


1 

5 p 

J ^ 

e 

3 

~T'”“ 

5 

4 

-j— -- - 

5 L... 

B 

Ilovajoa 

iooiH'/(s) 

4.5 


^20 

irc6 

C 

Ibcliis 

5001 i-KS) 

5.2 


^^21 

1906 

C 

Culiacan 

IfOOlH(G) 

10.7 


^^22 

19 u 6 

A 

Caborca 

300L¥(S) 

5.4 


^^23 

1991 

A 

Hernosille 

1801a^’(G) 

11.2 



1991 

B 

Novajoa 

i 5 oni'J(s) 

4.5 


d 25 

1991 

B 

Obregon 

i20ir'KG) 

10.6 


*^26 

1991 

C 

l’fc)ClTiS 

500Ii¥(S) 

5.2 


^27 

1991 

C 

Culiacan 

4001 M(G) 

10.7 


d28 

1991 

A 

Caborca 

3001IW(S) 

5.4 


*^29 

1996 

A 

Jlernosille 

240M(G) 

11.2 


<^30 

1996 

B 

Novajoa 

300 iiw(s) 

4.5 


^31 

1996 

B 

Obregon 

50LV(G) 

10.5 


^32 

1996 

B 

Obregon 

i^Ly(G) 

10.6 


*^33 

1996 

C 

I'bcliis 

2001IV/(S) 

4.5 


^34 

1996 

C 

Culiacan 

400I‘iy(G) 

10.7 


*^35 

1996 
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TABLE 25 

Import » Export Alternative 
Chos en a t second iteration 


Import to 
2tone 

~ll”!I^e’ of 

5 capacity 


Amount 

M.W. 

. IJPeriod during \;}Ticl: 
Straps fer occurs 

A 

Base 


25 

71 - 75 

A 

Peak 


15 

71 - 75 

A 

Base 


50 

96 - 2000 

A 

Peal^ 


90 

96 - 2000 

B 

Base 


^5 

71 - 75 

B 

Peak 


60 

:71 - 75 

E 

Base 


120 

71 - 75 

B 

Peak 


90 

96 - 2000 

B 

Peak 


60 

76 - 80 
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TABLE 26 


Transmission line alternatives 
Qiosen at second iterat ion 


Zones 

connected 

— IT 
0 
? 

Type 

Maxinun 5 

0 transmission 0 

5 canacitr (Ifl-0 5 

Year 

to 

start 

B - 

A 

220 KV 
S.C.R 

300 

1995 

C - 

B 

220 KV 
S.C.R 

300 

1971 
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TABLE 27 

Elrst iteration opcmtl on planninr: sub-problem re sults 


Decision f Estimated 5 Calculateci Decision 
variabledSaverage 0 average 5 variable 

5 plant 5 plant ^ 5 

S factor % 1} factor 5io Q 


li Estimated OCalculatecl 
0 average feverage 
5 plant folant 
j factor % gactor % 



55 

82 

dl7 

55 

66 

®2 

15 

?3 

'^IS 

15 

2?+ 

®3 

55 

85 

^19 

55 

85 

% 

15 

22 

^20 

55 

64 

'll 

55 

61 

d2i 

55 

68 

'^2 

15 

15 

*22 

15 

22 

'^3 

55 

81 

d, , 

u23 

55 

64 


15 

9 

d^ 

15 

17 


55 

72 

*^25 

15 

12 

■^6 

15 

17 

^26 

15 

20 


55 

73 

*^27 

15 

36 


15 

8 

<^28 

15 

26 


55 

63 

^29 

55 

16 

‘*10 

55 

60 

^30 

15 

27 

*11 

15 

25 

^31 . 

55 

60 

■*12 

55 

22 

d32 

15 

27 

*13 

15 

19 

^33 

15 

7 

*llf 

55 

63 

d^ 

55 

85 

*15 

15 

28 

^35 

15 

18 

d.,A 

15 

26 
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Since this represents extra unused capacity, in the next 
capacity expansion problen, we can delay these projects 
until 1976 - 80, time, period. Usin.?; the same reasoning 
constniction of the plant producing dg can be delayed until 
1981 - o 5 i and construction of plants producing d^^^ and 
until 1986 - 90. Prom figures 9,1 to 9.3 we get infor- 
mation regarding tlie time period of construction and can uti 

lise them for the next iteration of the capacity expansion 
prlblem. Similar interpretation can be giinsn for other 
values , 

The second iteration can nov; proceed for the capa- 
city expansion problen ^ath values of plant factors obtained 
from the first iteration of the operational planning prcblei’ 
Tiiis iteratiyo process can be continued till no further 
improvement in capacity installation sequences occur or the 
sequence converges. The total cost of installation at the 
first iteration with assumed values of parameters turns out 
to be ^896 million dollars. The second iterational capacity 
expansion totrl. cost is to the tune of ^225 million dollars. 
The second iteration uses different plant factor values 
(those obtained from first iteration operational planning 
model) and hence this difference in total costs are 


obtained 
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9.5 smiAEi A!J) conausioiB 

The purpose of this study uas to construct a 
niathcMatical model for uetemininp the best locations j tinej 
and size of future pover plants in a region. The economic 
objective vra.s the minimisation of total discounted cost, 
given a projected de-mand for electricity. 

The problem of planning the best capacity expan- 
sion of an electrical systen has been trsated in t’le litera- 
ture using qualitatiye approaches, plant by plant analysis, 
simulation and optimisation technj.ques, applied strictl 3 ’- to 
strategic (capacity expansion) aspect of the problem. 
Variable operating costs show to be a significant part of 
the total costs. Hence careful analjrsis of the operational, 
aspects of the problem is essential, if a truely minimum 
total cost capacity expansion policy is desired. A methodo- 
logy for tackl.ing both the capacity expansion problem and 
operational planning problem, in an iterative vay using 
adaptive modelling ideas at each iteration has been proposed. 
Considerable improvements in costs can be obtained by using 
the iterative methodology. 

The problem of dimensionality limits the sigie 
of the problem that can be solved by this methodcaogy. 

The integer code for capacity expaission sub-p3ro'KLei3 is 
based on Bala‘s algorithm Due to the tlae oost 

limitations , a trade off between of sdutioii 



227 


a.nc ccfmputation costs are necessary. However, tlic opera- 
tional planning L,P, problem is a systematic search procedure 
vrliich guarantees optimal solution. 

Ilodels arc simplificatioriS of reality and the one 
presented in tliis study is no exception. Howver, the mrj-n 
aspects of the electrical power system capacity expansion 
problem are included. Assumptions sucia as the deterministic 
character of the model v^ere made with the ulttoate goal of 
compromising optimality of solution \dth realism of the 
problem that can be solved. As more efficient algorithms 
are developed, bettor and more realistic prohLems can be 
solved -^rith the methodology presented in this dissertation. 

Planning studies such as the one presented l^ro 
are usually complicated hy the need to accurately project 
possible future economic conditions of the area. Th-is is 
a difficult task (if not iiiipossible to do) because of the 
many unexpected events tliat usually occur in a region. This 
problem can horover be circumvented by a repeated solution 
of the problem T > AT years (such as five years) and im- 
plement the decision for the first A? yaars of th© expan- 
sion policy obtained each time , This adaptive mof^iUlng 
process therefore makes use of both a planning hoxixon 
of length T an implementation period of A? 4. 
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9 .6 SCOPS FOR PCJRTIER STUDY 

The follovo.ng topics are recor.imended as fruitful 

areas for further study: 

1. Testing the methodology presented in tliis dissertation 
using the mixed integer capacity expansion model, uhere 
the continuous variables would represent the export, 
import decisions . This requires the development/or 
implementation of an efficient algorithm capable of 
solving realistic problems. 

2. Further investigation and numerical comparison of 
Bender's decomposition approach as modified by Geoff rion 
(■^8^) with the methodology suggested in this dissertation 
for larger problems appear specially fruitful, 

3. Development of stochastic decision making models for 
further investigation of the reliability aspects of the 
capacity expansion problem, 

4. Investigation of the dual of the operational planning 
problem and its role in the adaptive procedure used in 
this dissertation, 

. Development of models and systematic procedures or 
algorithms that combined with experience cc«ld be used 
for screening out altenKitives , This area is of 
special importance as the methoddogy of this disserta- 
tion assumes that a good screening Study bss been Bade 
before the capacity expansion problOB is solved • 


5 
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6. Development of models for finding the optiimim conouctive 
operation of the hydroelectric and other plants in each 
zone, 

7, Development of suo-models to the capacity expansion 
protlem that might help the plamxrs to systeriatically 
quantify some of the in’cangible aspects of the problems 
lil:e pollution and environr.ient. 
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